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Ici  rdgnait  aussi  cet  amour  qui  ne 
s'exprimc  plus  parce  qu'll  ne  participe 
pas  a  la  vie  de  ce  monde,  II  ne 
supporterait  peut-etre  aucune  4preuve, 

11  semble  a  chaque  Instant  trahl,  et 
la  moindre  amiti£  ordinaire  a  l'air 
de  la  vaincre,  et  cependant  sa  vie  est 
plus  profonde  que  nous-memes  et 
peut-etre  ne  nous  semble-t-il 
indifferent  que  parce  qu'il  se  salt 
reserve  pour  des  temps  plus  longs 
et  plus  surs, 

M.  Maeterlinck 

"Les  Avert! s"  du  "Tr£sor  des  Humbles" 


PREFACE 

The  nine  chapters  which  follow  represent  the  set  of  lectures  given  as  a 
final  year  one  semester  course  at  L'Universit£  Catholique  de  Louvain  for 
the  first  semester  of  the  1967-68  school  year.  Because  of  the  presence 
of  two  national  languages  with  the  lectures  given  in  a  third  it  was 
decided  to  record  the  material  as  covered  for  student  assistance  and 
availability  for  future  studies.  Also  the  material  often  records  in 
a  consistent  whole  unavailable  research  results,  and  puts  on  further 
record  the  nature  of  Joint  cooperation  between  our  associated  research 
groups  at  Stanford  and  Louvain. 

In  the  field  of  electrical  engineering  the  theory  of  state-variables 
has  raised  some  rather  paradoxical  situations.  On  the  one  hand  it  is 
often  claimed  that  nothing  can  be  achieved  with  state- variables  that  can 
not  be  done  with  more  classical  methods.  This  point  is  most  frequently 
raised  by  those  who  wish  to  construct  working  circuits.  On  the  other 
hand  the  mathematically  inclined  have  a  tendency  to  develop  rather  minute 
points  or  to  get  involved  in  the  elegence  of  the  theory  with  an  attendant 
sacrifice  of  the  practically  important  aspects.  As  a  consequence  the 
two  natures  of  theory  and  practice  tend  to  become  further  separated  when 
state-variables  are  involved.  Here  we  would  at  least  make  an  attempt 
to  resolve  this  paradoxical  situation;  that  is,  we  would  try  to  bring 
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theory  closer  to  practice  and  vice  versa.  This  is  done  by  presenting  a 
coherent  whole  with  emphasis  upon  thos_  aspects  of  the  theory  for  which 
use  can  almost  immediately  be  seen  or  which  have  proven  themselves  in 
practice.  Actually  the  subject  was  suggested  by  the  University  as  we 
felt  that  some  value  could  result  in  the  intended  types  of  treatment  we 
have  enjoyed  the  challenge  and  hope  that  the  venture  has  proven  profitable 
for  all  concerned. 

It  should  be  remembered  that  the  material  represents  lectures  and 
not  a  polished  book,  even  though  it  has  somewhat  the  form  of  a  book  for 
convenience  of  the  user.  As  a  consequence  of  its  lecture  form  as  well 
as  the  circumstances  of  its  construction,  there  is  much  omitted  which 
could  profitably  be  contained.  For  example,  there  are  points  of 
derivations  which  could  profitably  be  put  into  notes  for  completeness 
but  which  have  been  omitted  in  order  to  cover  the  material  desired  in 
the  allotted  time.  Of  equal  importance  is  the  scarcity  of  references; 
generally  only  a  single  reference  available  to  the  author's  students 
at  the  time  is  given  while  multiple  referencing  would  be  much  preferable. 
Likewise  there  are  some  topics,  as  topological  and  nonlinear  synthesis, 
which  have  been  almost  entirely  omitted  but  which  should  properly  not  be 
for  completeness.  Among  works  which  we  would  have  liked  to  add,  perhaps 
to  be  saved  for  a  revised  edition,  are  those  of  J.  Hiller  (active  theory), 
P.  Wang  (infinite  dimensional  theory),  H.  Watanabe  (nonlinear  theory), 

R.  Yarlagada  (topological  synthesis),  and  D.  Youla  (lumped-distributed 
synthesis).  A  list  of  symbols  and  an  index  is  appended  for  convenience. 

In  conjunction  with  our  belief  that  life  should  be  constructive 
and  associated  with  a  masculine  spirit  of  verse  which  enhances  its 
poetry,  we  incorporate  some  nontrivial  concepts  of  the  Flemish  writer 
in  French,  M.  Maeterlinck. 


R.  Newcomb 

Louvain,  January  1968 
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Om  dat  die  leeke  van  alien  zaken 
Rime  ende  dichte  willen  maken 
Gheijc  clerken,  dat  wonder  es, 

So  hebbic  mi  bewonden  des 
Dat  ic  nu  wil  bringhen  voort 
Wat  enen  diohter  toe  behoort, 

Die  te  rechte  sal  dichten  wel; 

Want  dichten  en  is  gheen  spel. 

Men  sal  ooc  voren  versinnen, 

Hoemen  dat  dicht  zal  beghinnen, 

Middelen  ende  daer  toe  enden. 

Jan  Boendale 

”De  Leke  Spieghel,  III" 
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In  the  words  of  M.  Maeterlinck  f'Les  Avertis"  du  "Tresor  des  Humbles" 

L’on  sent  que  c'est  l’heure  enfln 
d'affirmer  une  chose  plus  grave,  plus  humalne, 
plus  reelle  et  plus  profonde  que  l'amiti6, 
la  pitie  ou  1* amour;  une  chose  qul  bat 
mortellement  de  l’aile  tout  au  fond  de  la 
gorge,  et  qu'on  Ignore,  et  qu'on  n'a 
Jamais  dlte,  et  qu’ll  n'est  plus  possible 
de  dire,  car  tant  de  vies  se  passent 
a  se  taire!  ...  Et  le  temps  presse. 

for 

M.  A.  Gillett 
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Souvent,  nous  n'avons  pas  le  temps  de 
les  apercevoir;  ils  s'en  vont  sans  rien 
dire  et  ceux-la  nous  demeurent  a  jamais 
inconnus . 

M.  Maeterlinck 

"Les  Avertis"  du  "Tresor  des  Humbles" 


CHAPTER  I 

INTRODUCTION  -  THE  STATE 


A.  Summary 

Here  we  briefly  review  the  philosophical  nature  of  the  state  giving 
a  more  or  less  precise  mathematical  formulation  in  terms  of  system 
transformations  and  network  relationships.  An  example  concerning  the 
Brunt  structure  is  given  to  illustrate  various  points  of  the  theory  to 
be  considered. 

B.  The  State  -  Intuitively 

Intuitively  an  object  can  be  described  at  a  given  instant  by  a 
certain  set  of  conditions  which  in  fact  are  specified  by  the  object  be¬ 
ing  described;  these  conditions  are  often  referred  to  as  the  state  of 
the  object.  However,  in  scientific  discussions  the  state  is  usually 
taken  to  mean  that  set  of  conditions  which  when  specified  at  a  given 
instant  of  initiation  of  an  excitation  lead  to  a  predicted  response  over 
the  period  of  excitation.  Thus  the  concept  is  generally  applied  to 
causal  (that  is,  nonanticipatory  or  equivalently  antecedal)  systems  where 
it  is  possible  to  predict  the  output  to  a  given  input.  A  specification 
of  the  necessary  conditions  to  allow  determination  of  the  output,  that 
is  an  assignment  of  initial  conditions,  is  essentially  a  specification 
of  the  state.  The  state  then  is  that  entity,  described  through  a  set 
of  parameters  (perhaps  uncountably  infinite  in  number),  which  when  pre¬ 
scribed  initially  allows  a  unique  motion  of  the  entity  under  determinate 
excitations.  We  shall  soon  make  the  concept  precise  mathematically  at 
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which  time  we  will  see  that  a  slight  modification  is  of  interest  for 
treating  networks. 

C.  The  State  -  Uses 

Although  of  the  most  recent  development,  our  primary  interest  will 
be  the  use  of  the  state  for  design  or  synthesis.  For  synthesis  we  need 
to  develop  a  formulation  which  is  convenient  for  decomposition  and 
construction.  In  obtaining  a  suitable  development  we  shall  investigate 
analysis  methods  from  which  we  will  see  that  by  isolating  a  set  of 
state  variables  a  convenient  analysis  method  is  obtained.  The  method 
is  especially  convenient  for  digital  computer  formulation,  and  thus,  we 
will  obtain  several  methods  for  digital  computer  analysis  for  circuits. 
The  results  are  further  useful  for  investigation  of  the  transient  and 
frequency  responses  of  networks  as  well  as  for  the  determination  of 
natural  frequencies.  Similarly  a  useful  technique  for  investigating 
sensitivity  is  btained.-  Of  particular  Importance  is  also  the  means  of 
determining  "all"  possible  equivalents.  By  reversing  the  analysis 
process  one  is  led  to  several  design  formulations.  For  example,  given 
a  transfer  function  one  can  algebraically  set  up  a  canonical  set  of 
state  variable  equations,  by  a  means  suitable  for  digital  computer 
programming.  From  the  canonical  equations  one  can  revert  to  an  analog 
computer  realization,  the  result  being  of  considerable  use  for  integrated 
circuit  design  using  operational  amplifiers.  By  another  interpretation 
of  the  canonical  equations  one  can  obtain  an  alternate  minimal  capacitor 
synthesis  by  loading  a  gyrator- resistor  network.  By  proper  generaliza¬ 
tion  of  multivariable  functions  we  can  also  develop  a  synthesis  for 
lumped-distributed  circuits. 

Although  it  can  be  claimed  that  the  state  variables  are  nothing  more 
than  an  appropriate  choice  of  variables  for  ini  vial  conditions,  such  an 
outlook  is  rather  narrow.  In  fact  previous  results  obtained  from  an 
"initial  condition"  outlook  are  rather  weak  and  shallow  when  compared 
to  what  has  been  achieved  by  the  state  variable  outlook.  From  the 
previous  paragraph  we  can  summarize  the  results  of  state  variable  theory 
to  be  discussed  in  the  sequel  by  the  following  topics; 
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1.  Digital  Computer  Analysis 

a.  Formulation  of  canonical  equations 

(1)  Topological  means 

(2)  Reactive  extractions 

b.  Transient  analysis 

c.  Frequency  response 

2.  Analog  Simulation 

a.  Integrated  circuits 

b.  Filter  design 

3.  Equivalence 

a.  Minimal  realization  transformations 

b.  Nonminimai  (encirclements) 

4.  Sensitivity 

5.  Finite  Synthesis 

a.  Minimal  realizations 

b.  Loaded  n-port  theory 

c.  Lossless  synthesis  (hybrid) 

6.  Multivariable  Realizations 

a.  Minimal  realizations,  etc. 

b.  Lumped-distributed  synthesis 

c.  Noncommensurate  line  synthesis 

7.  Distributional  Generalizations 

a.  Representations 

b.  Time-variable  circuits 

8.  Inif inite-Dimensional  Extensions 

D.  The  State  -  Mathematical 

Let  us  consider  as  given  a  system  designed  to  map  inputs  u^  into 
outputs  If  we  know  all  inputs  applied  to  the  system  from  its  time 

of  construction  to  the  time  of  observation,  t,  then  y(t)  is  "uniquely" 
known  and  is  determined  through  a  knowledge  of  the  system  transformation. 
However,  it  is  more  frequent  that,  we  have  on  hand  a  given  system  which  we 
will  begin  to  use  at  time  t  ,  generally  without  a  knowledge  of  the  in¬ 
puts  applied  before  tQ.  We  will  assume  that  there  is  a  set  of  parameters 
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s(t  )  which  we  can  measure,  or  somehow  determine,  such  that  if  the 

w*  o 

input  u(t)  is  known  for  t  >  t  then  also  for  t  >  t  the  output 
— - —  —  o  —  o  - - — 

y(t)  is  uniquely  determined  [upon  a  specification  of  the  state  s(t  )]. 
«M  -  '*•  o 

Since  the  output  is  uniquely  determined,  there  exists  a  transformation 

T[  • ,  •  1  such  that 

*r 


y  =  T[u,s(t  )],  t  >  t  (1-1) 

Since  tQ  can  vary,  the  state  £  is  also  a  "function"  of  time  as  is 
of  course  reasonable  on  intuitive  grounds.  We  point  out  that  in  general 
y,  u,  and  particularly  s  are  multidimensional  quantities;  we  will 

m  *■'  *** 

take  u  as  on  m-vector,  y  as  an  n- vector,  and  s  as  a  k- vector  [for 
example,  k  will  often  be  the  number  of  capacitors  and  inductors  in  a 
circuit],  Pictorially  Eq.  ( I— 1 )  is  represented  as  in  Fig.  1-1. 


input 


transformation 
(with  state  s) 


O— > 


t  >  t 


Fig.  1-1.  SYSTEMS  REPRESENTATION. 


A  system  which  can  be  represented  by  a  transformation  of  the  form 
of  Eq.  ( I— 1 )  is  conveniently  called  a  Btate  determined  system.  One  can 
in  fact  make  a  detailed  study  of  the  general  types  of  state  determined 
systems  [1,  p.  67]  but  it  seems  more  important  for  our  purposes  to  proceed 
to  other  studies.  However,  we  define  a  few  useful  concepts.  First  is 
that  of  the  zero  state  9,  defined  through 

—  .  m 

0  =  T[O,0(t  )],  t  >  t  (1-2) 

o  —  o 

In  other  words  a  zero  state  is  any  state  which  gives  a  zero  output  for 
a  zero  input.  As  an  example  of  a  nonzero  zero  state  consider  the  balanced 
bridge  circuit  of  Fig.  1-2  where  the  capacitor  voltage  serves  as  the 
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state,  Si(t)  =  (v^t)},  and  we  take  the  applied  voltage  as  input  with 
the  source  current  as  output.  When  the  applied  voltage  is  zero  no  in¬ 
put  current  flows  as  is  seen  by  the  redrawing  shown  .'n  the  (b)  portion 
of  the  figure;  thus,  0  =  (v  (t)}.  We  observe  that  in  this  system  all 


Fig.  1-2.  NONZERO  ZEiiO-STATE  EXAMPLE. 


states  are  the  zero  state,  but  in  general  such  will  not  be  he  case. 

For  example  if  we  had  taken  i^  as  the  output,  the  output  would  only 

have  been  zero  if  V£  =  0,  that  is  for  this  new  system,  with  u  =  v, 

y  =  i^  the  state  s  =  (vc(t)}  is  only  tha  zero  state  when  it  is  zero; 

0  =  f0}. 

«•* 

With  the  concept  of  the  zero  state  on  hand  we  can  consider  the 
definition  of  a  linear  system.  A  system  is  called  linear  (with  respect 
to  inputs)  if  for  all  constants  k,  all  initial  states  s(t  ),  all 

haa  O 

zero  states  0(t  ),  and  all  inputs  u  and  u 

l[k<srs2,'2(V1  -  (I-3) 

We  observe  that  because  of  the  need  to  consider  the  state  there  is  a 

difference  between  a  linear  system  (in  its  mathematical  representation) 

and  a  linear  transformation.  An  immediate  consequence  of  this  definition 

of  linearity  is  the  fundamental  decomposition  obtained  by  taking  k  =  1, 

u,  =  u,  u„  =  0. 

»<1  ’ 
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T[u,srt)]  =  T[0,s(t  )]  +  T[u,0(t  )]  (1-4) 

That  is,  for  a  linear  system  the  total  response  can  be  broken  into  the 
sum  of  two  parts,  one  of  which  is  the  zero  input  response  and  the  other 
of  which  is  the  zero  state  response.  Thus,  superposition  not  only  holds 
with  respect  to  inputs,  as  Eq.  (1-3)  shows,  but  also  with  respect  to 
tne  response  from  initial  conditions. 


v 

c 


Fig.  1-3.  INTEGRATOR. 


As  an  example  of  the  decomposition  let  us  consider  the  integrator 
of  Fig.  1-3.  The  describing  equations  can  be  taken  as 


Ri 


i 


Cd[v  /K  -  v  J 
o _ o 

dt 


which  upon  simple  substitution  of  the  first  into  the  second  yields  the 
following  differential  equation  completely  in  terms  of  input  and  output 
variables. 


dv  v 

IlC(l-K)  o  o 

K  dt  +  K  "  Vi 


(I-5a) 


To  obtain  the  transformation  napping  the  input  into  the  output  this 
differential  equation  must  be  solved.  We  find  by  any  of  several  means 
(Laplace  transforms,  for  example) 
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We  see  that  Eq.  (1-4)  is  satisfied  and  that  s(t)  =  (v  (t)}  =  {y ( t) } 

"  1-K  ° 

is  a  suitable  choice  for  the  state.  Since  v  =  — — -  v  we  also  see 

c  K  o 

that  an  appropriate  (alternate)  choice  for  the  state  is  s^t)  =  { vc ( t) } . 

Perhaps  much  more  should  be  said  about  the  domains  of  definition 
of  the  various  quantities  but  such  discussions  can  also  get  lengthy. 

We  merely  mention  that  for  a  given  system  there  is  usually  some  restric¬ 
tion  on  the  type  of  inputs  allowed  as  well  as  the  range  of  outputs  for 
which  the  mathematical  transformation  T[ • , - ]  is  valid.  In  our  study 
we  will  most  often  assume  that  the  input  and  output  are  zero  before 
t  =  tQ  and  that  they,  along  v'ith  the  state,  are  real  valued. 

For  linear  systems  it  will  often  be  possible  to  find  a  description 
in  che  form 

A(t) s(t)  +  B(t)u(t)  (I-6a) 

C(t)s(t)  +  D(t)u(t)  +  E(t)  (I-6b) 

•M  **  »|  dt 

If  such  can  be  found,  these  equations  are  called  a  canonical  representa¬ 
tion  and  the  set 

R  =  CA,B,C,D,E} 

is  called  a  realization.  For  such  a  system  having  the  dimension,  k, 
of  the  state  finite,  we  ascribe  the  name  finite  or  differential  system. 
Likewise,  if  the  coefficient  matrices,  A^(t) ,  ...,  are  constant  then 
the  system  is  called  time- invariant  (actually  this  time-invariance  is  a 
special  case  of  a  more  general  definition  applicable  to  any  state 
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determined  system  [2,  p.  J).  In  most  situations  of  interest  the  useful 
information  about  the  system  is  contained  in  the  matrices  A,  B,  and 
C,  so  we  will  often  assume  that  either  E=Q  or  D=E=0.  Thus, 

**  •**  «*t  m  rf* 

most  of  our  concern  will  be  with  the  canonical  set  of  equations 


ds 

~  =  A  s  +  B  u  (l-7a ) 

y  =  C  s  +  D  u  (l-7b) 

mi  0*  m  m 

and  the  realization 

R  =  U-7c) 

It  is  possible  to  interrelate  the  canonical  equations  with  the  zero 
state  response,  Tju, ^(tQ ) ] ,  in  the  time-invariant  case  (a  similar 
development  holds  for  time-varying  systems).  When  the  realization  R 
is  constant,  Eqs .  ( 1—6 )  yield  a  continuous  transformation,  in  the  sense 
of  distribution  theory,  mapping  inputs  into  outputs  (in  the  zero  state). 
Consequently,  there  exists  a  matrix  Ji(t)  such  that  [3,  p.  J23] 

Ttu^Uo)]  =  h  *  u 


=  /  dT 


(1-8) 


where  *  denotes  convolution,  that  is,  the  integration  exhibited  (recall 
that  u(t)  is  zero  for  T  <  t  ).  The  n  X  m  matrix  h  consists  of 

Q  WH 

distributions  (functions,  impulses,  etc.)  and  is  called  a  distributional 
kernel ;  physically  it  represents  a  matrix  of  impulse  responses.  For 
Fig.  1-2  we  have,  for  example, 


( I-9a ) 
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while  for  Fig.  1-3  we  have 


h(t) 


K 

RC(l-K) 


RC(1-K}]| 


l(t) 


(l-9b) 


where  1 ( t )  is  the  unit  step  function  and  6(t)  =  dl(t)/dt  is  the  unit 
impulse.  By  taking  Laplace  transforms,  denoted  by  £[  ],  we  have  from 
Eq.  (1-8) 


=  H(P)  £[u]  ( I-10a) 

H(p)  =  f[h]  (1-lOb) 

where  JH(p ) ,  p  =  a  +  ja),  is  called  the  transfer  function  matrix  (it  is 
n  X  m  also).  By  taking  Laplace  transforms  in  Eq.  (l-6)  we  can  obtain, 
by  straightforward  substitution,  an  alternate  expression  for  the  transfer 
function  matrix 


H(p)  =  PE  +  D  +  C(plk  -  A)"1  §  (i-iOc) 

where  1.  is  the  k  X  k  identity  matrix.  One  of  the  problems  of  the 
theory  is  then  to  find  a  realization  R  =  (A,g,C,D,E)  given  a  transfer 
function  H(p)  since  then  the  canonical  equations  are  on  hand.  A  similar 
problem  is  to  obtain  the  canonical  equations  from  a  given  physical  struc¬ 
ture.  We  comment  that  Eq.  (i-lOc)  shows  that  the  transfer  functions  re¬ 
sulting  from  the  canonical  equations  are  always  rational,  when  k  is 
finite,  and  possess  at  most  a  simple  pole  at  infinity;  in  the  more  commonly 
treated  case  where  E  =  (),  H(p)  „ias  no  pole  at  infinity. 

We  will  illustrate  some  of  the  above  points,  while  exhibiting  a  set 
of  canonical  equations,  in  the  following  example  of  a  Brune  section.  How¬ 
ever,  first  we  comment  that  we  have  considered  a  given  construct  as  a  sys¬ 
tem  by  "orienting"  its  variables,  that  is,  by  specifying  inputs  end  outputs. 
Thus,  as  we  already  saw  in  Fig.  1-2,  a  given  construct  can  yield  several 
different  systems  by  having  different  inputs  and  outputs  assigned.  Never¬ 
theless,  the  state  will  generally  remain  invariant;  that  is,  given  a  con¬ 
struct,  there  is  an  associated  state  which  in  fact  can  be  used  with  all 
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svsteras  obtained  Irom  i he  construct.  Further,  a  network  has  been  de¬ 
fined  by  the  set  of  all  pairs  K.iJ  of  voltages  j/  and  currents 
allowed  at  its  ports  (4,  p.  7].  We  could  proceed  from  this  definition 
of  a  network  to  introduce  the  state  as  a  set  of  parameters  needed  at 
time  t  to  specify  allowed  pairs  [v,ij  for  t  >  t  .  But  for  our 

O  Mm  '  *«W  Q 

purposes  it  is  sufficient  to  orient  variables  at  the  network  ports  and 
work  with  inputs  and  outputs,  as  for  example  through  the  admittance  or 
scattering  matrices.  We  note,  though,  that  in  any  characterization 
there  is  a  minimum  value  for  the  size,  k,  of  the  state.  This  minimum 
size  is  often  referred  to  as  the  degree  6  of  the  system;  through 
Eq.  (i-lOc)  we  see  that  6  is  characterized  through  H;  thus  we  can 
write  5[H(p)]  or  (precisely  only  when  E  =  0) 

**  0k  0k 

5  =  min  k  =  &[h(p)]  =  system  degree 

We  will  later  see  how  to  calculate  6  directly  from  JH ( p )  but  for  now 
we  merely  comment  that  5  physically  represents  the  minimum  number  of 
integrators  necessary  for  an  analog  simulation  of  the  system  described 
by  the  canonical  equations  (l-7).  We  do  mention  that  it  is  sometimes 
of  interest  to  have  more  than  the  minimum  number  of  components  of  the 
state  present,  especially  for  the  determination  of  equivalent  realiza¬ 
tions  to  satisfy  some  specified  constraints  (as  for  example  the  desire 
to  incorporate  only  a  certain  type  of  transistor  in  a  design).  Fig¬ 
ure  1-2  has  already  illustrated  an  example  of  a  nonminimaV realization, 
where  we  define  a  minimal  realization  as  one  where  the  A  matrix  is 

_ _  -  .n- 

6x5,  that  is,  has  its  order  equal  to  the  degree  of  H(p).  In  this 
case  h(p)  =  1  / 1* ,  6[n  J  -  0,  and  we  see  that  the  system  of  Fig.  1-2  is 
equivalent  to  a  resistor,  the  situation  being  as  shown  in  Fig.  1-4, 
where  Fig.  I-2a  has  been  redrawn  in  the  (b)  portion. 

E .  The  State  -  Prune  Section  Example 

At  this  point  let  us  set  up  the  canonical  equations  for  the  non¬ 
reciprocal  Brune  section  of  Fig.  1-5  [5,  p.  ],  where  we  make  the 
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(a)  (b) 

Fig.  1-5.  NONRECIPROCAL  BP.UNE  SECTION  (a)  WITH  CAPACITOR  EXTRACTION  (b). 

particular  choice  of  input  and  output  (of  later  use  for  modeling  of 
filters  for  integrated  circuit  realization). 


In  order  to  analyze  the  Brune  section  to  obtain  the  canonical  state 
variable  equations  we  first  separate  the  dynamical  elements  by  removing 
the  capacitors  as  a  load  on  a  purely  resistive  4-port,  as  shown  in 
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Fig,  I-5b.  We  also  take  as  a  convention  for  the  gyrators  the  symbolism 
of  Fig.  1-6. 


g 

0 


Fig.  1-6.  GYRATOR  CONVENTIONS. 

By  summing  currents  at  the  nodes  marked  a,  b,  c  (in  Fig.  1-5)  and 
summing  voltages  around  the  loop  d,  respectively,  we  obtain 

"0  0 
0  0 

_gi  g: 

-l  l 

A  suitable  choice  for  the  state  io  generally  the  set  of  capacitor  voltages 
or  charges  and  inductor  currents  or  flux,  thus  we  let 


°  gi 

o  -g2 
0  0 
1  0 


LVU 


■i  o  i  on 

0  1-10 
0  0  0  1 

oooo 


(l-9b) 


s  = 


C1V3 


C2V4 


for  which  it  follows,  from  Fig.  I-5b,  that 


s  = 


l1! 


( I-9c ) 


( I-9d ) 


We  can  therefore  rewrite  Eq.  (l-9b)  to  specifically  exhibit  the  quantities 
of  interest  by  rearranging  the  columns. 
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0  o 

0  1 


-‘1  0 
-1  0 


-i 


-n 

0 


*2  ° 


0 

0 

0 

1/cl 


Kl/C2~ 

"l  o“ 

r-’i 

= 

“o" 

K2/c2 

LvJ 

-1  0 

L-.J 

0 

0 

0  1 

0 

0 

0  0 

0 

( I-9e) 


If  we  add  the  second  row  to  the  first  and  -g  times  the  last  row  to  the 
third,  we  can  isolate  y  fror.  &  to  get 


1 

H 

o 

1 _ 

“o 

-1  ~ 

o  (8a-B2)/c2 

"*  0 

0 

0  1 

u  + 

0 

0 

y  + 

0  VC2 

S  4 

-1 

0 

Vei  0 

0 

0 

MM 

0 

MM 

0 

1 

-1  0 

1 

0 

° 

0 

0 

(l-9f) 


Using  the  third  row  multiplied  by  -1  and  the  second  row  gives 


0  '"a/'a 

1 

o 

1 _ 

• 

s  = 

S  + 

MM 

.  e2-/ci  0  _ 

_erg2 

(l-9g) 


while  the  fourth  (by  -l)  row  and  then  the  first  give  the  desired  output 
equation 


1 

o 

H 

O 

H 

l  o” 

s  + 

MM 

o  (®1.~®2)/C2 

- 1 

H 

o 

_ 1 

(l-9h) 


These  last  two  equations  are  the  canonical  equations  for  the  Brune 
section. 
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Using  H(p)  =  D  +  C  (pi  -  A)  we  can  find  the  transfer  function. 


T  o“ 

M 

O 

_ 1 

P  82/c2' 

1 

o 

H(p)  = 

+ 

_  0  ^l~g2^C2 

_-b2/ci  P  _ 

ire2  °. 

0 


“I 

=7  p 

P  'e2/C2 

1 

o 

M  1 

1 

+  2  2/ 

P  +  K2/C1C2 

0  (erv 

82/°i  P 

grg2  0 

z 

-  — 

—  — 

+  glic  1C2 


2 

P  + 


glg2 

C1C2 


(gl~g2) 


2 

P 


c  c 
1  2 


(l-9i) 


We  comment  that  one  of  the  alternate  choices  available  for  the  state  is 


and  that  for  this,  or  any  other  choice  for  the  state,  we  obtain  the  same 
transfer  function.  In  fact  we  observe  that  there  is  a  nonsingular  trans¬ 
formation  mapping  one  choice  for  the  state  into  another,  that  is, 


( i-io) 
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As  we  will  later  see,  any  minimal  realization  is  related  to  any  other 
through  a  nonsingular  transformation  on  the  state  as  in  Eq.  (i-lO).  In 
this  case  6  [tjj  =  2,  and  thus  the  realization 


f 

1 

o 

1 

|  <n 

0  1~ 

1 

o 

_ 1 

1 

0 

\ 

C2 

o 

_C1 

t 

g  -  g  0 

81  K2 

» 

-i 

o  Ure2) 

L  C2  J 

l 

0 

1 

( 

is  minimal. 

F.  Discussion 

By  way  of  introduction  (or  review,  depending  upon  previous  back¬ 
ground),  we  have  considered  the  meaning  of  the  state  and  given  the 
primary  equations  related  to  our  further  studies.  For  differential 
systems  the  equations  of  most  interest  are  the  canonical  ones. 

i.=As+B&  (i-lla) 

y  =  C  s  +  D  u  (i-llb) 

with  the  associated  transfer  function  yielding  the  output  y  in  terms 
of  the  input  _u,  when  initially  in  the  zero  state  ^(tQ)  =  0^( t q ) , 
through 


y  = 


h  *  u 


given  by 


£[h]  =  H(p)  =  D  +  C  (pi,  -  A)_1B 


(i-llc) 


(i-lld) 


We  observe  that  in  this  differential  system  case  the  state  is  that  set 
of  parameters  for  which  a  matrix  set  of  first  order  differential  equations 
can  be  set  up  in  terms  of  the  transfer  function  and  its  realization.  The 
matrix  h  is  the  impulse  response  matrix  with  its  Laplace  transform 
£[h]  being  the  transfer  function.  From  the  expression  for  H(p)  in 
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terms  of  the  realization  jl  =  {A.B.C.D)  matrices,  it  is  clear  that  the 
poles  of  H(p)  are  zeros  of  the  determinant  of  pj.^  -  A,  that  is,  the 
natural  frequencies  of  the  system  are  eigenvalues  of  the  matrix  A. 

We  also  observe  that  if  we  have  two  realizations  R  =  {A.B.C.D,} 

**  ty  ** 

and  R  =  {A,B,C  D)  related  through 


A  =  T_1A  T,  B  =  T_1B,  C  =  C  T,  D  =  D 

MM  m  M  ^  M  M  M  AM 


(1-lle) 


with  T  nonsingular,  then  the  two  transfer  functions  are  identical. 
Thus  we  have 

Hi, p)  =  D  +  C  (pi  -  A)  B  =  D  +  C  T  (pT  T  -  T  A  T)  T  B 

'  '  AM  AM  '  *•»!{  •«  '  Aft  K  m  — •  '  MM 


=  D  +  C  (pi,  -  A)-1B  =  H(p) 

W-S  AM  »M*  I'M  AM 


consequently  we  can  investigate  equivalent  systems  by  manipulating  the 
state  variable  e  jations  through  methods  associated  with  the  transforma¬ 
tion  of  Eq.  (i-lle),  which  in  fact  can  be  interpreted  in  terms  of  the 
state  as  a  basis  change  in  the  state  space  through  s^~  2*.  it*  We  are  tlien 
led  to  observ j  that  there  is  a  k -dimensional  space,  the  state  space,  in 
which  we  have  introduced  (Cartesian)  coordinates  against  which  the  com¬ 
ponents  of  ^s  for  the  canonical  equations  are  measured.  The  actual 
state,  for  a  given  input.  u,(t)  and  an  initial  state  s^t^),  traverses 
the  state  space  on  a  trajectory  js(t),  this  trajectory  giving  the  "motion" 
or  behavior  of  the  system,  as  verified  by  Eq .  (l-lla,b). 

Our  primary  interest  will  be  with  linear  networks  considered  as 
systems  through  the  transformation  formulation  so  far  discussed.  One 
could  consider  the  more  general  nonlinear  case  described  by  the  matrix 
differential  equations 


S  =  f  (s,  u,  u) 


( I-12a) 


=  S  (s,  u,  u) 


( I-12b) 


However,  very  little  is  available  in  the  way  of  synthesis  for  such 
equations,  so  we  have  chosen  to  concentrate  on  the  linear  case.  We  also 
choose  to  devote  efforts  primarily  to  the  continuous-time  case  since  it 
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is  of  most  interest  for  network  studies.  But  because  our  treatment  will 


generally  be  of  an  algebraic  nature,  the  results  are  almost  all  valid  for 
discrete-time  systems,  which  in  fact  have  considerable  practical  impor¬ 
tance,  for  example,  through  the  theory  of  automata. 

In  our  treatment  we  have  not  proceeded  in  the  most  rigorous  manner 
possible  since  we  wish  to  bring  out  only  the  basic  and  most  important 
points  for  our  later  use.  Once  the  concepts  we  have  treated  are  grasped 
in  principle,  the  more  detailed  works  are  available  to  those  interested 
[lj.[6].  However  we  have  not  wished  to  sacrifice  completely  the  rigor  of 
the  theory  so  have  proceeded  in  a  rather  precise  manner  for  the  detail 
given.  Although  most  of  our  emphasis  will  he  upon  networks,  we  have 
given  a  somewhat  general  systems  formulation  in  order  not  to  overly  limit 
the  treatment.  As  a  consequence  we  will  most  frequently  work  with  a  net¬ 
work  in  an  input-output  situation,  as  for  example  through  the  admittance 

matrix  where  the  input  u  is  the  set  of  port  voltages  v,  and  the 

01  ** 

output  y  the  port  currents  i  (in  which  case  m  =  n).  Since  such  a 

MM 

(port)  description  tells  very  little  about  the  internal  structure  we  will 
use  the  state  to  discuss  internal  operation  and  construction  of  the 
network.  A  network  is  a  system  with  electrical  inputs  and  outputs. 


It  is  of  interest  to  know 


Li»g  wiC  LauGiuCB.  equations 


so  we  next  turn  to  a  discussion  of  the  setting  up  of  state  variable  equ¬ 


ations. 
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H.  Exercises 


1.  Set  up  the  canonical  equations  for  the  Hazony  section  of  Fig. 
EI-1.  Do  this  for  the  input-output  variables  of  Eq .  (l-9a)  as 
well  as  for  the  admittance  and  impedance  matrices  as  transfer 
funct ions . 

c 

Fig.  EI-1.  HAZONY  SECTION. 

*2.  Given  the  canonical  equations  for  the  admittance  matrix  (as  the 
transfer  function)  and  those  for  the  impedance  matrix,  find  the 
relations  between  the  two  realization  set  matrices.  Repeat  for 
the  scattering  matrix  and  the  admittance  matrix  given. 

3.  A  given  network  has  the  canonical  equations 


,  =  [j  -g  . *  CO  “ 


(El-lb) 


*4. 

5. 


6. 


Find  the  transfer  function. 


...... .O.C.  *  .w...  1 

b.  Find  the  zero  input  response  for  s(t  )  =  I  I.  Plot  the 

,  .  °  L-2J 

trajectory  s^t;  in  state  space. 

Discuss  a  formulation  for  "transfer  functions"  in  terms  of  the 


realization  matrices  for  time-variable  networks. 


As  we  have  mentioned,  the  state  applies  to  much  more  than  sci¬ 
entific  or  physical  systems.  Investigate  the  concept  in  terms 
of,  for  example,  language  formation  or  motion  picture  production. 
Consider  any  network  of  interest  and  set  up  the  appropriate  state 
space  equations.  From  these,  investigate  the  minimality  of  the 
realization  as  well  as  other  sets  of  canonical  equations  yielding 
the  desired  transfer  function. 
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Mais  d'autres  s',att  ardent  un  peu ,  nous 
regardent  en  souriunt  att ent ivement , 
seinblant  sur  le  point  d'avouer  qu'ils  ont 
tout  compris. 

M.  Maeterlinck 

"Les  Avert  is"  du  "Tresor  des  Humbles" 


CHAPTER  II 

FORMULATION  OF  CANONICAL  EQUATIONS 


A.  Summary 

By  the  tise  of  appropriate  replacements  and  capacitor  extractions  a 
simple  method  of  equation  formulation  suitable  for  digital  computer  use 
is  presented;  the  method  is  described  in  terms  of  the  admittance  descrip¬ 
tion  but  can  be  used  in  other  situations.  This  method  is  followed  by 
the  outline  of  a  topological  one  which  exhibits  a  more  general  set  of 
equations. 


B.  Capacitor  Extractions 

Let  us  consider  as  given  a  finite  circuit,  that  is,  a  connection  of 
a  finite  number  of  resistors,  capacitors,  indue. ors,  transformers,  gyra- 
tors,  and  devices,  such  as  transistors,  which  cun  be  modeled  by  the 
above  elements.  (We  assume  linear  but  perhaps  time  variable  and  active 
elements  at  this  point;  that  is,  negative  as  well  as  positive  element 
values  which  may  vary  with  time  are  allowed.)  To  illustrate  the  method, 
we  search  for  the  canonical  state  variable  equations  for  the  admittance 
matrix  as  transfer  function  [lj.  To  concentrate  on  fundamental  concepts, 
we  replace  all  inductors  by  the  capacitor-loaded  gyrator  equivalent 
shown  in  Fig.  II-l. 

After  making  such  a  replacement  we  extract  all  capacitors  into  a 
separate  network  which  loads  a  multiport  described  completely  by  alge¬ 
braic  constraints.  If  the  admittance  matrix  is  n  X  n  and  if  there  are 
c  capacitors  extracted,  the  situation  is  shown  in  Fig.  I I -2 ,  where 
the  "resistive"  (n  +  c)-port  is  loaded  by  a  capacitive  c-port. 
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(b) 

Fig.  II-l.  INDUCTOR  EQUIVALENT. 
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Here  we  have  taken  the  matrix  ^  as  the  c  X  c  diagonal  matrix  of 
capacitance  values;  any  capacitive  coupling  we  assume  to  have  been  taken 
into  account  through  transformers  absorbed  into  the  resistive  (n  +  c)- 
port .  We  also  assume  £(t)  to  be  nonsingular.  (Any  singularity  can 
actually  be  accounted  for  again  by  a  change  in  the  (n  +  c)-port,  but 
we  omit  discussion  of  this  rather  tricky  point  in  order  to  clearly  pro¬ 
ceed.)  Substituting  the  (b)  and  (c)  portions  of  Eq .  (lI-2)  into  the  (a) 
one  yields 
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The  second  set  of  (c)  equations  gives  the  derivative  portion  of  the 
canonical  equations,  while  the  first  set  of  n  equations  gives  the 
output  portion.  Thus, 


£  =  -.222  £  3.  -  Jfel  2, 


( ll-4a) 


1  ■  Z12  £  4  +  2.H  & 


( ll-4b) 


We  have  obtained  the  realization 


R  ={“^22  ~_1,  12V  ~12  2*1’  «2ll| 


( II-4c) 


in  a  simple  manner.  It  is  worth  mentioning  that  if  time-variable  elements 
are  present  the  realization  matrices  are  functions  of  time,  in  which  case 
we  have  succeeded  in  setting  up  the  canonical  state-variable  equations  for 
time-variable  circuits.  In  the  time-invariant  situation  we  observe  that 
the  method  proceeds  only  when  there  is  no  pole  in  the  (n-port)  admittance 
matrix  at  infinity;  we  will  later  (Sec.  C)  obtain  a  graph  theory  con¬ 
dition  for  no  pole  at  infinity  such  that  a  test  can  be  directly  made  on 
the  circuit  graph.  In  any  case,  time-variable  or  not,  the  method  proceeds 
if  and  only  if  the  coupling  admittance  matrix  Y  exists;  the  existence 
of  Y  is  equivalent  to  the  existence  of  the  inverse  of  the  $  matrix  in 
the  general  description,  v^  =  ®  i^,  for  the  (n  +  c)-port  coupling  network. 

As  an  example,  let  us  consider  the  2-port  of  Fig.  II-3,  which  is  a 
subportion  of  the  nonreciprocal  Brune  section,  useful  for  its  own  sake 
(since  it  is  equivalent  to  a  series  inductor  in  cascade  with  a  trans¬ 
former)  . 


Fig.  I 1-3.  SHUNT -CAPACITOR  LOADED-GYRATOR  CASCADE. 
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By  extracting  the  capacitor  as  shown  in  the  (b)  portion  of  the  figure, 
we  can  obtain  the  appropriate  equations.  First  we  write  the  general 
description  for  the  3-port  coupling  structure  (by  respectively  summing 
currents  at  node  a  and  then  writing  i^  and  i^  in  terms  of 
through  the  gyrator  relationships). 
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( II-5a) 


The  coefficient  matrix  of  the  currents  is  nonsingular,  being  a  permuta¬ 
tion  matrix,  and  thus  on  premultiplying  Eq.  (ll-5a)  by  its  inverse  we 
find 
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(ll-5b) 


where  we  have  made  the  partition  appropriate  to  Eq.  (il-l).  Note  that 
is  skew-symmetric,  =  (where  means  transpose),  as  ex¬ 

pected,  since  it  is  constructed  solely  from  gyrators. 

Equation  (lI-3)  is  directly 
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( II-5c ) 


where  we  have 
Thus  we  have, 


partitioned  the  last  c  =  1  equations  to  be  split  off. 
by  such  a  split,  the  canonical  equations  directly  as 


0-s  +  [gx  -g2J 


(ll-5d) 
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We  comment  that,  since  the  nonreciprocal  Brune  section  itself  has  a  pole 

at  infinity,  no  Y  exists  for  it.  However,  on  removal  of  the  pole  at 
**c 

infinity,  Eqs.  (lI-5)  result;  hence;  the  canonical  equations  for  the 
admittance  description  of  the  Brune  section  are  merely  obtained  from 
Eqs.  (ll-5d,e)  by  adding 


to  the  right  of  Eq.  (ll-5e).  Note  also  that  the  canonical  equations 
previously  found  for  the  Brune  section  were  for  a  different  set  of  input- 
output  variables  (that  is,  a  different  system).  Still  the  same  method 
was  applied  at  that  point. 

We  also  comment  that,  upon  adding  suitable  ports  and  ignoring  vari¬ 
ables  of  no  interest,  we  can  use  the  same  method  to  find  almost  any 
input-output  canonical  set  of  state-variable  equations,  perhaps  also  after 
simple  transformations  on  the  variables.  This  result  is  directly  seen  by 
setting  up  equations  in  hybrid  form. 

Since  the  steps  carried  out  are  easily  programmed,  the  procedure  is 

a  very  convenient  one  for  use  in  setting  up  canonical  equations  on  a 

digital  computer.  For  such  purposes  one  needs  a  method  for  obtaining 

Ulr  coupling  admittance  Y  on  the  computer.  Perhaps  the  most  co  e- 

Mow!  method  is  to  reduce  the  indefinite  admittance  matrix  [2,  p.  sj 

for  the  resistive  coupling  network  to  obtain  Y  ;  several  programs  are 

•sc 

available  for  finding  the  indefinite  admittance  matrix,  but  a  program  is 
also  very  easily  written  from  scratch.  An  alternate  and  almost  equally 
useful  method  is  to  use  the  topological  methods  which  we  now  discuss. 
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C.  Topological  Formula l ion 

Let  us  again  consider  a  finite  circuit  lor  which  the  equivalence  of 
Fig.  II-l  is  used  to  replace  inductors;  again  this  replacement  is  not 
necessary  but  is  convenient  for  simplification  of  already  complicated 
expressions.  Also  we  will  assume  that  the  admittance  description  is 
desired  for  which  voltage  sources  have  been  placed  at  the  ports. 

By  replacing  each  circuit  element  branch  by  a  line  segment,  with  an 


arbitrarily  assigned  orientation,  as  shown  in  Fig.  II-4,  we  obtain  an 
oriented  graph  to  represent  the  circuit,  the  branches  of  which  we  can 


Fig.  I I -4.  EXAMPLE  GRAPH  REPLACEMENTS. 


number  in  some  useful  manner.  A  graph  associated  with  a  network  or 
circuit  structure  will  be  called  a  network  graph. 

In  order  to  proceed  we  introduce  the  following  somewhat  standard 
nomenclature  associated  with  a  network  graph: 

node=vertex  a  dot  on  the  graph  (=  a  terminal  of  a  circuit 

element  branch) 

branch  a  line  connecting  two  nodes  (=  a  circuit  element 

branch) 

path  a  sequence  of  branches  and  associated  nodes 

connected  graph  a  network  graph  in  which  every  node  is  connected 
to  every  other  node  by  a  path 

separate  part  a  maximally  connected  subgraph  (that  is,  a  sub¬ 
graph  for  which  all  branches  are  connected  to  all 
other  branches  in  the  subgraph  and  to  no  others) 

tree  a  maximally  connected  subgraph  of  a  separate  part 

which  contains  no  closed  path 

forest  a  collection  of  trees  of  a  graph,  one  for  eack 

separate  part 

coiree  the  set  of  branches  (in  a  separate  part)  which 

remain  when  a  (fixed)  tree  is  deleted 

link  a  branch  of  a  (fixed)  cotree 
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Although  these  definitions  are  not  completely  rigorous  (for  example 
"connected"  and  "closed  path"  are  not  made  precise),  they  should  be  in¬ 
tuitively  clear,  perhaps  after  an  example,  and  are  sufficient  for  our 
purposes.  To  further  proceed  we  introduce  the  following  symbols: 

b  =  total  number  of  branches 
l  =  total  number  of  links  (cotree  branches) 
s  =  number  of  separate  parts 
t  =  total  number  of  tree  branches 
v  =  total  number  of  nodes 

Here  l  and  t  are  formed  by  summing  over  all  i.rees  in  a  forest.  For 
each  separate  part  the  number  of  tree  branches  is  o  le  less  than  the 
number  of  nodes  while  it  is  also  clear  that  b  =  l  +  t.  Thus  we  can 
directly  predict  the  number  of  tree  branches  and  links,  without  express¬ 
ly  exhibiting  a  tree,  through 

t  =  v  -  s,  i=b-v  +  e  ( II-6) 

As  an  example,  let  us  consider  the  2-port  of  Fig.  II-5  which  has 
been  closed,  as  mentioned  above,  on  voltage  sources  (as  will  be  appro¬ 
priate  to  setting  up  the  canonical  equations;  note  that  this  network  is 
identical  in  port  behavior  to  the  nonreciprocal  Brune  section  of  Fig. 
1-5).  A  possible  network  graph  is  shown  in  the  (b)  portion,  with  other 
graphs  resulting  by  different  choices  of  branch  orientation  and  number¬ 
ing.  Note  that  by  simple  count  b=8,  s  =  2,  v  =  5,  and  thus,  by 
Eq.  (lI-6),  t  =  3  and  i  =  5;  these  numbers  are  checked  from  the  graph 
where  a  possible  choice  for  a  tree  is  shown  in  boldface  (note  that  there 
are  other  choices  for  a  tree,  but  that  in  a  given  analysis  only  one  at  a 
time  is  need). 


(a)  (b) 

Fig.  II-5  EXAMPLE  GRAPH  FROM  CIRCUIT  STRUCTURE. 
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Since  we  will  wish  to  sum  currents  .it  the  imdis,  vie  have  also  labeled 


them.  We  observe  that  lor  a  node  analysis  we  wish  to  choose  tree  branch 
voltages  as  independent  variables  while  lor  a  loop  analysis  we  wish  to 
choose  link  currents.  In  setting  up  the  state-variable  equations  we 
actually  will  work  with  both  types  of  variables. 

Next  we  introduce  the  following  (column)  vector  variables: 

i.  =  vector  of  branch  currents  (b  X  l) 

-b 

v  =  vector  of  branch  voltages  (b  X  l) 

«*b 

i,  =  vector  of  link  currents  (i  X  l) 
v  =  vector  of  tree  branch  voltages  (t  X  l) 

;v  =  vector  of  port  (source)  voltages  (n  X  l) 

=  vector  of  port  currents  (n  x  l) 

Along  with  these  variables  we  assume  the  polarity  of  a  given  branch’s 
variables  in  conjunction  with  the  given  branch  orientation  as  shown  in 
Fig.  I I -6. 

i 

Fig.  11-6.  POLARITY  OF  VARIABLES. 


Now  we  introduce  the  cut  set  and  tie  set  matrices  from  which  the 
analysis  can  truly  begin.  For  a  given  circuit  we  pick  a  fixed  forest. 
The  cut  set  matrix  c  is  defined  by  considering  the  tree  branches  in 
numerical  order;  for  each  tree  branch  a  circle  (or  similar  curve)  is 
drawn  such  that  of  all  the  tree  branches  only  the  prescribed  one  is  cut 
by  the  circle.  The  (oriented)  set  of  branches  cut  by  any  one  circle  is 
called  a  cut  set .  For  any  one  cut  set  all  the  currents  entering  the 
circle  on  the  cut  set  branches  must  sum  to  zero  by  Kirchhoff's  current 
law;  considering  all  t  cut  sets  we  obtain 


0  =  (7  i 

Mi  AMD 


(II-7) 


where  ('  is  the  t 
*\ 

As  an  example.  Fig. 


X  b  cut  set  matrix  (consisting  of  0  or  ± 

I I —7  shows  the  cuts  for  the  particular  graph. 


I's). 

The 
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resulting  cut  set  matrix  is  given  as  the  coefficient  matrix  in  the 
equation 


cut  I  -* 
cut  II  -* 
cut  III  -* 


1  0 
0  1 
0  -1 


( I 1—8) 


III 

cut  set  I  =  branch  1 

cut  set  II  =  branch  2  (out)  and  branch 

3  (out) 

cut  set  III  =  branch  2  (in)  and  branch 

4  (in) 


II 


Fig.  II-7.  EXAMPLE  COTS  FOR  C. 


The  tie  set  matrix  ‘l  is  defined  in  a  somewhat  dual  manner.  Again 

DM 

a  forest  is  chosen.  On  removing  all  links  a  particular  link  is  rein¬ 
serted;  the  (oriented)  branches  forming  a  closed  path  under  this  rein¬ 
sertion  are  the  associated  tie  set.  Ordering  all  tie  sets  according  to 
the  numerical  order  of  the  links  defines,  through  Kirchhoff's  voltage  law 
(applied  to  each  loop  of  lie  set  branches), 


0  =  3  v 

MM  ** 


( II— 9) 


where  ^  is  the  £  X  b  tie  set  matrix  (again  consisting  of  0  or 
±  l's).  For  example,  Fig.  II-8  has 


( 11-10) 
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Fig.  1 1 -8. 


Lie  set  I  =  branch  2  (+)  and  branch  1  (+) 

tie  set  II  =  branch  3  (-t),  branch  1  (  +  ) 
and  branch  4  (-) 


EXAMPLE  TIES  FOR  J . 


We  also  claim  that  it  is  possible  to  write  (recall  that 
transpose) 


V. 

=  C  v 

( Il-lla ) 

Mtb 

•m  nnt 

«ib 

=  i  > 

**  m  JO 

( II-llb) 

For  the  plausibility  of  Eq .  (il-lla),  say,  let  us  argue  as  follows.  By 

Kirchhoff's  voltage  law  it  should  be  clear  that  the  tree  branch  voltages 

determine  ail  link  voltages;  hence  there  is  a  linear  transformation  to 

eive  Xh  =  where  A  is  some  b  X  t  matrix  in  fact  consisting  of 

zeros  and  (+  or  -)  ones.  If  we  consider  the  graph  as  a  closed  system 

then  the  total  input  power  is  zero,  P,  =  i  =  0.  Thus  v.  i.  = 

^  in  **0  b  **b 

v  A  i  =  0.  Since  the  tree  branch  voltages  can  be  arbitrarily  assigned 
y*  £  "*  ***  b 

(when  the  graph  is  considered  as  an  abstract  object),  we  must  require 
^  =  JJ.  In  other  words  if  we  choose  Eq .  (il-lla),  then  Eq.  (lI-7) 

follows  as  a  possibility.  [Of  course,  a  proof  requires  that  we  argue  in 
reverse,  but  this  can  be  done  by  beginning  with  Eqs.  (il-ll)  at  first.] 
For  convenience  of  notation  we  next  choose  a  numbering  of  branches 
such  that  all  the  branches  occur  first;  thus 

(il-llc) 
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in  which  case  Eqs.  (ll-lla,b)  show  that  the  cut  set  and  tie  set  matrices 
can  be  partitioned  as 


c  = 


[t  !  l-j 


(il-lld) 


where  ,£  and  T  are,  respectively  t  x  I  and  l  X  t  matrices;  ^  is, 
of  course,  the  t  X  t  identity.  We  observe  that 


£  = 


(il-lle) 


since  again 


ihi 


b*b 


0  =  iJ1*  £] 

t  >-t  — 


=  It[T  +  C]JLX  =  0 


and  v^  and  can  be  arbitrarily  assigned. 

Our  next  step  is  to  place  all  voltage  sources  in  tree  branches. 

(We  remark  that  we  are  only  considering  the  presence  of  voltage  sources; 
if  current  sources  are  present,  only  simple  modifications  are  necessary, 
or  one  can  use  the  equivalence  of  Fig.  11-9.)  Next  we  place  as  many 


1 


Fig.  I I -9.  CURRENT  SOURCE  EQUIVALENT. 


as  possible  of  the  capacitors  in  tree  branches — any  left  over  are  some¬ 
how  "excess";  but  the  need  for  considering  these  excess  capacitors  is  in 
fact  the  reason  for  our  treating  this  topological  method.  It  follows 
that  if  there  is  a  capacitor  link  then  the  path  formed  by  the  associated 
tie  set  branches  consists  entirely  of  voltage  sources  and  capacitors-- 
such  gives  rise  to  a  pole  at  infinity,  for  example,  in  the  admittance 
matrix.  Let  us  now  further  fix  our  numbering  of  branches  such  that 

v  and  vfl  take  the  form 
**>t  'mx> 
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V 


— t 


(11-12) 


where  the  subscripts  c  and  r  refer  to  capacitor  and  resistor  portions 
of  the  graph. 

At  this  point  we  can  begin  the  real  procedure  for  setting  up  the 
state-variable  equations  [3j.  If  we  partition  the  matrix  ^F,  of  Eq. 
(il-lld) ,  using  C  =  -J,  we  find  for  Eq.  (lI-9),  i2,  =  Jxb»  and  for 
Eq.  (II-7),  0=0^ 
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(ll-13a) 
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(ll-13b) 


In  these  equations  T  =  0  since  if  there  is  a  capacitor  in  a  link  then 
there  is  no  resistor  in  the  tree  branches  of  the  associated  tie  set.  We 
also  write  i  for  the  source  current  and  note  i  =  -i,  where  i  is 

— S  **«  M  — 

the  port  current. 
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The  circuit  element  constraints  are  next  assumed  to  be  of  the 


form 


*-ct 

d£At/dt 

^ci’ 

d£lSe«/dt 

•irt 

St-rt 

1 

•Hf 

1 

In  actual  fact  this  form  places  some  restrictions  on  the  types  of  cir¬ 
cuits  allowed  since  no  coupling  between  tree  branch  and  link  resistive 
(gyrator)  elements  Is  allowed;  for  example,  the  circuit  of  Fig.  11-10 
is  ruled  out.  Of  course  a  more  general  treatment  is  possible  by  using 


*"*rt 

St 

St  i 

1 - 

•*-> 
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>1 

i _ 

i  9 
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Set 

Se 
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>i 

_ i 

( II-14b) 


but  as  we  will  see,  the  result  is  already  comol  i<~at  od  enough  irt  notation. 


note:  1  and  6  are  required  tree  branches 

Fig.  11-10.  EXAMPLE  OF  RESISTIVE  COUPLING  BEWEEN 
TREE  BRANCHES. 

Our  next  job  is  to  make  appropriate  substitutions,  etc.  Through 
the  various  equations  Indexed  as  shown  we  can  write  the  right  side  of 
Eq.  (ll-14a)  as 
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and 
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( 1 I-16a) 


Equation  (ll-log)  is  the  same  as 


s  =  A  s  +  B  u  +  F  u 

im  Mt  ^  mm  mm 


where  we  have 


s  =  C  v  u  =  v 

**»  ■*«t'»ct ,  <■»  *•» 


(  1 1  - 1  fib) 


£  =  [let  + 


-n"M  ^  d££Mi2St 


T  N 


12  dt 


+  <X22£i[”«iri  +  ^23^rt+  £t  l23^«W  £t  ^2^.  ^22^t 

r  ~  -n1!  ~  d£i 

~  =  [let  +  J  |"*Il2  dt  *11 

+  Z,22>U  [~ir£  +  **23^irt  +  £t  I23^l23^  -t  ^23^]  ~21  j 
£~  ~  [ict  +  ~12>Ull2*t  ]  Il2£illl 


(11-16. 


For  the  output  equations  we  can  use  the  source  portion  of  Eq.  (ll-13b) 


to  get 


-i  =  T, ,  i  „  +  T  i  „ 


(11-17) 


But,  from  Eq .  (ll-14a),  =  d£^|/dt  ir  £  =  ZgZrt''  the  ici) 

part  can  be  evaluated  from  Eq .  (ll-15a),  and  the  Jl  ^  part  from  Eq .  (il- 
15f ) .  Thus  we  find 


v=Cs+Du+Eu 

m*  9m  9m  9m  **  ** 


( II-l 8a ) 


( II-18b) 


~  d£^Ii2£t  ~  -i 

& =  «ln  dt  +  In££«i2£t  £  "  l2i^ 


X  ["iri  +  «23^1rt  +  £t  ^23^23^  ~t 
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Nevertheless  if  C*  =  0  then  F  =  0  and  E  =  0,  and  thus,  when  there 
are  no  cepacitor-source  tie  sets,  we  obtain  the  canonical  equations.  It 
should  be  observed  that  the  results  are  valid  for  time-variable  elements 
and  that  the  only  real  restriction  on  the  result  is  the  requirement  that 
there  be  no  resistive  coupling  between  tree  branches  and  links,  that  is, 
zero  and  in  Eq .  (ll-14b). 

Even  in  the  time-invariant  case  where  there  are  no  capacitor-source 
tie  sets,  where  considerable  simplification  occurs,  the  equations  still 
remain  rather  messy.  Thus  we  observe  that,  although  the  formulation  is 
important  for  illustrating  the  general  nature  of  network  state-space¬ 
like  equations,  the  approach  is  not  the  most  useful  to  be  taken  for 
normal  analysis  or  synthesis. 

As  an  example,  let  us  reconsider  the  nonreciprocal  Brune  section  of 
Fig.  1-5.  This  is  redrawn  in  Fig.  11-11,  where  the  appropriate  tree  is 
shown  with  the  numbering  requested  by  the  theory. 
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The  tie  set  and  cut  set  matrices  are  found 
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(  1 1  - 1 9c ) 
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( I I-19e) 


The  result  is  checked  by  calculating  the  transfer  function  matrix  in  the 
time-invariant  case. 


D.  Transformation  to  Canonical  Form 

Because  of  the  presence  of  JF  in  the  resultant  topological  equations, 
it  is  of  interest  to  find  a  transformation  to  eliminate  the  derivative  of 
the  input  in  the  differential  equation  for  the  state.  Fo~  this  let  us 
assume  that  we  have  on  hand  a  set  of  equations 

X  =  Ax  +  Bu  +  Fu  ( II-20a) 

y  =  Cx  +  Du  +  Eu  (ll-20b) 

NO  w*  M  Vt*  M 

The  transformation 

(ll-20c) 

leads  to  the  canonical  set 

s  =  A  s  +  (B  +  F  -  F)u  (ll-20d) 

y  =  C  s  +  (D  +  f)u  +  E  u  (ll-20e) 

*«W\  **  «***»*"* 

v.'e  observe  that  such  a  transformation,  for  which  the  input  becomes  part 
of  the  state,  loaves  the  A.,  (2,  and  matrices  unchangtd. 

E.  Combination  of  Methods 

If  one  applies  the  topological  method  to  e  purely  resistive  structure, 
the  results  are  considerably  simplified.  In  the  cases  where  there  is  no 
coupling  betveen  tree  branches  and  links,  one  merely  has  that  the  admit¬ 
tance  is  given  by  JO  of  Eq .  (ll-18c).  We  poini  out  that  the  operations 


x  =  s  +  F  u 

•»  Mm. 
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to  obtain  J3  are  in  this  resistive  case  relatively  easy  to  set  up  on  a 
computer.  Hence  if  capacitor  extractions  are  first  made  and  then  a 
topological  analysis  carried  out  on  the  resulting  resistive  structure,  a 
very  convenient  method  of  setting  up  state-variable  equations  via  the 
computer  results.  The  method  is  also  quite  easily  extended  to  cover 
those  cases  where  there  is  resistive  coupling  between  links  and  tree 
branches . 

By  first  setting  up  the  graph  of  the  circuit,  the  topological  ap¬ 
proach  can  be  used  to  check  the  circuit  for  capacitor-source  tie  sets  to 
establish  the  existence  of  the  E  matrix.  If  there  Bre  such  tie  sets, 
the  topological  formulation  to  calculate  ]£  can  actually  be  carried  out-- 
the  last  of  Eq .  (ll-18c) — since  this  calculation  in  itself  is  not  too 
dil f icult . 

F.  Discussion 

Because  we  feel  it  important  to  understand  somewhat  more  fully  how 
state-variable  equations  can  arise,  as  well  as  more  of  their  meaning,  we 
have  presented  two  convenient  methods  of  setting  up  the  canonical  equa¬ 
tions.  Although  both  methods  cover  most  situations  of  inte'  st  and  have 
been  presented  for  tie  time-varying  case,  neither  one  is  in  itself  com¬ 
pletely  general.  The  capacitor  extraction  method  is  lacking  in  that 
there  can  be  no  capacitor-source  tie  sets  in  the  circuit,  while  the 
topological  method  needs  to  be  extended  to  cover  the  case  where  non- 
dynamical  (that  is,  resistive)  portions  have  coupling  between  the  tree 
branches  and  the  links.  The  capacitor  extraction  method  has  the  advan¬ 
tage  of  simplicity  while  the  topological  method  has  the  advantage  of 
proceeding  directly  from  the  circuit  structure.  When  the  two  methods 
are  combined  by  applying  the  topological  techniques  to  the  nondynamic 
portions  resulting  from  the  capacitor  extractions,  an  excellent  method 
appropriate  for  computer  analysis  of  networks  results. 

To  this  point  we  have  not  commented  upon  the  existence  of  various 
inverses  needed  in  the  topological  approach.  To  investigate  these  would 
cause  an  inappropriate  diversion  so  we  merely  mention  that  in  the  case  of 
passive  time-invariant  circuit  elements  all  inverses  are  known  to  exist 
[3,  p.  511 J . 
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In  many  applications,  especially  lor  integrated  circuits,  one  meets 
voltage-controlled  voltage  sources.  By  changing  somewhat  the  theory, 
these  can  be  handled  directly,  but  for  our  purposes  it  is  worth  observing 
that  the  topological  theory  presented  applies  if  one  is  willing  to  use 
the  equivalence  of  Fig.  11-12,  for  which  each  of  the  cascade  portions 
posseses  a  conductance  matrix. 


1  1 


Fig.  11-12.  CONTROLLED  SOURCE  EQUIVALENT. 


Since  the  topological  method  is  in  itself  a  bit  complicated  in  end 
results,  it  is  of  interest  to  note  that  the  results  are  almost  identical 
to  those  obtained  by  Bryant  [4]  by  very  similar  means. 

Our  next  step  will  be  to  reverse  the  procedure  and  set  up  a  physical 
realization  from  a  state-variable  realization. 
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H.  Exercises 


1.  Set  up  canonical  state-variables  equations  for  the  filter 
circuit  of  Fig.  EII-1. 


‘1 

■MV" 


u  =  vA 


X 

— 

- - 

X 

)( 

C1 

C2 

=  y 


Fig.  EII-1.  DEGREE  TWO  FILTER. 


2.  Set  up  the  canonical  state-variable  equations  for  the  classical 
degree  two  feedback  section  of  Fig.  EII-2. 


Fig.  EII-2.  DEGREE  TWO  FEEDBACK  SECTION. 


*3.  Develop  a  method  for  the  analysis  by  topological  means  of  the 
general  resistive  structure  coming  from  the  capacitor  extrac¬ 
tion  method  such  has  been  proposed  by  E.  Bailey  [5j. 

4.  Set  up  the  canonical  state-variable  equations  for  the  inte¬ 
grated  circuit  integrator  of  Fig.  III-4b  and  investigate  var¬ 
ious  transformations  on  the  resultant  equations. 

*5.  Investigate  the  existence  Of  the  inverses  needed  to  form  A, 

•*» 

JB,  £  of  Eq .  (ll-16c).  From  such  an  investigation,  exhibit  an 
example  of  a  circuit  with  no  canonical  set  of  state-variable 
equations.  Further,  investigate  the  set  of  equations  needed 
to  be  discussed  such  that  all  circuits,  active  or  passive  but 
with  differential  equation  descriptions,  are  covered. 

6.  Set  up  the  canonical  equations  by  the  topological  method  with¬ 
out  using  gyrator  replacements  when  only  inductors  and 
capacitors  (as  well  as  sources)  are  present. 
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A  la  hate,  sagement  et  rainutieusement , 
ils  se  preparent  a  vivre. 

Et  puis,  vers  la  vingtieme  annee, 
s'eloignent  a  la  hilte,  en  etouffant  leurs 
pas,  comme  s' ils  venaient  de  d^couvrir 
qu'ils  s'etaient  trompes  de  deraeure  et 
qu'ils  allaient  passer  leur  vie  parmi  des 
hommes  qu'ils  ne  conne.issaient  pas. 

M.  Maeterlinck 

"Les  Avert  is"  du  'Tresor  des  Humbles" 


CHAPTER  III 

INTEGRATED  AND  ANALOG  CIRCUIT  CONFIGURATIONS 


A .  Summa ry 

The  canonical  equations  are  convenient  for  system  simulation,  espe¬ 
cially  through  the  use  of  integrated  circuits.  Here  we  discuss  the  con¬ 
cepts  of  interest  in  terms  of  appropriate  integrated  circuit  configurations 
In  the  development  special  operational  amplifier  circuits  are  considered 
to  illustrate  some  of  the  points  associated  with  integrated  circuit  struc¬ 
tures  . 

B .  Canonical  Equation  Simulation  -  Block  Diagram 

Let  us  consider  the  canonical  equations  of  the  form 

i  =  &  &  +  3  a  (I-7a ) 

( I-7b) 

where  the  dot  has  been  used  to  denote  time  differentiation.  If  we  integrat 
these  canonical  equations  while  denoting  the  (zero  state)  integral  operator 
as  l/p,  that  is, 


y  =  C  s  +  D  u 

A  m  FI  •>  »* 


1 

P 


J  d~ 


( Ill-la ) 
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then  we  arrive  at  the  useful  equations  for  analog  simulation 


&  =(jUk)  + 1  &J  dii-ib) 

r  Ss  +  n  (in-ic) 

where,  as  before,  is  the  k  >  k  identity,  the  state  £  being  a 

vo. 't or.  For  any  input  the  system  can  be  simulated  from  a  given 
real  iTai  ion  K  =  such  that  the  output  is  determined 

by  the  block  diagram  of  Fig.  1II-1.  Note  that  since  the  various  sub¬ 
systems  are  multidimensional,  the  separate  blocks  have,  in  general, 
multiple  inputs  and  outputs. 


Fig.  III-l.  BLOCK  DIAGRAM  FOR  CANONICAL  EQUATIONS. 


Several  things  can  be  noted  concerning  Fig.  I I 1—1 : 

1.  Positive  feedback  is  used  and  hence  for  (asymptotic)  stability 
we  require  ^A  to  have  all  of  its  eigenvalues  negative. 

2.  Except  for  the  integrators,  all  blocks  consist  simply  of  gain 
elements.  Such  multidimensional  gain  blocks  can  be  constructed 
by  interconnecting  one-dimensional  gain  blocks,  as  shown  for 
example  in  Fig.  III-2  for  the  2-input,  3-output  case.  We  shall 
later  see  a  method  of  summing,  with  gain,  many  inputs  using  a 
single  amplifier,  but  at  this  point  remark  that  the  gain  blocks 
as  well  as  summers  need  consist  only  of  operational  amplifiers 
and  resistors. 

3.  All  integrators  are  uncoupled  and  of  unity  gain.  In  practice, 
and  especially  with  integrated  circuits,  nonunity  gain  integra¬ 
tors  must  be  used,  necessitating  a  scale  change.  Since  It  is 
most  convenient  to  construct  all  components  identical  with  in¬ 
tegrated  circuits,  it  is  practically  more  useful  to  simulate  t V. ». 
system  through  the  equations 
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H 


( 1 1 1-2a ) 


J  =  £  a  +  fi  a. 


(lll-2b) 


where  A  is  an  appropriate  gain  constant.  A  simulation  of 
these  latter  equations  is  just  as  for  the  previous  ones  except 
that  the  integrator  and  .A,  £  blocks  are  scaled  in  gain. 

4.  For  most  practical  simulations  it  is  customary  to  use  voltages 
as  varieties,  in  which  case  all  gains  are  for  voltage  transfer 
elements . 

5.  Time-variable  realizations  R  are  allowed,  in  which  case  it  is 
of  interest  to  observe  that  our  use  of  p  is  as  a  differential 
operator  and  not  as  the  Laplace  transform  variable. 


Fig.  III-2 .  THE  2-INPUT,  3-OUTPUT  GAIN  BLOCK. 


C .  Integrators  and  Summers 

In  order  to  simulate  the  canonical  equations  we  see  that  it  is  of 
Interest  to  have  gain  blocks,  integrators,  and  summers.  A  glance  at 
Fig.  III-2,  as  well  as  the  manner  in  which  summation  occurs  in  Fig.  Ill-)  , 
shows  that  the  gain  portions  can  be  incorporated  in  the  summers.  Con¬ 
sequently,  we  concentrate  upon  one-dimensional  integrators  and  multiple- 
input,  single-output  summers  with  emphasis  upon  structures  suitable  for 
integrated  circuits. 

The  basic  building  block  is  the  operational  amplifier.  For  inte¬ 
grated  circuits  one  likes  to  use  symmetrical  structures  with  equal  resis¬ 
tors,  with  quantities  ol  interest  determined  by  ratios  of  resistors  in 
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place  of  absolute  values,  where  possible.  Likewise  one  generally  avoids 
pnp  transistors  where  possible  because  of  processing  problems  associated 
with  making  both  npn  and  pnp  transistors  simultaneously.  One  is 
therefore  led  to  consider  the  basic  operational  amplifier  structure  of 
Fig.  III-3,  on  which  many  refinements  are  made  to  obtain  various  types 
of  improvements,  as  higher  gain  by  cascading  of  input  amplifiers.  For 
reasonable  values  of  R,  larger  than  the  T emitter-base  resistance 
(say  R  es  3  kii),  the  gain  of  the  device  is  roughly  [l,  p.  J 

K  V.  U  10  V.  at  room  temperature]  ( II 1-3) 

4kT  b  b 


where  q  =  electron  charge,  k  =  Boltzmann's  constant,  T  =  absolute  tempera¬ 
ture.  We  observe  that  a  differential  amplifier  is  obtained,  this  being 
convenient  for  summers  which  both  add  and  subtract.  On  the  circuit 


diagi-am  some  of  the  dc  voltages  have  been  indicated  for  convenience  with 
the  input  voltages  v+  and  v  assumed  held  at  zero  volts  dc  by  exter¬ 
nal  circuitry  under  the  application  of  no  Bignal.  The  zener  diode  is 


Fig.  II I -3.  BASIC  DIFFERENTIAL  OPERATIONAL  AMPLIFIER. 


One  can,  of  course,  use  the  standard  capacitor  feedback  structure 
for  integration,  as  shown  in  Fig.  III-4a,  but  if  a  completely  integrated 
device  is  desired,  which  includes  integrated  capacitors,  then  it  is  most 
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R 


(a)  (b) 


Fig.  I I 1-4.  POSSIBLE  INTEGRATORS. 

convenient  to  use  the  integrator  of  Fig.  III-4b,  which  in  fact  gives  a 
slightly  larger  gain  constant  also.  Note  that  as  with  most  su.h  opera¬ 
tional  amplifier  circuits  we  desire  infinite  gain,  K  =  <»,  in  the  basic 
amplifier  itself,  in  which  case  the  grounded  amplifier  configuration  gives 


Vo  =  RpC  Vi 


( 1 1 1~4  ) 


Concerning  summation,  the  diagram  of  Fig.  III-5  yields  a  convenient 
circuit  for  integration  which  has,  for  K  =  °°,  the  input-output  relation- 
ship  [l,  p.  ] 


v 

o 


m 

+ 

j=i 


RG 


+  1  r+  + 

-  G  .  v  . 

,+  J  J 


m 


2  "Vj 

j=l 


(lll-5a) 


Through  this  relationship  any  values  of  the  coefficients  can  be  obtained 
through  a  solution  of  simultaneous  equations  since,  for  the  resistance 
R*  and  R.  we  have  the  necessary  conductances  defined  as 


G 


+ 


G*  =  "7  ,  g”  =  — -  (lll-5b) 

J  R.  J  R. 
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Fig.  111-5.  SUMMER. 


However,  it  should  be  observed  that  inconvenient  values  for  construction 
through  the  use  of  integrated  circuits  can  occur  and  thus  a  cascade  of 
components  may  sometimes  be  necessary. 


D .  Scalar  Degree  Two  Realizations 

The  most  practically  met  situations  are  those  of  scalar  transfer 
functions.  In  such  cases  the  transfer  function  can  be  written  as  the 
product  of  oeg.*ee  one  and  two  factors,  having  real  coefficients  if  we 
assume  that  the  original  transfer  function  is  rational  with  real  coeffi¬ 
cients.  For  sensitivity  reasons  it  is  most  useful  to  construct  the  trans¬ 
fer  function  through  its  factors  instead  of  in  one  complete  form.  Thus 
we  exhibit  a  structure  for  the  transfer  function 

Sp  +  ci 

T  (p )  =  d  +  - - 2  ( III-6a ) 

p  +  2^0)  p  +  Ob 


where  we  assume  for  stability  reasons  that  the  undamped  natural  frequency 

CD  and  the  damping  ratio  £  are  nonnegat  ive .  We  ’’omark  that  the  qual.it  v 
n 

factor  Q  can  be  defined  by 


(lIl-Gb'' 


SEL-67-1 1 0 


48 


and  that  degree  one  transfer  fund  ions  are  simply  realized  (and  hence 
loft  as  an  exercise). 

We  claim  that  a  realization  of  the  general  degree  two  transfer  func¬ 
tion  is  given  by  [application  of  Eq.  (l-lOc)  gives  Eq.  (lll-6a)  I 


0 

1 

f 

II 

CDf 

0 

2 

1 

a 

-Uj 

n 

l 

.  j 

C  =  [ci  c2j  ,  D  =  d 

(lll-6c) 


Assuming  nonnegative  c^ ,  c^,  d,  a  circuit  diagram  suitable  for  integra¬ 
tion  would  be  as  shown  in  Fig.  III-6  where  the  values  of  resistance  can 
be  adjusted  for  available  capacitance  ranges.  The  presence  of  feedback 
loops  can  readily  be  seen,  as  well  as  an  appreciation  gained  for  the 
complications  attendant  on  going  to  the  complete  simulation  of  higher 
degree  transfer  functions  (without  the  initial  factorization).  We  observe 
that  the  minimum  number  of  capacitors,  two,  is  used  for  Fig.  III-6. 


Fig.  II1-6 .  POSSIBLE  DEGREE  TWO  SCALAR  SIMULATION. 
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Canonical  Equation  Simulation  -  Admittances 


In  Section  Il-b  we  saw  that  the  state-variable  equations  could  be 
set  up  for  admittance  "transfer"  function  (matrices)  by  extracting 
capacitors.  Here  we  can  reverse  the  procedure.  Thus  consider  the  resis¬ 
tive  (n+c)-port,  assumed  time-invariant,  described  by 

ill  iXl2 

Xe  =  {lII";a) 

Jai  1*22  _ 

and  loaded  in  its  final  c-ports  by  c  unit  capacitors,  aa  shown  in 
Fig.  III-7.  We  calc  late  for  the  input  admittance 

il„  -  il2(pic  +  i22)'i21  (lII-7b! 

If  we  compare  the  result  with  that  of  the  transfer  function 

T(p)  =  g  +£(plk  -  A)_1B  (i-lOc) 

we  see  that  the  identification 

To  -Cl 

■  u  j 


is  possible,  with  the  dimension  of  the  state  chosen  as  the  number  of 
capacitors,  k  =  c.  Consequently,  given  a  minimal  (or  even  nonminimal) 


Fig.  II 1-7.  CAPACITOR  LOADED  STRUCTURE . 
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real  izatio”  R  =  tA,jg,£,J>),  we  can  construct  a  circuit,  when  the  t  runs 
fer  function  is  an  admittance  matrix,  by  synthesizing  of  Eq.  \  1 1 1—7*.  i 

and  loading  in  k  =  c  unit  capacitors.  But  Jf  being  a  constant  matrix 
is  realized  through  the  use  of  (positive  and  negative)  resistors  and 
gyrators.  Later  we  will  show  how  can  be  transformed  to  become  posi¬ 

tive-real,  if  the  original  transfer  function  admittance,  jf  ,  is  positive 
real  but  such  requires  the  development  of  more  theory.  We  can  remark,  lir-v 
ever,  that  if  the  state-variable  equations  have  a  term  E  u  added  to  thf 
output  equations,  this  term  can  be  synthesized  by  a  transformer  network 
(constructed  from  gyrators  if  desired)  loaded  in  unit  capacitors  with  the 
result  connected  in  parallel  with  that  of  Fig.  III-7. 

To  synthesize  Y  Itself,  we  can  proceed  by  decomposing  it  into  its 

"c 

symmetric  and  skew-symmetric  parts, 

r»c  *c  ay  +  Xc  sk  (iII-Ho 


where 


2Y  =Y+Y,  2Y  =  Y  -  Y  (lll-8l>) 

""c  sy  ■"•c  »c  ’■c  sk  ’■c  -c 

and  again,  the  super  tilde  denotes  transposition.  The  skew-symmetric 

part  is  immediately  constructed  from  gyrators,  one  for  each  nonzero  entrv 
for  example.  The  symmetric  part  can  be  further  decomposed  is 

Y  =  G[l_  +  (-1  )]G  (ill -Be; 

where  +  denotes  the  direct  sum  of  two  matrices.  The  righi  side  ot 
Eq.  (lll-8c)  can  be  synthesized  by  loading  a  gyrator  coupling  network  of 
admittance  matrix 


0 


G 


Y 

*g 


(in  t;.i 


r>] 


:: r,i, -  f>  j- 1  m 


in  r+  unit  positive  resistors  and  r_  unit  negative  resistors,  as 
shown  in  Fig.  III-8  [recall  that  a  formula  similar  to  Eq.  (lll-7b) 
appli^r].  The  coupling  structure  itself  results  as  a  parallel  connection 
of  the  circuits  for  the  symmetric  and  skew-symmetric  parts  of  Y  . 


As  an  example  of  the  method,  let  us  confide?  the  degree  two  low-pass 
admittance 


yin(|,)  ' 


+  2£p  +  1 


(lll-9a) 


We  observe  that  this  admittance  is  not  positive  real  (as  l/y  has  a 

in 

double  pole  at  infinity)  in  which  case  active  devices  must  be  incorporated. 
Combining  the  realization  of  Eqs  (lll-6c)  with  Yfi  of  Eq.  (lll-7c) 
yields 


Y 

n»c 


-1  0 
0  -1 


1  2CJ 


(lll-9b) 


which  has  the  symmetric  and  skew-symmetric  parts 


■  0 

-i/2 

-1/2 

0 

1/2" 

0 

II 

x: 

(0 

u 

>* 

'  0  *l/2 

1/2  0 

-1/2" 

-1 

.  1/2 

0 

2C- 

_l/2  1 

0  . 
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To  diagonalize  ihe  symmetric  part  we  can  add  -\/4rt  times  the  last  row 
to  :he  first  a:  d  then  add  -4?  times  the  first  row  to  the  second.  In 
terms  of  elementary  matrices  this  gives 


p 

0 

o' 

p 

0 

-l/< 

1 

0 

0 

1 

0 

- 

Lo 

c 

1. 

p 

0 

1  . 

1  o  0  l  o  I  -  |  o  ?r.  o 

o  lj  Lo  o  lj  [  o  o  21 


-l  /«'  o  n 


On  multiplying  out  the  inverses  of  the  transformation  matrices  (which 
are  easily  found  by  changing  sign  on  the  off-diagonal  terms),  we  arrive 
at 


-l/8£  0  0  "I  T 1.  4C,  0 

0  2£  0  0  10 

0  0  2(J  [l/4C  0  1 


( Ill-9d ) 


We  observe  that  the  diagonal  matrix  is  not  quite  in  the  form  used  in 
Eq.  (lll-8c),  but  this  is  not  crucial  since  we  merely  use  nonunit  resis¬ 
tors  with  the  negative  one  placed  first  (the  other  form  can  easily  be 
obtained  by  using  some  additional  steps).  We  then  wish  to  load  a  pyrn- 
tor  6-port  described  by 
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in  one  negative  and  two  positive  resistors  to  obtain 
is  shown  in  Fig,  III -9a. 


Y 

mc  sy 


The  result 


J  . 


(a) 


(b) 


Fig.  III-9 ,  CIRCUITS  FOR  Y  (a)  AND  Y  ,  (b). 

"^c  sy  mc  sk 


The  circuit  for  ^  is  similarly  obtained  and  is  shown  in  the 

C  S  K 

(b)  portion  of  Fig.  III-9.  The  two  portions  of  this  figure  are  connected 
in  parallel  with  the  final  two  ports  loaded  in  unit  capacitors  to  obtain 
the  desired  input  admittance  at  port  1. 


F.  Discussion 

Using  the  canonical  state-variable  equations,  analog  configurations 
can  easily  be  set  up  using  a  block  diagram  representation  of  the  equa¬ 
tions;  the  resulting  components  are  realized  through  summevs  and  inte¬ 
grators,  the  latter  being  obtained  through  the  use  of  operational  amplifier 
circuits.  Since  integrated  operational  amplifiers  have  proven  extremely 
practical  and  since  the  only  other  elements  needed  are  resistors  and  ca¬ 
pacitors,  both  of  which  can  be  integrated,  the  method  is  quite  useful 
for  integrated  circuit  designs. 

It  is  of  interest  to  observe  that  exactly  as  many  capacitors  are 
used  as  there  are  state  variables,  and  in  fact  no  fewer  can  ever  be  used. 
Since,  of  the  components  required  here,  capacitors  are  the  most  difficult 
elements  to  make  in  integrated  circuits,  the  method  is  about  as  conve¬ 
nient  as  could  ever  be  hoped  for.  As  a  consequence  we  have  introduced 
some  basic  configurations  particularly  suited  for  integration. 
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It  should  be  mentioned  that  In  integrated  circuits  the  ratios  of 
resistors  are  rather  accurately  obtained,  wn  .vroas  absolute  values  are 
extremely  hi'rd  to  fix  accurately.  If  we  observe  the  coefficients  for  the 
summer  multipliers,  Eq.  (lll-5a),  we  see  that  indeed  these  coefficients 
depend  only  upon  ratios  of  resistances.  The  situation  is  somewhat  dif¬ 
ferent  for  the  integrator  where  both  resistance  and  capacitance  are 
involved.  In  fact,  since  RC  products  only  quite  far  away  from  unity  are 
available  in  integrated  circuit  form,  it  is  important  to  introduce  an 
integrating  scale  constant  in  the  state-variable  equations,  the  X  of 
Eq.  (lll-2b) . 

We  observe  that  although  the  equations  allow  time-variable  coeffi¬ 
cients  and,  in  fact,  the  circuit  representations  hold  for  such  coeffi¬ 
cients,  it  is  practically  quite  difficult  to  perform  time  variations  on 
the  operational  amplifier  structures. 

Although  we  have  not  discussed  the  possibility,  it  is  actually  more 
convenient  to  perform  time  variations  by  use  of  the  capacitance  extrac¬ 
tion  method.  But  we  have  discussed  how  the  previous  analysis  method, 
through  capacitor  extraction,  can  be  carried  over  to  synthesis  to  create 
a  resistive  coupling  structure  by  specifying  the  admittance  coupling 
matrix  Y£  in  terms  of  the  realization  R  =  •  In  conjunction 

with  this  we  have  given  one  method  of  synthesis  of  in  terms  of  gyra- 

tors,  which  can  be  integrated  [2],  and  positive  and  negative  resistors. 
Since  the  negative  resistors  cause  some  concern  for  practical  integration, 
it  is  of  perhaps  more  practical  interest  to  point  out  that  can  be 

obtained  as  "<n  interconnection  of  voltage-controlled  current  sources  and 
that  such  sources  are  relatively  easy  to  integrate  [l,  p.  ]. 

Of  the  two  methods  presented,  the  first  probably  has  the  advantage 
in  scalar  situations  of  allowing  for  smaller  sensitivities.  To  obtain 
these  sensitivities  of  small  size  it  is  important  to  decompose  the  trans¬ 
fer  function  into  degree  one  or  two  portions  and  cascade  the  resulting 
sections.  However  it  is  worth  mentioning  that  a  good  sensitivity  analysis 
of  the  second  (capacitor  extraction)  method  has  as  yet  not  been  made. 

Here  we  really  only  treated  the  synthesis  of  voltage  transfer  func¬ 
tions  (by  the  operational  amplifier  techniques)  or  of  admittance  matrices 
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(by  the  capacitor  extraction  methods).  However,  by  the  use  of  voltage- 
to-current  or  current  converters,  other  specifications  can  equally  well 
be  realized. 
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H.  Exercises 

1.  Set  up  the  state-variable  configuration  using  integrated  opera¬ 
tional  amplifiers  for  the  degree  one  transfer  function.  Compare 
with  results  obtainable  with  simple  RC  circuits. 

2.  Discuss  modifications  needed  in  the  theories  if  terms  of  the 
form  E  u  are  present.  Explain  why  these  are  avoided,  where 
possible,  in  the  operational  amplifier  techniques. 

3.  Complete  the  example  of  Section  E  by  drawing  the  final  overall 
circuit.  Compare  with  alternate  methods  and  discuss  advantages 
and  disadvantages  of  the  method. 

4.  Synthesize 


s'in(|,)  *  ~ 


12. 


p  +  2£p  +  1 


5.  Discuss  circuits  fc-r  obtaining  the  gyrators  and  negative  resis¬ 
tors  needed  in  Fig,  III-9. 

6.  Investigate  methods  of  obtaining  practical  realizations  for  the 
case  of  time-variable  structures. 
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* 
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soupjonner ,  et  cependant  leurs  yeux  ont 
une  certitude  si  profonde  qu'll  faut  qu'ils 
sachent  tout  et  qu'ils  a  lent  eu  plus  d'un 
scir  le  temps  de  se  dire  leur  secret. 

M.  Maeterlinck 

"Les  Avert  Is"  du  "Tr4sor  des  Humbles" 


IV.  MINIMAL  REALIZATION  CREATION 


A.  Summary 

By  conversion  of  a  high  order  differential  equation  to  a  set  of  first 
degree  ones  a  minimal  realization  is  relatively  easily  obtained  in  the 
scalar  case.  For  matrix  transfer  functions  the  algebraic  method  of  Ho 
is  presented  for  obtaining  minimal  realizations. 

B.  Scalar  Minimal  Realizations 

Previously  we  have  seen  how  a  given  circuit  can  be  analyzed  to  ob¬ 
tain  an  appropriate  set  of  canonical  equations.  Likewise  we  have  seen  how 
a  circuit  can  be  obtained  when  a  realization  is  on  hand,  that  is  when  the 
canonical  equations  are  on  hand.  Here  we  complete  the  picture  for  time- 
invariant  structures  by  giving  an  algebraic  procedure  for  finding  a  min¬ 
imal  realization  from  a  given  transfer  function.  We  begin  with  the  scalar 
case  for  which  the  result  can  be  easily  given. 

We  therefore  first  begin  by  assuming  as  given  the  scalar  transfer 
function 


T(P) 


.  6  .  6-1  , 
dp  +  d^p  +  . . .  +  d2p  +  dx 


6  6-1 

p  +  a^.  p  +  . . .  +  a2p  +  a^ 


(tv-1) 


of  degree  &•  If  we  treat  p  as  the  differential  operator  d/dt  this 
transfer  function  defines  the  differential  equation 
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[pb  +  appb  1  +  ...  +  a2p  +  a1]y  =  [dp°  *  dg  p6  1  +  . . .  +  dgp  +  d^u 

( IV-2a) 

We  can  now  introduce  some  changes  of  variables  beginning  ’vith 

y  =  sT  f  du  ( IV-2b) 

which  results  in 

[p°_1  +  afipU~2  +  ...  +  a2^Ps1  +  a1s1  =  [(db-aPjd)p0-1  +  ♦  (d^a^ju 

Next  letting 

P>1  =  s2  +  (dg~a^d)  (lV-2c) 

results  in 

r  5-2  6-3  -1 

Lp  +  a6p  +  • a3JP«2  +  a2a2  +  al3l  = 

0(d6-ra6-ld)  '  a5(d5“a5d)^p8"2  +  +  Kdra!d>  "  a2(d6'a6d)^U 

Continuing  by  letting 

Ps2  =  S3  +  J(d5-rVld)  '  a5(Va5d)^  (lV_2d) 

etc.,  results  in  the  final  equation 

ps6  +  a6s6  +  a6_186_1  +  ...  +  a2s2  +  alSl  =  b&u  (lV-2e) 

where  b_  is  a  combination  of  the  a,  and  d.  coefficients.  We  have 
o  ll 

then  obtained  the  canonical  equations  which  can  be  summarized  as 
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We  observe  that  the  realization  is  minimal  since  A  is  &  x  6  and 

H 

T(p)  has  degree  &.  Also,  the  same  procedure  holds  for  the  time-varying 
case  with  these  however  being  additional  derivatives  of  coefficients  in 
the  B  matrix. 

From  Eq.  (lV-3)  many  other  (in  fact  all)  minimal  realizations  can  be 
obtained  by  use  of  nonsingular  transformations  on  the  state,  that  is  by 
introducing 


3  =  T  s 


59 


SEL-67-110 


C.  Matrix  Minimal  Realizations 


The  matrix  case  is  much  more  difficult  to  pursue.  We  folio*  the 
algebraic  procedure  of  Ho  [l]  by  first  introducing  a  nonminimal  realiza¬ 
tion  which  is  reduced  to  be  minimal. 

We  begin  by  observing  the  form  of  the  transfer  function  matrix  in 
terms  of  the  realization  matrices.  Assuming  the  realization  to  be  min¬ 
imal,  that  is  the  state  of  minimal  dimension  5,  k  =  6,  we  obtain  on 
expanding  the  inverse  of  pi  -  A  about  p  =  oo, 

i*  o  *** 

T(p)  -  D  ♦  £ 


where  A*  is  the  ith 
T(p )  itself  about  p 


Since  .T(p)  is  rational  we  can  equate  term  by  term  in  the  last  two  ex¬ 
pressions  to  obtain  that  R  =  {A.g.fi.gl  is  a  realization  if  and  only  if 

£=!_!=  X(°°)  ( IV-4c) 

£l  =  S  £l£  .  1  =  0.  1.  •••  ( IV-4d ) 

Our  job  is  to  hunt  for  an  ^A,  B,C  which  satisfy,  Eq .  (lV-4d);  we  com¬ 
ment  that  this  last  equation  holds  no  matter  if  the  realization  is 
minimal  or  not,  but  that  we  are  actually  searching  for  a  minimal  one. 

Since  T(p)  is  rational  there  is  a  relationship  among  the  A  of 
Eq.  ( IV-4b) .  To  obtain  this  relationship  we  can  find  the  least  common 
denominator  polynomial 

g(p)  =  Pr  +  arpr_1  +  ( IV-5a) 


,  n  "  i+1  ■» 

i=0  p 


<IV-4a) 


power  of  A.  By  making  a  direct  expansion  of 
<»  yields  the  coefficients  A^  for  the  series 


®  A 

Kp)  -  4.,  *  1  -7TT 

i=0  p 


( IV-4b) 
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of  T(p)  which  next  allows  us  to  write  the  transfer  function  as  a  matrix 
polynomial  divided  by  g(p).  Thus 


J(P) 


P\+l  +  vSr  +  ‘  ‘  +  PS2  +  Hi 

r  r-1 

p  +  afp  +  ...  +  ax 


( IV-5fc) 


As  a  consequence  the  product  g(p)£(p)  is  polynomial  and  on  using  the 
expansion  of  Eq .  (lV-4b)  we  have 


Equating  those  coefficients, 


of 


-l 

P 


to  zero  we  find 


""k  aj^k-r+j-l  ’  k  -  ( IV-5c) 

As  we  saw  in  Eq.  (lV-3)  the  A  matrix  was  the  companion  matrix 
determined  solely  by  g(p).  As  a  consequence  we  introduce  its  general¬ 
ization,  for  which  we  recall  that  T(p)  is  an  n  X  m  matrix.  Thus,  the 
generalized  (rn  X  rn)  companion  matrix  for  g(p)  is  defined  by 


( IV-6a) 


where,  as  before,  1  is  the  n  ■■  n  identity  matrix.  To  accompany  this 
m  n 

companion  matrix  we  nerd  the  generalized  ( rr.  nr.)  Hankel  matrix. 
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A 


r-1 


■®r 


“2 


A 

—  r 


A 


2r-2 


(lV-6b) 


From  Eq.  (lV-5c)  we  observe  that  acts  to  shift  rows,  or  columns,  of 

when  the  two  are  multiplied,  that  is 


A 


2 


A 

*~r 


n  s 

—  n-r 


=  Sfi 

~  r-  m 


4 


4 


A 


r+1 


A  , 
’~2r-l 


( IV-6c) 


where  the  superscript  tilde  denotes  matrix  transposition.  As  a  conse¬ 
quence  premul tiplication  of  by  J3*  brings  A^  to  the  (l,l)  position 

of  the  result.  In  order  to  isolate  this  position  we  define  the  p  x  7 
matrix 


1 

~P>7 


i 

i 

I 


( IV-6d) 


for  which  the  first  p  columns  are  the  identity  matrix  with  the  remain¬ 
ing  columns  zero.  Then 


a.  =  i  (nis  )  1 

“i  ~n ,  rn  -n-r  -m.rm 


( IV-6e) 


A  possible  realization  is 


B 


=  S  1 
—  r  — m,rm 


C 


1 

«-n,rn 


P  =  l(ro) 


( IV- 7 ) 
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for  note  that  Eq.  (lV-6e)  is  just  A.  =  C  A^B  which  is  as  required  by 
Eq.  (lV-4d).  This  realization  however  is  not  generally  minimal,  having 
k  =  rn  which  is  generally  larger  than  the  minimum  size,  &,  required; 
as  we  will  see,  Eq .  (lV-13),  this  latter  is  given  by  &  =  rank  S^.  As 
a  consequence  let 


£>  =  rank  S 
-  r 


( IV**8a) 


in  which  case  one  can  readily  find  nonsingular  matrices  M  and  N  to 
bring  _Sr  to  diagonal  form 


M  S  N  =  1R  1- 
—  T  *—  -&,rn  ~o,im 


( IV-8b) 


In  terms  of  the  matrices  defined  to  this  point  we  can  now  exhibit  a  min¬ 
imal  realization.  Our  result  is:  a  rational  n  X  m  transfer  function 
matrix  T(p),  finite  at  infinity,  has  a  minimal  realization  given  by 


acts  as  a  pseudo-inverse  for  since  direct  calculation  gives 


s  =  s  s"  s 

-r  —  r  t  _  r 


s  s# 

"  r  ■”  r  “  r  "*r 


( IV-lOb) 


Next  consider  the  following  sequence  of  operations  which  begins  from 
Eq.  ( IV-6e) . 
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a  =  i  (n1  s  )  i  =1  n1  s  s  I 

—  l  ->n , rn  — n  — r  -m , rm  ‘-n ,  rn  -n  -  r  -  r  -r  -m.rtn 

i  ti  U  *  u 

—  1  g  A  g '  <  SI  —  1  S  S^  S  q  A  gTT  g  j 

' n ,rn  ~r  -m  ~r  -r  -m,rm  -  -n,rn  -r  -  r  -  r  “‘m  -r  ~r  -m , rm 

=  (l  S  N  I-  )(l-  M  S  n1  N  1R  )(l_  MSI  ) 

~n,rn  t  -  -o,rm  -o , rn  -  -  r  -m  ~  -o,rm  -o,rn  -  —  r  —  m , rm 

=  C  (l.  M  n  S  N  1~  )*  B  =  C  AiB 

-  -o,m  -  ~n  -r  —  ~o,rm  -  —  —  — 

Here  the  next  to  the  lest  step  is  justified  by  iteration  of  the  result 


rn  M  Sr  ^  N  =  1R  M  &  S  S  ft  N  L 

~o,rn  • - r— m-  —  o,rm  -6  ,rn  ~  n  —  r  — r  ~  r  *~m  -  -  6 , 

=  lt  M  n  S  N  1 


rm 


If.  M  n  S  N  1-  1_  M  n  S  N  lc 

'■o.rn  —  ~n  r  —  —  &,rm~&,rn--— n~r——6, 


rm 


As  a  consequence  a  realization  has  been  obtained  and  it  only  remains  to 
show  that  it  is  minimal. 

For  this  latter  demonstration  let  us  introduce  the  ordinary  observa¬ 
bility  and  controllability  matrices 


P  —  [C i  A  C ,  . . . ,  A  C]  ,  Q  —  [ B ,  A  B ,  . . . ,  A  B]  (iV* 1 1 ) 


r-1 


Then  for  any  realization,  since  =  C  we  find  from  direct  multi¬ 

plication  that 

_Sr  =P  Q  (IV-12) 

Now  suppose  that  there  exists  a  realization  having  A  of  size  k  X  k 

with  k  <  6  =  rank  S  .  We  have  a  contradiction  since 
—  r 

rank  S  S  min  [rank  P,  rank  Q]  5  k  <  6  =  rank  S  (lV-13) 

—  r  —  *—  — -r 

where  the  middle  inequality  follows  from  P  and  Q  being  of  sizes 
k  X  rm.  We  conclude  that  the  realization  is  the  smallest  possible  with 
6  being  what  we  have  previously  called  the  degree. 
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0.  Examples 


Consider  the  transfer  function 


T(p) 


1 

P  +  1 


1 

(p+l)(p+2) 


( IV-l^a ) 


One  procedure  would  be  to  connect  a  degree  one  realization  between  the 
inpul  and  first  output  and  a  degree  two  realization  between  the  input  and 
seco”  i  output.  However  the  iinal  re.  ull  would  have  a  3-dimensional  state, 
which  would  not  be  minimal  since,  as  we  next  show,  two  dimensions  suffice. 
Hence  we  proceed  to  apply  the  theory  of  the  previous  section. 

The  least  common  denominator  is 

ft(p)  =  p2  +  3p  +  2  =  p2  +  a2p  +  aL  (lV-14b) 


Thus  we  have 


m  =  1 


n  =  2 


r  =  2 


( IV-14c) 


and  for  we  must  calculate  the  expansion  of  T(p)  about  infinity  up 
to  Ag .  We  find  by  simply  dividing  the  denominators  into  the  numerators 
beginning  with  the  highest  powers  of  p 


i  rji  i  r-n  i  r 1 

i(p)  ■  5  „  *  T  ,  *  —  ,  *  ••• 

H  L°J  p  L  ij  p  L-3  J 


i  .  i  i 

■*  A  ,  +  ~  A  +  — r  A.  +  —z  A_  + 
*1-1  p  »i0  2  3  »"2 

P  P 


(lV-14d) 


The  Hankol  matrix  can  then  be  formed 


Sir- 


j"*0  «ll 

[il 


a  ta 

[:]  u 


(lV-14e) 
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and  one  finds  by  the  use  of  elementary  operations  Ilia! 


The  final  matrix  necessary  for  Eqs.  (lV-9)  is  the  companion  matrix  asso¬ 
ciated  with  gfp) . 


( IV-14h) 


We  can  then  calculate  the  minimal  realization  using  Eq .  (lV-9) 


0 

0 

l 

0 


( IV- 15a  ) 


( I V  - 1 5b  ) 


( IV- 15c 
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(lV-15d) 


One  can  easily  check  that  T(p)  results  from  this  realization  through  the 
calculation  of  D  +  c(p  1  -  A)  ^B. 

By  physically  constructing,  as  in  Chapter  III,  the  canonical  state 
variable  equations 


( IV-15e) 


( IV-15f ) 


one  can  obtain  a  device  with  the  given  transfer  function  and  which  uses 
the  minimum  number  of  dynamical  elements  (capacitors,  say,  for  integrated 
circuits) . 

Next  let  us  consider  the  minimal  realization  of  the  general  degree 
two  scalar 


T(p)  =  d  + 


2 

P 


C2P  ♦ 

or  2 

+  2&D  p  +  ci) 
nr  n 


(lll-6a) 


which  was  previously  considered  (Fig.  III-6).  We  have 


m  =  n  =  1  ,  r  =  6  =  2 


(lV-16a) 


and  Eqs.  (lV-7)  already  give  a  minimal  realization,  as  do  Eqs.  (lV-3)  as 
well  as  Eqs.  (lll-6c).  For  Eqs.  (lV-7)  we  have 


A 


0  1 

-CD^  *  2£co 
n  ^  n 


(lV-16b) 


which  follows  on  identification  of  terms  from 

2  2  2 

g(p)  =  p  +  2£o)nP  +  CU^  =  p  +  a2p  +  ax  (lV-16c) 
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Likewi se 


S  - 
—  r 


c  -24.  •  c 

I  n  2 


24u.  c,  -  i  1  t-14' 
n  1 


it  2 


( IV- Hid) 


which  follows  irom  I he  expansion  of  T(p)  about  infinity 


t(P  )  =  ci  ♦  y 


20>.  e  -  (lolC2)  '2c 
n  1  n  2 

:i 

P 


( IV-lftc) 


Equations  (lV-7^  give 


0 

i 

C2 

A  =  0  = 

-  -1 

2 

-ti. 

n 

-24 

n 

,  B  =  S  T 

-  r  -  in ,  rrn 

e.  -  24..  c 

1  •  n  2 

('  -  1  -  il  il|  ,  D  -  ti 

—  —  n ,  rn  — 


(iv-lfil ) 


We  observe  that  the  calculations  for  Eq .  (lV-9)  are  sometimes  unneces¬ 
sarily  burdensome,  as  for  example,  in  this  case  M  and  N  are  not  even 
needed.  Also  from  the  simplicity  of  Eq ,  (lll-6e)  «'hich  has  C  -  I  c  ,  e  j 
we  see  that  perhaps  there  is  a  more  convenient  method  (as  yet  undiscovered) 
for  finding  minimal  realizations. 


E.  Discussion 

Using  a  basic  equation,  (lV-6e),  for  a  decomposition  of  the  matrices 

A  obtained  by  expanding  the  t vansf or  function  T(p)  about  infinity  a 
—  i 

generally  nonminlmal  realization,  Eq ,  (lV-7),  is  easily  found  from  which 
simple  but  ingenious  manipulations  lead  to  a  minimal  realization,  Eq . 
flV-fO,  The  matrix  ense  is  seen  to  lie  somewhat  a  generalization  of  the 
scalar  situation  where  a  minimal  realization  is  relatively  easily  obtained 
by  converting  a  higher  order  differential  equation  to  a  set  ol  first  order 
ones.  Because  the  method  proceeds  in  an  algebraic  manner  directly  lrom 
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the  transfer  function  it  is  quite  suitable  for  computer  synthesis  of 
systems,  although  as  yet  we  are  unaware  of  such  a  program  being  carried 
out.  In  fact  it  appears  that  it  is  worthwhile  looking  for  improved  meth¬ 
ods,  since,  as  the  last  example  has  shown,  there  are  sometimes  situations 
when  easier  calculations  than  those  called  for  by  the  general  theory  can 
be  used. 

There  are  of  course  other  methods  of  obtaining  minimal  realizations. 
One  such  is  to  augment  T  such  that  m  =  n,  make  appropriate  frequency 
shifts  and  constant  additions  such  that  it  is  positive  or  bounded-real, 
and  then  give  a  minimal  reactive  synthesis  of  the  result  [2j.  Other 
methods  exist  which  work  in  the  time  domain  from  impulse  response  ma¬ 
trices  £3  J  [4  ] .  But  for  the  time-invariant  case  the  procedure  of  Ho, 
given  here,  presents  the  most  promising  because  of  its  possibilities  for 
computer  synthesis  of  systems.  Nevertheless  we  will  later,  Chapter  IX, 
briefly  look  at  the  time-domain  for  time-variable  synthesis  procedures. 

At  this  point  we  have  on  hand  the  basic  portions  of  the  important 
theories.  We  have  seen  how  to  set  up  the  canonical  equations  from  a 
circuit,  and  now  from  a  transfer  function,  and  we  have  shown  how  to  ob¬ 
tain  a  circuit  from  the  canonical  equations  and  thus  from  a  transfer 
function.  As  a  consequence  our  remaining  topics  are  all  associated  with 
improvements  and  extensions  of  the  basic  results.  We  first  look  into  meth¬ 
ods  of  finding  equivalents,  which  require  more  knowledge  of  the  concepts 
of  observability  and  controllability. 
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G.  Exercises 


1.  Find  j)  in  compact  l'orm  lor  Eq .  (lV-3). 

2.  For  the  general  degree  two  scalar  transfer  function,  exhibit 
completely  Eq .  (lV-3)  and  compare  with  the  several  other 
results  available. 

3.  Insert  the  modifications  required  for  Eq .  (lV-3)  to  hold  for 
time-variable  circuits. 

4.  Find  a  realization  for 

1  1 
p  +  1  ’  (p+l)(pt2) 

and  compare  with  the  results  of  Eq .  (lV-15). 

5.  Find  a  realization  for 

T(p)  = 


for  an  arbitrary.  What  is  the  nature  of  the  result  when 
a  =  1? 

*6.  Investigate  the  realization  of  T(p)  by  factorization  into 

degree  one  or  two  parts  and  the  realization  in  minimal  form  of 
each  part . 

7.  For  Eq .  (lV-15d)  find  _M  and  _N  and  determine  a  minimal 

realization  using  the  general  theory  associated  with  this  M 

and  Compare  with  the  realization  of  Eq .  (lV-15f). 

8.  Find  a  realization  for  T(p)  =  l/(p+l)  and  one  for  T(p)  = 

l/( p  +l)(p+2)  and  "connect"  the  two  to  obtain  a  realization 

for  the  text  example  of  Eq .  (lV-14a).  Compare  the  result  with 
that  of  the  text  and  discuss  with  specific  reference  to 
minimality. 


P  ~  * 

p  +  1 
2 

p  +  a 
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II  passe,  entre  duex  etres  que  se  rencontrent 
pour  la  premiere  fois,  d'6tranges  secrets 
de  vie  et  de  mort;  et  bien  d'autres  secrets 
qui  n'ont  pas  encore  de  nom,  mats  qui 
s'emparent  imm^diatement  de  notre 
attitude,  de  nos  regards  et  de  notre 
visage. 

M.  Maeterlinck 

"Les  Avertis"  du  "Tresor  des  Humbles" 


CHAPTER  V 
EQUIVALENCE 


A.  Summary 

Through  the  use  of  various  transformations  on  the  canonical  state- 
variable  equations  one  can  generally  find  all  canonical  equation 
representations  for  a  given  transfer  function.  When  the  realizations 
are  minimal  this  occurs  through  nonsingular  transformations  on  the 
state.  When  it  is  a  question  of  nonminimal  equivalents,  decompositions 
involving  the  "encirclement"  of  controllable  and  observable  portions 
result. 

B.  Minimal  Equivalents 

Given  a  transfer  function  matrix  T(p)  which  is  rational  and  having 
T(“)  well  defined  we  have  seen  in  the  last  chapter  how  to  obtain  a 
canonical  set  of  state  variable  equations 

s  =  As  +  Bu 

y  =  Cs  +  Du 

»  •<*'  *»•* 

such  that  the  state  has  minimal  dimension,  6,  and  with 

T(p>  >  D  ♦  Cfplj-Ar'g  (I-lld) 


(I-lla) 

(I-llb) 
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One  problem  of  equivalence,  and  that  which  we  treat  here,  is  that  of 
finding  other  realizations  R  =  {A,  B,  C,  D) ,  perhaps  nonminimal,  for 
which  the  above  equations  are  Line.  Here  in  fact  we  will  find  all  such 
realizations.  However  we  first  show  how  to  find  all  minimal  realizations. 
Let  us  consider  as  on  hand  two  minimal  realizations  R  =  (A,  B,  C,  D| 

A  A  /■  ' 

and  U  zr  f  m  B,  C,  D|  of  a  tfiv<n  transfer  function  n>  in  matrix  T(p). 

•A  Ml  A* 

A 

IV"  define  the  observability,  P  aiul  P,  matrices  and  controllability, 

A 

Q  and  Q,  matrices  as  before,  Fq.  ( I V— 11);  then  v/e  find 


S  -  PQ 

=  PQ 

(V-la) 

■mm  2"  *»  M 

We  also  recall 

that 

P  and  Q  have 

o  rows  and 

are  of  rank 

b,  in 

which  case  QQ, 

**  * 

PP, 

and  the  same  expressions  in 

terms  of  £ 

VMA 

and  § 

are  b*b  nonsingular  mai rices.  If  we  premultiply  S  by  P  we  obtain 

Q  =  [(PP)_1P^]Q  =  TQ  (V-lb) 

M*  m  .tot 

which  serves  to  define  the  transformation  matrix  T  =  (PP)  *PP  which 
is  nonsingular  by  the  fact  that 

•v  A  A  A 

S  --  PTQ  =  PQ 
**  r  ■**  Mi  aa 


has  rank  b  and  T 
by  the  transpose  of 


QQ, 

IMr 


is  bxb.  Postmultiplying  both  sides  of  this  latter 

A 

Q  gives,  on  cancellation  of  the  nonsingular  matrix 


P  =  PT  (V-lc) 

mm  mm 


Since  the  first  m  columns  in  Q  are  B  we  conclude  from  Eq.  (V-lb)  that 

A  -v 

B  =  TB.  Likewise  since  the  first  n  rows  of  P  are  C  we  have  from 

(Ml  v**  A"* 

A 

Eq.  (V-lc)  that  C  =  CT.  The  canonical  state  variable  equations  are  then 

A*1* 


y 

<A*\ 


As  +  TBu 

/V  -] 

CT  s  +  Du 

AA  Avk-  M 


A 

s 


y 

*V*> 


A  A  A 

As  +  Bu 

M*  Am  M 


AA 

Cs  i-  Du 

Am,  rm.  •*  •’X 


(V-ld) 

(V-le) 
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It  Is  then  reasonable  that  the  Identification 


s  =  Ts  (V-lf) 

f»  M 

should  be  made.  In  which  case  T  ATs  =  As.  As  any  Initial  state  Is 
allowed  we  can  cancel  the  s  to  conclude  that  any  two  minimal  realiza- 
tions  are  related  through  a  nonsingular  transformation  by  the  relation¬ 
ships 


A  -1 

A  =  T  AT, 

9**  m  W.M 


B 

WK 


B 


> 


t 


=  CT 

m  mr 


(V-2a) 


In  other  words,  any  two  minimal  realizations  are  given  one  in  terms  of 
the  other  through  Eqs.  (V-2a)  where  in  fact 


T  =  (PP)_1PP 

to*  AM 


(  V-2b) 


By  letting  T  run  through  all  nonsingular  6x6  matrices  we  obtain  all 
minimal  realizations  from  a  given  one. 

We  comment  that  previously  we  checked,  at  Eq.  (I-lle),  that  this 
transformation,  Eq.  (V-2a),  does  leave  the  transfer  function  invariant. 

As  an  example  let  us  reconsider  the  Brune  section  of  Chapter  I  for 
which 


s  = 

MM 

0 

-Vc2l5  ♦[  Mi 

(I-9g) 

•  V0! 

0  J  °. 

y  = 

MM 

0  is+fl  ol  u 

m  Pm 

(l-9h) 

0 

(grg2)/c2J  1°  1- 

If  it  is  desired  to  have  a  skew-symmetric  £  matrix  we  can  set 


T  = 


“ll 
L  t21 


t12 
*22  - 


and  examine  the  set  of  (nonlinear  in  t^j)  entries  in 
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T  '.AT  such  that  the  result  is  skew-symmetric.  Wo  find  a  suitable  T 

**  m  .*■» 

IS 


Thus,  wo  find 


which  wo  know  yields  an  equivalent  structure  to  the  original  Brune 
section. 

C.  Controllability  and  Observability 

In  order  to  proceed  to  notminimal  equivalents  it  is  necessary  to 
introduce  the  concepts  of  controllability  and  observability  which  we 
have  already  seen  enter  into  t he  theory  of  equivalence  through  the 
matrices  P  and  (}. 

fm  ** 

To  be  somewhat  precise  we  say  that  an  initial  state  ^(t^)  is 

control lablc  if  there  exists  a  finite  time  t  and  an  input  u(t), 

t  <  t  <  t  ,  such  that  s  (t)  =  ()  for  t  >  t  ,  that  is,  such  that 

the  state  can  be  brought  to  zero  (which  is  also  the  origin  of  the  state 

space).  By  beginning  on  a  trajectory  of  a  controllable  state  starting 

at  t  we  sec  that  later  values  of  time  yield  controllable  initial 
o 

states  and  lienee  we  can  work  with  controllable'  states  is(t)  in  which 
case  we  can  decompose  t  lie  state  spare  into  the  so!  of  controllable  s(al«s 


71 


SEL-67-1 in 


and  those  which  are  not,  the  uncontrollable  states  (this  requires  also 
letting  t1  tend  to  infinity). 

On  the  other  hand  an  initial  state  s. (t  )  is  observable  if  there 

<wl  o  — -  - 

exists  a  finite  time  t  and  a  zero  input  output  y(t),  t  <  t  <  t  , 
such  that  a  knowledge  of  y(t)  determines  s  (t  ).  Again  we  extend 

|MA  O 

the  concept  to  all  times  and  hence  can  decompose  the  set  of  states  into 
those  which  are  observable  and  nonobservable. 

Unfortunately  the  background  concepts  needed  to  derive  useful 
results  from  these  definitions  are  rather  complicated  so  we  will  state 
ome  of  the  results  omitting  to  some  extent  noncrucial  proofs.  As  back¬ 
ground  we  recall  that  a  vector  jc  is  in  the  nullspace  of  a  matrix  ,M 
if  Mx  =  0.  Considering  the  time-invariant  case,  a  state  s(t  )  is 
controllable  if  it  is  not  in  the  null-space  of  [1,  p.  409] 


W(t  ,t  )  =  [  e~(to_t)BBe^(t_to)dt  (V-4a) 

**  o'  1  J  «*»■ 


Likewise  a  state  is  observable  if  it  is  not  in  the  null-space  of 


M(t  ,t  )  =  f  (A(t o  t)CCe«(t_to)dt  (V-4b) 

vm  O’  1  J  “** 


One  can  see  the  validity  of  this  latter,  for  example,  by  noting  that 
the  zero  input-output  is 

y(t)  =  Ce^*  *o^s(t  ) 

*  I**  O 

If  we  multiply  by  exp[A(to~t) ]C  and  integrate  we  have 


f  e^to  t^Cy(t)dt  =  M(t  ,t  )s(t  ) 
J  o  1  «  o 


from  which  s(t  )  can  be  determined  if  it  is  not  in  the  null-space  of 
o 

M(t  ,t  ). 

»  o  1 
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From  the  similarity  in  iorm  and  associated  statements  of  W  and  M 
we  see  that  the  controllability  and  observability  properties  of 

s  ^  As  h  Du  (I-lla) 

Art  t»*  frr 

y  =  Cs  i  Du  (I-llb) 

ff,  _  «••*  W*  — 

are  respectively  the  observability  and  controllability  properties  of  the 
transposed  system 


(V-5a) 


y,  -  Bx  +  Du. 

**  t.  W  t"  m  *  t 


(V-5b) 


This  result  is  customarily  referred  to  as  the  principle  of  system 
duality  and  esscntinllv  mean:-'  that  nc'd  to  consider  (>nly  oir  n|  Mr 
two  concepts  (controllability  or  observability)  as  independent. 

Actually  the  matrices  M  and  W  are  rather  difficult  to  work  with 

r*\  v** 

and  have  been  only  introduced  to  obtain  the  principle'  of  duality  which 
links  the  concepts,  Equivalent  results  are  expressed  in  terms  of  the 
observability  and  controllability  matrices 


P  =  [C,  AC, 

ft  AA 


Q  -  [B,  AD, 

rrn  M* 


Ak_1Bl  ( V-6) 


where  k  is  the  order  of  A.  Thus  the  set  of  eon  1 rollahle  (initial) 
states  is  the  space  spanned  by  the  columns  oi  ()  while  the  set  of  non- 
observable  states  is  perpendicular  to  the  space  spanned  by  tin  columns 
of  [2,  pp,  500,  50'1 1 .  These  criterion  are  easier  lo  apply,  as  com¬ 

pared  to  those  for  M  and  W.  We  note  that  if  P  and  Q  have  rank 

Ml  M  «■ 

k  then  all  states  are  controllable  mil  observable;  in  liiis  situation 

it  is  actually  time  that  the  realization  is  minimal,  k  =  o  (as  S 

•"T 

of  Eq.  (IV-12)  has  rank  k). 


I).  Nonminimal  Equivalenls 

Al  lhiz  point  we  can  turn  lo  the  general  result.  From  two  sect  ions 
previous  we  know  how  lo  find  all  minimal  equivalents  so  wi  are  interested 
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in  the  eases  where  the  dimension  k  of  the  state  is  larger  than  the 
minimum  size  6.  Such  can  occur  when  there  are  either  uncontrollable  or 
nonobservable  states  present,  or  both.  Consequently  it  is  convenient  to 
partition  the  state  vector  s  into  vr "ious  canonical  components,  as 


r ~cn  -co  ~un  ~uo, 
s  =  Is  ,  s  ,  s  ,  s  I 

**c  *  *«  •*  ’  — 


(V-6) 


where  the  superscript  indices  have  the  following  meaning: 


c: 

controllable 

• 

o: 

observable 

u: 

uncontrollable 

n: 

nonobservable 

Thus,  for  example, 

uo 

s 

MW 

is  the  set  of  uncontrollable  but  observable 

states. 

To  accompany  the  partition  of  the  states  we  can  partition  a  given 
realization  to  obtain  the  canonical  equations  in  the  form 


r  •  cn  -| 
s 

•  CO 
s 


•  un 
s 


•  uo 
s 


y 


■-\l 

-12 

A13 

~21 

—22 

-23 

A31 

-32 

^33 

-41 

-42 

A<J 

[£l  —2 

C„  C .]  r  s' 

-3  «-4  1  — 

CO 

s 

un 

s 

uo 

L  s 


A 

cn 

rs  1 

+ 

rB,i 

-'14 

-24 

*34 

■M 

CO 

s 

— 

un 

s 

— 

uo 

52 

53 

-44- 

_s 

L*4J 

+  Du 


(V-7a) 


(V-7b) 


In  order  to  have  the  state  s  partitioned  in  the  form  given  by 
Eq.  (V-6)  generally  requires  that  a  transformation  be  performed  upon 
the  state.  But  once  such  a  partition  has  been  performed  we  see  from 
the  physical  meaning  of  controllability  and  observability  that  B3,  B^, 

Cg  are  zero.  Also  since  there  should  be  no  way  for  the  input  to 
couple  to  the  uncontrollable  states,  A  ,  A  ,  A  and  A  are  also 
zero.  Since  also  the  nonobservable  states  should  not  be  seen  at  the  out¬ 
put  even  after  coupling  through  observable  states  we  find  A  ,  A  and 
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arc  also  zero. 


Thus  we  ran  obtain  111*'  decomposition 


-12  -l.T  -14 


4»s  4* 


-  cn  I  r  -1 
s  +  B,  u 

CO 

2  »2 

sUU  0 


(V-8a) 


r«  C2  o  Cj]  [scn  I  +  Du 


(V-8b) 


I'dii.i l  1  one  (V-8)  -iv  .t  canon  1  <  .1 1  i*iim  •nr  real  i/.n  1  inns  *>l  a  given  1  rior- 1  <•»• 
Lime t ion  T(p)  when  (lie  state  lias  nonminimal  size,  They  can  lie  obtain*  <1 
from  any  other  1  eal  i /.*  •  i  <>r  l»v  *  1  rnnslonr.at  ion  'I  applied  l<*  the  slat** 
s  p.  1 72 1 


S  r  T  S 

—  —  r—  <  * 


(  V-8c) 


To  aclnally  find  T^  llni*  are  fixed  proceuin.s.  bin  we  remark  Mint  the 

dimensions  of  the  four  subcomponen  I  s  of  s  can  be  de  t  *  -rained  from  1* 

and  Q  in  which  case  one  can  solve  lor  T  by  hunting;  for  a  canonical 

realization  I?  -  {A  ,  B  .  C  ,  1)1  ,  that  is,  one  of  .‘he  form  ol  Kqs. 
c  *-c  *"C  •»  c  — 

(V-8a,b),  in  terms  of  a  given  one  11  -  (A.  B.  C,  1)1  by  applying  the 
resul t 


T  A  =  AT  ,  TUB,  C  CT  (V-8d) 

«  C**C  -i  "-C'  ***>' 


which  holds  since  Eq.  (V-8c)  is  valid. 

Tin*  import nnt  point  to  *>bsei-ve  is  that  only  the  matrices  oi  a 

minimal  realization  1!  lb,  1)1  enter  into  T(p)  ,  that  is 

111  -22-2-2-  •*- 


T(p)  -  I)  +  C(pl-A)  ’ll  I)  1  (  (pi  -A  )  ~ 1 B  (V-'M 

ml  M  «•*  *  K  V  ~  JL  \  J 
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the  other  nonzero  entries  of  R  are  com,  letely  arbitrary.  Thus,  given 

c 

any  minimal  realization  we  can  find  all  other  realizations,  nonminimal 
or  not,  by  "encircling"  the  minimal  one  arbitrarily  but  as  required  by 
Eqs.  (V~8a,b)  and  then  transforming  by  arbitrary  (nonsJngular)  Tc  as 
required  by  Eq.  (V-8d).  This  being  the  case  we  can  derive  any  realization 
ft  from  any  other  R  as  shown  in  Fig.  V-l  [4]. 


R 


- >— 

b'.»sis 
r nange 
T 

"C 


n  v  n 

\ 

A  A 

c  m 

minimal 

basis 

m  c 

7 

basis 

extraction 

change 

encirclement 

change 

T 

* 

c 

Equivalence  for  Two  Realizations  R  and  ft 


Figure  V-l 


Of  most  practical  interest  to  us  is  the  derivation  of  nonminimal 
realizations  from  minimal  ones.  Since  we  can  readily  find  a  minimal 
realization  the  procedure  of  encirclement  is  convenient  for  taking  a 
given  transfer  function  T(p)  and  finding  all  realizations.  Note  that 
Eq.  (V-9)  shows  that  minimal  realizations  have  all  state  components 
controllable  and  observable. 

As  an  example,  the  circuit  of  Fig.  V-2  has 


s  =  -s+u  (V-lOa) 

y  =  s  (V-lOb) 


If  for  some  reason  one  were  to  want  a  configuration  using  two  capacitors, 
perhaps  to  be  used  jointly  for  some  other  purpose,  but  with  only  observable 
portions  one  could  proceed  from 
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Oin  tan  easily  check  that  these  two  sets  ol'  canonical  equations  viilt! 


the  same  transfer  function.  To  obtain  the  most  general  real i/a t  l>.  o 
the  type  required  one  next  can  apply  tin*  transl.in.iiil  imi  >n  hq. 


1 


Example  Circuit 
Figure  V-2 


In  E<|S.  (V-10)  we  comment  that  or,  ft,  y  are  arbitrary  constants. 
However,  If  y  -  0  t lien  sU°  is  not  observable  so  that  there  is  some 
further  constraint  placed  on  the  arbitrariness  of  the  matrix  C.  ;  this 

w-l 

constraint  we  believe  remains  to  be  determined  hut  should  he  ••xpressalile 
in  terms  of  the  observability  mutrix  P. 

From  Section  Ill-d)  we  know  that  for  u  =  and  y  =  i  the 

equations  of  Eq.  (V-10)  can  be  physically  realized  by  loading  a  circuit 
realization  of  the  coupling  admittance  matrix 


(V-lOe) 


in  two  unit  capacitors.  To  obtain  the  output  as  a  voltage  one  can  then 
insert  a  resistor  and  its  negative  in  series  with  the  source  to  convert 
y  =  i  to  y  =  v  ,  as  shown  in  Fig.  V-3.  Such  gives  an  alternate  but 
not  too  practical  realization  scheme. 


i 


Augmentation  to  Convert  to  Voltage  Output 
Figure  V-3 
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E. 


Discussion 


Given  one  set  of  canonical  state  variable  equations  we  have  shown 
in  this  chapter  how  to  find  all  others  possessing  the  same  transfer 
function.  Since  we  have  previously  seen  how  to  find  one  realization,  in 
fact  minimal,  from  a  transfer  function,  we  are  now  in  a  position  to 
find  all  state  variable  realizations  from  a  given  transfer  function. 

In  some  sense  then  we  have  found  all  equivalents. 

However  in  another  sense  we  have  not  completed  the  picture  since 
we  have  not  shown  how  to  find  all  physical  circuits  yielding  a  given  set 
of  canonical  equations.  To  be  sure  there  are  several  since,  for  example, 
we  can  give  an  analog  simulation  or  we  can  synthesize  a  resistive  coupl¬ 
ing  network  to  load  in  capacitors  and  indeed  these  two  methods  yield 
different  structures.  However,  one  can  apply  the  standard  theory  of 
Howitt  [5]  to  generally  find  all  physical  resistive  coupling  circuits, 
the  ones  containing  operational  amplifiers  usually  being  included  in 
the  result. 

The  theory  has  been  given  for  time-invariant  systems.  The  primary 
reason  for  excluding  time-variable  ones  at  this  point  is  that  one  can 
not  generally  expect  the  decomposition  of  the  state  into  the  components 
s  =  rsCn,  sC°,  s'n,  sU0l  to  hold  for  all  time  unless  there  is  some 
restriction  placed  upon  the  system.  Of  course  time-invariance  is  a 
sufficient  restriction  in  which  case  a  constant  transformation  exists  to 
bring  the  realization  into  canonical  form.  Nevertheless  much  can  be 
said  about  the  time- variable  case  where  the  use  of  proper  transformations, 
which  may  be  time-variable  even  in  the  time- invariant  case,  yields  a 
different  canonical  form  [6],  Perhaps  the  flow  pattern  of  Fig.  V-4  is  of 
interest  in  depicting  the  structure  of  the  actual  decomposition. 

The  somewhat  complete  nature  of  the  equivalence  results,  which  have 
not  been  obtained  by  other  means,  should  give  sufficient  justification 
for  the  existence  and  study  of  state  variable  theory.  Nevertheless  the 
concepts  of  controllability  and  observability  can  be  expressed  in  terms 
of  cancellations  in  [pL.-A]  and  C[pl.-A]  respectively  [1, 
pp.  389,  408],  Likewise,  if  internal  variables  are  considered  in  an 
(i v  =  Si  description  the  concepts  can  be  expressed  in  terms  of  the  (i_ 
and  matrices  [7]. 
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Flow  Pattern  for  Canonical  Realization 
Figure  V-4 

In  summary,  using  the  dual  concepts  of  controllability  and 
observability  we  have  been  able  to  obtain  a  feeling  for  the  internal 
structure  of  time-invariant  systems  through  the  form  of  canonical 
realizations.  Using  the  results  we  have  also  been  able  to  obtain  all 
canonical  state  variable  equations,  thus  allowing  a  designer  maximum 
freedom  of  choice  to  obtain  a  desired  circuit  configuration. 
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G.  Exercises 

1.  Complete  two  syntheses  of  the  canonical  equations  of  Eq.  (V-10) . 
Compare  the  results  and  discuss  relationships  between  them. 

2.  Find  all  canonical  equations  using  two  state  parameters  for 
equivalents  to  the  circuit  of  Fig.  V-2.  What  chants  if  an  arbitrary 
number  of  capacitors  are  allowed? 

3.  Suppose  that  it  is  possible  to  find  a  time-variable  transformation 

T  (t)  to  bring  the  state  to  the  canonical  form  of  Eq.  (V-6). 

•*c 

Discuss  the  changes  in  Eqs.  (V-8)  and  Fig.  V-l. 

4.  Discuss  why  the  basis  change  T  for  Fig.  V-l  could  actually  be 
omitted  from  the  figure. 

5.  Show  how  Tfi  can  be  created,  at  least  to  a  great  extent,  directly 
from  P  and  Q  [4,  p.  374], 

6.  Find  all  equivalents  for  the  integrator  of  Fig.  III-4a)  and  discuss 
factors  Influencing  the  choice  of  one  over  another. 
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Rein  n'est  visible  et  cependant  nous 
voyono  tout,  ils  ont  peur  de  noi'.s,  parce 
que  nous  les  avertissons  sans  cesse  et 
malgr4  nous;  et  a  peine  les  avons-nous 
abordls  qu'ils  sentent  que  nous 
rlagissons  contre  leur  avenir. 

M.  Maeterlinck 

"Les  Av«rtis"  du  "Tresor  des  Humbles" 


CHAPTER  VI 

SENSITIVITY  AND  TRANSITION  MATRICES 


A.  Summary 

Using  the  canonical  equations  transfer  function  sensitivity  can  be 
conveniently  expressed,  this  being  done  here  for  scalar  transfer  func¬ 
tions.  Time  domain  calculations  can  also  be  made  in  which  case  conve¬ 
nient  methods  of  computation  for  fundamental  matrices  are  also  presented. 


B.  Scalar  Transfer  Function  Sensitivity 

In  terms  of  precent  changes  it  is  of  interest  to  know  how  much  a 
transfer  function  changes  with  a  given  change  in  some  parameter.  Thus, 
for  transistor  circuits  it  is  of  interest  often  to  know  the  effects  of 
replacing  one  transistor  by  another  one  having  the  same  characteristics 
except  for  a  different  current  gain,  B.  Or  alternatively  with  integrated 
circuits  one  would  desire  to  know  how  the  overall  performance  is  affected 
by  a  change  in  temperature.  To  study  such,  the  sensitivity  of  a  (scalar) 
transfer  function  T(p)  to  a  parameter  x  has  been  defined  as  [l] 


T(p)  _  _x_  ^T(p) 
x  T(p)  dx 


(VI-1) 


Note  that  in  this  definition  the  sensitivity  is  a  complex  valued  function 
of  a  complex  variable  p.  In  most  cases  of  interest  one  really  desires 
i-i  know  the  behavior  of  the  magnitude  of  the  transfer  function  for 
sinusoidal  signals,  that  is  the  actually  desired  quantity  is  jcu)  | . 
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However  this  latter  is  analytically  difficult  to  work  with  and  one  does 
have  the  relationships 


,T(p)  _  s  |T(p)  |  T(p) 

x  -  x  J  r« 


(VI-2a) 


and 


,s^p)|  >  |SIT<P>I| 


( Vl-2b) 


both  of  which  are  relatively  easy  to  check. 

The  sensitivity  can  l>o  evaluated  in  terms  of  a  state  space  realiz- 
tion  through  lif fcrenliation  of 

T(p)  =  D  +  C(pl  ,-A)-1B  (VI-3a) 

If  for  any  matrix  G  we  realize  that 

=  -G~ 1  G_1  (VI-3b) 


then  we  obtain 


dT 


dc,  .  . . -1  ,  ..-1  3a.  ,  . -1 

^ +  T5'pik-i>  a +  a 

(VI-4) 

+  C(pl  -A)_1 
—  —  k  —  tx 


We  observe  that,  except  for  the  derivations,  the  only  operations  involved 
are  those  already  used  in  forming  the  transfer  function  from  the  realiza¬ 
tion.  Consequently,  this  method  of  determining  the  sensitivity  is  quite 
applicable  to  computer  analysis  of  circuits  where  we  have  previously 
seen  that  there  are  convenient  methods  of  obtaining  the  realization 
R  =  {A,  B,  C,  D}  from  the  circuit  diagram.  We  observe,  for  example, 
that  if  the  realization  is  set  up  in  the  special  form  of  Eq.  (IV-3) 
where  C  =  [1,  0,  ...,  0],  then  3C/3x  =  0  while  M/^x  also  takes  a 
simple  form  (having  only  nonzero  entries  in  the  last  row) . 
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where  p^  is  a  pole  of  the  transfer  function  T(p) .  In  general  the 
poles  of  T(p)  art  eigenvalues  ol  A^  or,  what  is  the  same,  zeros  of 
the  determinant  (pi  of  pi, -A.  If  we  assume  that  p  is  a  simple 

“  K  K  - 

eigenvalue  of  A  i hen  wc  can  evaluate  the  pole  sensitivity  s  ^  for 
p^  with  respect  to  x  as  follows.  We  have,  which  serves  to  define 
the  polynomial  K(p)  , 


^(p)  -  (p-Pk>K(p),  K(pk>  ^  0 


on  differentiation 


A(p) 

x 


(p-pk) 


K(p) 

x 


^K(p) 


(VI- 7a) 


( VI-7b) 


Solving  for  ep  /  x  on  letting  p  =  p  gives,  on  noting  that  K(p  ) 

K  K  K 

=  A(p)/ep  evaluated  at  p  =  p^, 


where 


pPh  _  Pk  _  _  A(p)/  x  - 

Px  “  x  r\(p)/-  p  I 

P  =  Pk 


( VT-7c) 


A(P)  =  det (pi  -A)  ( VT-7d) 

**  K  *" 

As  a  consequence  the  pole  position  sensitivity  is  relatively  easily 
evaluated  in  terms  of  the  A  matrix  and  with  the  us<  of  a  computet  |  2|. 

To  illustrate  the  situation  let  us  again  consider  the  transler 
function  of  Eq,  (VI-5n) ,  we  have 

\(p)  -  det  (pi  ^-A)  =  p2  f  2^,^p  t-  ^  (VT-8a) 

and  thus 


cV^(p)/-  *  =  2'; ^p ,  V\(p)/  P  =  2p  t  2'.  ( VI-81)) 

There  are  two  poles  of  T(p) ,  let  us  consider 


*7 


SEL-67-1 10 


(VI-8c) 


(VI-8d) 


D.  Time-Domain  Variations 

In  many  situations  the  quantity  of  most  importance  is  the  actual 
output  change  as  a  function  of  time  due  to  a  parameter  change.  In  such 
situations  the  canonical  state  variable  equations 

s  =  As  +  Bu  ( VI-9a) 

y  =  Cs  +  Du  (VI-9b) 

04 


can  advantageously  be  used. 

Again  let  us  consider  a  parameter  x,  as  well  as  constant  (in 
time)  realization  matrices  A,  B,  C,  D,  the  last  one  being  a  scalar  by 
virtue  of  our  treatment  of  single  input  single  output  systems.  Then  we 
find  on  differentiation  with  respect  to  x 


=  A(-U*)  +  [4»s  +  4^u] 

"  r  X  r'X 


(VI-10a) 


7x 


cC 

+  -*s  + 

r'X- 


pD 


(VI-10b) 


To  determine  r>y/$x  we  can  first  solve  Eq.  (VI-9a)  for  and  then 
Eq.  (Vl-lOa)  for  r's/.'x.  The  important  thing  to  observe  is  that  the 
same  matrix  A  occurs  in  the  two  situations,  only  the  forcing  functions 

differ  being  Bu  in  the  first  case  and  (.""W/  x)s  +  (  B/«'x)u  in  the 

^  *•*  •* 

second. 

The  problem  in  this  case  is  one  of  solving  the  differential  equation 
z  =  Az  +  f  with  f  known.  Such  solutions  are  obtained  in  a  straight- 
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forward  manner,  and  are  in  fact  conveniently  obtained  on  a  digital 
computer,  as  discussed  in  the  next  section.  Consequently,  1  he  variations 
in  the  output,  ey/Ax,  as  a  function  of  time  are  conveniently  obtain*  d. 
Of  course  they  can  also  be  normalized,  as  for  the  transfer  function,  to 
give  percent  changes  If  so  desired. 

E.  Transition  Matrix  Evaluation 

Theoretically  it  is  a  relatively  simple  matter  to  solve  the 
differential  equation 


z  =  Az  +  f  (Vi-iia) 

ia  M  m 

where  f  is  a  known  forcing  function  independent  ol  x  and  A  is  a 

on  A* 

square  kxk  matrix,  also  independent  of  x  but  perhaps  not  of  time. 

To  solve  Eq.  ( VI-1 la) ,  which  is  the  type  of  equation  appearing  in  Eqs. 
(VI-9a,  10a),  we  first  solve  the  equation 

Z  =  AZ,  Z(to)  «  lk  (Vl-llb) 

which  is  the  original  one  with  the  k- vector  x  replaced  bv  the  kvk 
matrix  Z,  without  the  forcing  function  and  with  the  identity  matrix  for 
initial  conditions.  The  solution  to  the  latter  equation  ran  be  denoted 
by  <J'(t,to>  and  is  called  the  transition  matrix  for  the  system.  In  the 
case  where  A  is  constant  in  time  this  transition  matrix  can  be  explic¬ 
itly  evaluated  as 


<T>(t,t  )  =  e-^  *o\  constant  A  (VT-12a) 

where  the  exponential  of  a  matrix  is  defined  precisely  by 

2  i 

e**  =  1  +  At  +  A2  +  . .  .  +  A1  x  , , .  (VI-12b) 

In  fact  one  can  directly  check  that  the  exponential  transition  matrix  of 
F,q.  (\T-12a)  does  solve  the  unforced  differentia]  equation  of  fVT-Ilh). 


89 


SE1/-87-]  1  0 


As  an  example,  if  as  in  Eq.  (V-3b)  we  have 


A  =  afo  -1 
1  0 

then 


-1 

°l» 

.3  3r  „ 

A  =  a  0 

W* 

l], 

A4  \ 

-  -  a  l%> 

A5  =  a5 A  (  VI- 13  b) 

0 

-1 

-1 

0 

(VI-13a) 


in  which  case  Eq.  (VI-12b)  gives 


At 

e- 


2  2  4  4 

1-a  t  /2'.4a  t  /4 !+. .. 

at-a3 13/3 +a5t5/5 

cos  jat  -j  sin  jat 

j  sin  jat  cos  jat 


3  3  5  5 

-at+a  t  /3!-a  t  /5 ! +. . . 

2  2  4  4 

1-a  t  /2 !  +a  t  /4!  , 

(VI-13c) 


In  the  case  of  the  zero  input  situation  with  k-vector  jz  we  simply 

multiply  Z(t  )  by  z(t  )  to  get 
<w  o  «*  o 

z(t)  =  <J>(t,t  )z(t  )  (VI-14) 

—  —  O  o 

which  yields  the  zero  input  response.  If  f  /  0,  then  by  treating  z 
as  the  output  we  can  apply  the  fundamental  decomposition  of  Eq.  (1-4). 

In  the  time-invariant  case  we  then  wish  to  convolute  the  impulse  response 
l(t,0)l(t)  with  £(t),  where  l(t)  is  the  unit-step  function.  Thus 
the  general  solution  of  interest  to  Eq.  (V-lla)  is 

t 

z ( t)  =  e-(t_to)z(t  )  +  [e£(t"T>f(T)dT,  t>  t  (VI-15) 

t 

o 

One  can  check  that  this  latter  is  a  solution  by  direct  substitution  in 
the  original  differential  equation. 
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Several  points  of  observation  are  worth  observing.  We  see  that  In 
the  time-invariant  case  the  transition  matrix  is  found  by  summing  an 
infinite  series.  Since  the  series  is  always  uniformly  convergent  one 
can  use  the  series  summation  as  a  method  for  finding  the  transition 
matrix  on  a  digital  computer.  Such  a  method  involves  only  summation 
and  matrix  multiplication  and  the  error  after  a  finite  number  of  terms 
are  considered  is  relatively  easily  determined  [3],  Alternate  methods 

result  from  noting  that  exp  At  is  the  inverse  Laplace  transform  of 

-1  ** 

(pi  -A)  as  Eq.  (Vl-llb)  shows.  Consequently,  all  entries  in  exp  At 

nk  **  ■* 

are  exponentials  or  time  multiplied  exponentials;  these  can  be  determined 
from  a  partial  fraction  expansion  of  (pl^-A)  1  where  in  fact  iterative 
methods  can  be  used  to  replace  evaluation  of  his  inverse  by  simple 
matrix  multiplications  [4]  [5],  If  also  f  has  a  rational  Laplace 
transform  the  final  z(t)  for  Eq.  (VI-15)  can  be  relatively  simply 
found  by  inversion  of  Laplace  transforms.  Alternatively  the  needed 
convolution  can  be  carried  out  directly,  though  less  conveniently,  on 
the  computer. 

F.  Discussion 

In  terms  of  the  realization  matrices  several  types  of  sensitivity 
have  been  discussed  and  evaluated,  all  for  scalar  transfer  functions  of 
time-invariant  networks.  Both  transfer  function  and  pole  position 
sensitivity  are  relatively  easily  evaluated  while  time  domain  variations 
require  a  solution  of  the  canonical  equations  to  find  the  transition 
matrix  exp  At. 

Actually  to  determine  the  variations  in  the  output  y(t)  due  to  x 
parameter  changes,  Sy/'tx,  requires  two  solutions  of  the  equations 
z=Az+f,  first  with  f  =  Bu,  with  z(t)=s(t),  and  then  with 

pft  ft*'  ft*  ft-  ftn()**0 

f  =  (? A/r*x)s  +  (rWrbc)u  subject  to  ,^s(t  )/c*x  =  z(t  ),  this  latter 
*1  «.—  ««'  *■>  o  —>o 

often  being  taken  as  zero.  Typical  results  in  the  somewhat  unrealistic 
situations  where  x  =  a  are  plotted  in  [2,  p.  341], 

Because  changes  in  responses  due  to  circuit  element  variations  can 
be  disturbing  it  is  often  desirable  to  try  to  find  circuitry  which 
minimizes  such  variations.  One  can  see  from  the  formula  T(p)  -  D  + 
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C(pl. -A)  B  that  if  the  entire  transfer  function  is  obtained  by  a  single 
realization  then  the  feedback  supplied  by  the  configuration  will  generally 
mean  that  each  circuit  element  can  possibly  stn  ugly  interact  with  all 
other  components  resulting  in  relatively  high  sensitivity.  On  the  other 
hand  if  the  transfer  function  is  broken  into  degree  one  or  two  factors 
c  b 

as  T(p)  =  :[d.  t - -  ]itr  D  +  C.(pl  -A.)  B  1  then  those  circuit  elements 

i  i  P+»A  ^  .1  r.j  -2-j  -l 

occurring  in  a  given  portion  only  relatively  strongly  Interact  with  those 
components  associated  with  t h<  appropriate  degree  one  or  two  realization. 
Consequently  there  is  practical  value  in  designs  based  upon  the  factor¬ 
ization  of  transfer  functions  into  small  degree  sub-portions. 

Finally  we  mention  that,  as  with  most  oihoi  state- variable  techniques, 
the  theory  of  sensitivity  is  made  practical  tor  the  use  of  digital  com¬ 
puters  through  the  techniques  discussed. 
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H. 

I. 

2. 


Exercises 


Exhibit  a  formula  for  y/'x  in  terms  of  the  realization  matrices 
and  the  initial  state  and  input. 

Show  [4]  that 


k-1 


(pl.-A) 

m  K  * 


-1 


_  Yu 


k-i-l 
dfn'  —  1 


t 
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where 


k  J  k-1 


d(p)  =  p  +  +  ...  +  d^p 


+  d. 


and 


B  =  1.  ,  B,  =  B  A  +  d,l.  ... 

»o  »k  -1  *0-1  l»k 


Ik-1  =  Ik-2*  +  “k-lik'  I  -  Ik-1*  +  \h 


3.  Find  the  sensitivity  of  the  Brune  section,  Fig.  1-5,  to  variations 
in  the  two  gyrators.  From  this  determine  which  gyrator  should  be 
most  stably  constructed. 

4.  Discuss  the  actual  programming  involved  in  setting  up  Eq.  (Vl-15) . 
Give  a  flow  chart  for  a  program  to  determine  ^iy/Sx  on  a  digital 
computer. 
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II  se  peut  qu'il  n'y  lit  aucune  arriere  - 
pens4e  entre  deux  hommes.  mais  il  y 
a  des  choses  plus  imp4rieuses  et 
plus  profondes  que  la  pens4e.  J'ai 
4t4  plusieurs  fois  t4raoin  de  ces 
choses,  et  un  Jour  je  les  ai  vues  de 
si  pres  que  Je  ne  savais  plus  s'il 
s'agissait  d'un  autre  ou  de  moi-meme  ... 

M.  Maeterlinck 

"Les  Avertis"  du  "Tresor  des  Humbles" 


CHAPTER  VII 

POSITIVE-REAL  ADMITTANCE  SYNTHESIS 


A.  Summary 

The  results  of  the  Positive-Real  Lemma,  whose  proof  is  merely  outlined, 
are  applied  to  obtain  a  transformation  which  yields  a  positive-real  coupl¬ 
ing  admittance  to  load  in  capacitors  such  that  a  passive  circuit  synthesises 
a  positive-real  admittance  Y(p) . 

B.  Introductory  Remarks 

Previously,  Section  III  D) ,  we  saw  that  if  an  admittance  matrix 
Y(p)  has  a  state- variable  realization  R  =  [A,  B,  C,  D)  then  a  physical 
structure  yielding  Y(p)  as  the  input  n-port  admittance  results  from 
loading  a  resistive  coupling  (n+k)-port  structure  described  by  the 
admittance 


Y 

*“>C 


(III-7  :) 


in  k  unit  capacitors.  Here  Y(p)  is  an  nyn  matrix  while  k  is 
the  size  of  the  state;  conveniently  k  is  taken  as  the  minimal  value 
b,  this  being  the  degree  of  Y(p).  The  structure  is  as  in  Fig.  VII-1 
which  is  Fig.  III-7  repeated  for  convenience. 
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Realization  Structure 
6  =  B[Y(p)]  for  Minimal 

Figure  VII-1 

However,  even  when  Y(p)  can  be  obtained  through  the  use  of  only  passive 
circuit  elements,  this  method  may  require  other  than  passive  elements 

since  Y  may  not  be  obtainable  without  the  use  of  active  elements. 

«-c 

Consequently  we  recall  that  all  minimal  equivalents  can  be  obtained  by 
transformations  performed  upon  Y^,  thus  all  minimal  capacitor  structures 
result  by  allowing  T  to  vary  in 

Y  =  '  B  -Cl'  1  (VII-1) 

»*• 

t”ab  -t-1at 

Our  interest  here  is  to  search  for  a  proper  choice  of  the  transformation 
T  such  that  the  new  coupling  admittance  matrix  can  be  realized  by 

passive  resistors  (and  gyrators,  recall  Fig.  III-8), 

We  recall  that  the  condition  for  a  given  rational  nyn  matrix  Tf(p) 
to  be  the  admittance  matrix  of  a  passive  n-port  constructed  of  only 
passive  circuit  elements  is  that  Y(p)  is  positive-real  [1,  p.  240]. 

By  definition  a  matrix  Y(p)  is  positive-real  if 

a)  Y(p)  is  holomorphic  in  Re  p  >  C 

b)  Y(p*)  =  Y*(p)  in  Re  p  >  0 
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c)  The  Hermitian  part  Y„(p),  2T„(p)  =  Y(p)  +  Y*(p),  is 

nonnegative  definite  in  Re  p  >  0, 

where  the  superscript  asterisk  denotes  complex  conjugation.  If  _Y(p) 
is  positive-real  and  rational  we  will  call  it  PR  for  convenience. 

Since  it  is  known  that  any  rational  positive-real  matrix  has  a 
passive  synthesis  in  the  form  of  Fig.  VII-1,  it  is  then  a  matter  of 
searching  for  a  suitable  transformation  T  to  make  $c  positive-real 
when  Y(p)  is.  The  purpose  of  the  next  sections  is  to  obtain  the  desired 
T. 

m 

C.  The  PR  Lemma 

First  we  recall  that  any  PR  matrix  Y^p)  can  be  decomposed  into 
the  sum  of  two  matrices 


Y (p)  =  Y.(p)  +  Y  (p)  (VXI-3) 

where  Y  ,  the  lossless  part  has  all  its  poles  on  the  ju)  axis  [and 

satisfies  Y  (p)  =  -Y  (-p) ]  while  Y  (p)  has  poles  only  in  the  open 

***  L 

left  half  plane;  both  Y.  and  Y  are  PR  while  the  decr.-roosition  can 

~L  -  o 

be  obtained  through  a  partial  fraction  expansion.  Since  the  poles  of 

Y.  and  Y  can  net  coincide,  a  minimal  realization  for  Y.  can  be 
-L  -  o  —L 

"added"  to  a  minimal  realization  for  Y  to  obtain  one  for  Y.  As  a 

■w  o  — 

consequence  we  will  first  obtain  properties  of  these  separate  realizations 
and  then  show  how  to  combine  them  to  give  the  proper  meaning  to  the  word 
"added."  For  convenience  we  assume  Y(°°)  =  Y  (») . 

m  «0 

The  basic  result  in  the  theory  is  as  follows  [2]. 

The  PR  Lemma:  Let  Y(p)  be  an  nxn  rational  matrix  with  real  co¬ 
efficients  and  with  no  poles  in  Re  p  >  0,  and  let  R  =  (A,  B,  C,  D} 
be  a  minimal  realization.  Then  Y(p)  is  PR  if  and  only  if  there 
exist  matrices  W^,  L,  and  a  (unique)  positive  definite  (symmetric) 
P  satisfying 


PA  +  AP  =  -LL 

mw*  ■■  — 


(VII-3a) 


PB  =  C  -  LW 

W  M  Mm  00 


(VI 1-3 b) 
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f 


w  w 

H  •  '*W 


=  D  4  I> 


(VII-ic) 


Outline  ni  Demons t vn  1 1  on :  is  the  steps  in  tie  proof  are  eh  tailed  and 
Involved  f2’  we  merely  outline  i he  main  ideas  with  emphasis  upon  those 
points  of  interest  for  actual  calculations. 

To  see  that  if  Eqs.  (VII-3)  hold  then  Y(p)  is  Pit  is  strnifiht- 
ioiw.ml  from  the  followiim  cal enl at  ions,  since  Y(p)  i  w  .<■  mim  <1  holo- 
morphic  in  Re  p  .•  h  anil  has  real  coefficients. 

2Y„(pl  =  Y ( p  ’ )  *  Y(p)  =  I HD  4  li(p‘]f-A)'1(i  C(pJ;  '  R 
-  W  W  4  U|  (p*l.  -A)_1P  +  P(plt  -A)_1  10 
4  B(p*l.-A)'1LWi  4  WJL(plf-A)  1 B 
=  WJV^  +  B(p*lf  -A)_1[P(p+pv)4-LL  )(pj;  -A)  ’j) 

4  B(p*l.  —A ) ~ 1 1 AV  +  W  L(pl.-A)"1H 


I  !L+e(e %  -£> II  h  "i)  I 


4  B(p *1,  -A)  *[  (p+p*)J>  |(pl; -A)  1  R 


(VI 1-1) 


This  last  shows  that  Y  (p)  is  positive  semidefinite  for  all  p  with 
pip  0,  that  is  in  tin.*  ri^ht  half  plane,  since  P  ean  also  be 
tactored  into  P  -  P1  ^P1^2  with  the  square  loots  also  symmetric. 

To  show  that  Y(p)  is  PR  only  if  Eqs.  (VII-.'!)  hold  is  more 
difficult.  We  first  iincl  a  W(p)  satisfying 


Y (p)  +  Y(-p)  -  W (-p)W(p) 


(VII -5a) 


will  Me  tnvt.lii-i  W(p)  is  hoi  omorpli  ie .  top;.  Uier  with  its  ripht  inverse 

•.V  ,  Mi  111'  riitlit  halt  plain  .  Sueh  \V(o)  inn  In'  lound  eon  wn  h  n  I  I  v . 

-  r 

Im;  l  lie  calculations  can  become  involvid  |'  j,  The  use  t  ii  i  ■  nirti  ail.n 
W(p)  is  used  'o  "iia  >  a  alee  ine  in  i  n  l  "l  1 1  l  t  y  o|  U  <  j  > )  .  eliia  i:ip|ei 
I  ac  tori  za  I  ions  as  the  mu  ol  r.aie...  I  i  n  LuSl  .au.H1  J:1,  i<-.  ."I  ■.  an  t . . 

i  NOT  w  v.s .-wjii. 

<)7  .  !•  I  —  ( >  V  —  1  in 


One  can  then  show  that  W(p)  has  the  minimal  realization  R  a 
[A,  3,  L,  W  j  which  serves  tc  define  L;  note  that  the  matrices  A 

m  m  7  m co  »m. 

and  B  are  identical  for  Y(p)  and  W(p).  We  then  transform  the 
minimal  realization 


R1  = 


'  A 

o 

j 

‘  B 

.LL 

-A  . 

LWm 

0*mO$ 

»  [WJL,  -B],  WJT,  I  (VII-5b) 


of  W(-p)W(p)  through  Eq.  (V-2a)  using 


T  = 


h  2 
'-p  i*  J 


to  get  the  equivalent  realization  for  W(— p) W Cp> 

3  1,  [W  L+BP,  -B],  W  W 


R2  = 


A  0 
0  -A 


PB  +  LW  J 

+  m  m*  mca 


(VII-5C) 


(VII-5d) 


Here  P  is  the  unique  positive  definite  solution  of  the  equation 

PA  +  A?  -  LL  (VTJ-3|a> 


Next  we  note  that  a  realization  for  Y(p)+Y(-p)  is 


D  +  D 


(VII-5e) 


On  noting  the  conditions  for  equivalence  and  identification  of  realizations 
we  obtain  R2  = 

On  noting  that  almost  all  of  the  previous  holds  except  that  W  =  0,- 
and  hence  L  =  0,  when  Y  is  lossless  and  zero  at  Infinity,  we 

la  m  m 

conclude  that  in  the  lossless  case  there  exists  a  positive  definite 
(symmetric)  JP  such  that 


R,  and  the  PR  Lemma  follows,  Q.E.D, 
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0 


(VII-6a) 


£iAl  + 

*a  =  k 


(VII-6b) 


where  (A^,  C^,  0}  is  a  minimal  realization  of  the  lossless  PR 

admittance  which  is  zero  at  infinity.  As  a  consequence  we  can  replace 
the  conditions  of  the  PR  Lemma  to  allow  simple  poles  on  the  jw  axis, 
none  though  at  infinity ,  if  we  use 


P  =  PT  +  P  ,  A  =  A  +  A 

•>  *»L  •'O  —  **L  <~o 


B  = 

V 

y 

C  =  [CT,  C  ], 
m  L*L7  ■*0J, 

II 

'  0 

■?o- 

L 

"•o 

where  the  subscript  zeros  refer  to  the  realization  of  Yo,  that  portion 
of  Y(p)  with  only  open  left  half  plane  poles.  Note,  however,  that  now 
P  is  no  longer  unique  by  virtue  of  the  presence  of  P  . 

^  L 

In  conclusion,  if  Y(p)  Is  PR  with  no  pole  at  infinity  then 

Eqs.  (VII-3)  hold  with  the  various  matrices  obtained  using  Eq.  (VII-7) 
upon  decomposing  Y(p)  into  the  sum  of  a  lossless  part  Y^tp)  and  a 
nonlossless  part  YQ(p).  The  calculations  are  theoretically  very 
straightforward  but  the  computation  for  W(p)  with  the  proper  holomorphic 
inverse  gives  considerable  difficulty  in  practice.  However  once  such  a 
W(p)  is  found  Eqs.  (VII-3a)  can  be  solved  for  Pq  in  a  very  straight¬ 
forward  manner  as  a  set  of  linear  algebraic  equations  subject  to  the 
positive  definite  constraints.  As  it  stands  the  method  does  not  allow 
the  direct  treatment  of  poles  at  infinity  and  these  must  therefore  be 
extracted  separately  as  an  added  term  pC^,  for  the  right  ^ide  of  Eq. 
(VII-2),  to  be  independently  considered  for  synthesis  purposes. 


D.  PR  Admittance  Synthesis 

We  assume  as  given  an  nxn  PR  admittance  matrix  which  we  can,  as  a 
consequence,  decompose  into 
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Y(p)  =  Y  Cp)  +  Y.  (p)  +  pC 

«o  *«L 


(VII-8) 


where  Y  (p)  is  holomorphic  in  Re  p  >  0,  all  the  poles  of  Y  (p)  are 
on  the  Joj  axis  and  simple  with  none  at  infinity,  and  all  three  terms  on 
the  right  of  Eq.  (VI-8)  are  separately  PR.  The  term  pC^  is  separately 
synthesized,  using  for  example  only  capacitors  loading  transformers 
[1,  p.  204];  the  resulting  network  for  pc  is  connected  in  parallel 


with  one  of 


Y„  +  Yt  . 

»o  — L 


To  synthesize  Y  +  Y,  we  find  any  minimal  realization 

**o  — L 

R  =  {[A^+Aq],  I" Bl,  [Cl,Cq],  and  t*len  determine  a  desired  JP  = 


’h' 

B 

*-  ~oJ 


P. +P  as  for  Eq.  (VII-7).  Since  P  is  positive  definite  we  find  its 

*L  "°  -  1/2 

(unique)  positive  definite  square  root  P  .  Thus 


P  =  P1/2P1/2 


(VII-9a) 


In  actual  fact,  since  P  is  in  direct  sum  form  we  can  also  write  P 
in  direct  sum  form  as 


=  pt1/2  ;  p1/2 

—  L  — o 


(VII-9b) 


Next  we  apply  the  theory  of  equivalence  of  Chapter  V,  choosing 


T  =  P 


(VII-9c) 


-1/2  1/2 

where  P  is  the  inverse  of  P  [note  that  the  P  of  Eq.  (V-2b) 

has  a  different  meaning  than  the  P  of  Eq.  (VII-9c)  whereas  the  T's 

are  the  same].  We  then  have  a  realization  ft  =  {P1/,2AP-1,/2,  P1//2B, 

-1/2 

CP  ,  D}  derived  from  the  original  R  having  its  entries  as  given  by 
Eqs.  (VII-7).  As  a  consequence,  by  our  introductory  comments  and  Eq. 
(VII-1)  we  can  form 


p1/2b  -p1/2ap-1/2 

MM  Ml  W«M 


(Vll-lOa) 
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-i(  Vll-lOb) 


Y  (»)  +  Y  (^) 

•»L  »o 

L*nJ 


(pt1/2;p1/2 


~i  c  c  )(P-1/2 ;  P-1/2, 

1  *L“o  -L  — o 

-(pJ/2+P1/2)[A  +A  ](P"1/2;P"  /2) 
~L  -o  ~L  ~o  -L  -o 


By  our  previous  reasoning  Y(p)-pC  results  from  loading  the  resistive 

m  ^03 

coupling  network  having  the  admittance  matrix  Y  in  6  unit  capacitors, 

—  c 

where  6  is  the  degree  of  Y(p)-pC  .  Our  claim  is  now  that  1?  is  PR 

if  Y(p)  is,  such  that  a  circuit  structure  from  Y ^  need  use  only 

*  "“C  -1/2 
5jriatcrs  and  positive  resistors.  That  is,  the  choice  T  =  P  has 

allowed  a  completely  passive  synthesis  of  a  PR  admittance  matrix. 

To  see  that  $  is  PR  we  merely  need  to  cheek  to  see  if  it  has 
~c 

a  positive  semideflnlte  Hermitian  part.  Thus  we  form 


Y  +  9  = 

»*c  -  c 


D  +  D 
%<  -» 

p1/2b  -  r1/2c 


-CP"1/2  +  BP1/2 

-p1/2ap'1/2  -  p"1/2ap1/2 


(Vll-lla) 


(1n  ;  P~1/2) 


W  W  -C  +  BP 

k  00  >00  —  ~  «- 

L  PB  -  C  -PA  -  AP  J 


(1  +  P"i/2>  (Vll-llb) 

**n 


(1  ;p:1/2;p-1/2> 

•.n  — 'L  —o 


W  W  0 

M  00—00  » 


0  0 


-L  W  0 


-W  L  -I  (1  +p"1/2;p_1/2  (VII-llc) 
"•=o*o  |  »n  *iL  ~o 

0 


L  L 
*o-o  -i 


where  we  have  used  the  fact  that  P1//2  is  symmetric,  p1//2  =  P1//2,  as 
well  as  Eqs.  (VII-3)  in  their  extended  form  valid  for  the  inclusion  of 
lossless  parts,  Eq,  (VII-7)  .  That  is,  is  that  W(°°i  which  corre¬ 
sponds  to  Y  (p)  while  L  =  [0,  L  1 ,  If  W  has  rank  r,  that  is  if 

r  is  the  rank  of  Y  (p)  +  Y  (p) ,  then  we  can  rewrite  Eq.  (VII-llc)  as 

m  O  O 


mC 


»«“  C 


i  [w  ,  o,  -i*  p"1/21 

—  r  L-cc7  '  —  o*"o 


(Vll-lld) 
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As  shown  by  Section  III  D),  Yc  can  now  be  synthesized  by  gyrators  and 
r  positive  resistors.  For  instance  Fig.  III-8  applies  to  synthesize 
the  symmetric  part,  which  is  one-half  of  Eq.  (Vll-lld),  with  r  unit 
resistors  and  a  gyrator  coupling  network  described  by  the  gyrator  conduc¬ 
tance  matrix 


G 


0 


(Vll-lle) 


We  comment  that  zeros  in  G  which  designate  rows  and  columns  of 

zeros  in  the  symmetric  part  of  9^  are  as  expected  since  they  are 

associated  with  the  lossless  part  YT (p)  for  which  no  resistors  are 

*>L 

necessary.  In  fact  since 


r  =  rank[Y(p)+Y(-p) J  (Vll-llf) 

and  since  this  rank  corresponds  to  the  minimum  number  of  resistors  pos¬ 
sible  in  a  synthesis,  we  see  that  besides  using  a  minimum  number  of 
capacitors  this  method  uses  the  minimum  number  of  resistors.  In  fact  in 
the  case  where  the  original  Y(p)  is  lossless,  j3  of  Eq.  (Vll-lle) 
reduces  completely  to  zero.  Of  course  the  vanishing  of  the  P 
portions  of  Eq.  (Vll-lld)  does  not  mean  that  P,  never  enters  into 
consideration;  for  example  P  occurs  in  the  skew- symmetric  portion 
which  acts  through  gyrators  to  couple  the  capacitors  to  the  input  ports 
in  a  lossless  manner. 


E.  Example 

Let  us  apply  the  method  to  the  PR  scalar 


y(p) 


4p3  +  2p2  +  18p 
3  2 

p  t  2p  +  4p  +  8 


(VII-12a) 
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f 


•Ip 

p+2 


( VII-12b) 


2p 


p2  *  . 


The  latter  split  Rives  the  decomposition  into  lossless  and  nonlossless 

parts;  thus  y  (p)  =  2p/(p2+4),  y  (p)  =  4p/(p+2). 

L  o 

For  v.  and  y  appropriate  realizations  R  and  R  are  obtain- 
'  I,  '  o  L  o 

ed  from  Eq.  (IV-3)  as 


\ 

0  ll, 

( 

2 

) 

.  -4  0j 

.  n . 

!  1,  0  i ,  i 0)  j 
( 


(VII-12c) 


R  =  {[-2  !,  i  -81.  f  1  !.  !  1  n 


( VII-12d) 


For  y  we  have 
o 


y(p)  +  y  (-p)  = 


-Bp 


=  ( vii-i2e) 

-pj-2  p+2 


-p  +  '1 

We  observe  that  w(p)  is  unique  to  within  a  minus  siRn;  we  choose 


W  (p)  = 


2  s!  2p 
p+2 


(VII-12f ) 


for  which  a  realization  following  Eq.  (IV-3)  is  ([-2],  (-4/2],  [1], 
[2/2]}.  We  thus  desire  to  choose  a  transformation  T  -  1/  /  2  to  bring 
this  B  of  -4/2  to  T  =  -8.  Thus  we  have  as  the  appropriate 

realization  R  for  W 

w 

Rw  =  ([-2],  [-8],  [  1//  2.1 ,  |  2 /2]1  (VII-12g) 

We  have  at  this  point  Lq  =  \/J~2  and  V.'^  =  2/2.  T!ie  transformation 

P  is  found  from 

o 


P  A  +  A  P  =  -4P  =  -L  I,  =  -1/2  (VII-12h) 

00  00  o  00 


or 


P  =1/8,  P7  =  1/2/2 

0  0  v 


(VII-12i) 
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To  find  T  we  observe  that  y  (p)  +  y, (-p)  =  0  in  which  cace 
L  L  Li 

1,^  =  0  and  we  simply  solve  for  a  positive  definite  satisfying 

-PL*L  +  «^L-L  =  —2’  that  18 


pli 

P12 

'  0 

1 

+ 

0 

-4  ' 

P11 

P12 

— 

’  0 

0  ■ 

-  P12 

P22- 

.-4 

0 

1 

0 

P12 

P22  . 

0 

0 

(VII-12J) 


The  (1,1)  and  (2,2)  entries  require  p  =  0  while  the  (1,2) 

L  m 

entry  gives  p^  =  4 p22  as  does  the  (2,1)  term.  Positive  definiteness 
merely  requires  p  >  0  while  P  B  =  C  requires  8p  =  1.  Thus 

•**  1  ^  1 1  **  1 1  /  Z 


II 

1/2  0 

P1/2  - 

’  -L 

‘  i//*2  o 

-  0  1/8. 

.  0  1/2/2*. 

(VII-12k) 


.1/2 


Now  the  original  coupling  admittance, 
is 


before  the  application  of 


4 

-1  0 

-1  ‘ 

( VII-121) 

2 

0  -1 

0 

0 

4  0 

0 

-8 

0  0 

2  - 

which  is  not  PR 
We  then  form 


as  can  be  seen  by  the  principal  middle  submatrix 


-V 


[i  ;  p1/2]y  [i  ;  p'1/2j 

w»  Jp~  C  v-  J 


where 


1/2 


_l/2  •  „ 
=  P  +  p 
"L  o 


or 


(VII-12m) 
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$ 

>-c 


-J2  0 


0-2  0 
2  0  0 

0  0  2 


-1 

0 

4 

0 


0 

-1 

0 

0 


0 

0 

0 


0 

0 

0 


Finally  we  have  for  ?  =  ?  +  $ 

*c  -c  sy  —c  sk 


Note  that  $  takes  the  form  predicted  by  Eq.  (Vll-lld).  The  final 

^  c  sy 

circuit  diagram  is  shown  in  Fig.  VTI-2. 

In  the  figure  we  observe  that  y  and  y  arc  separately  realized 

Li  O 

and  then  connected  in  parallel.  In  all  situations  Y  will  be 

wc  sy 

associated  only  with  yQ(p)  but  in  this  case  the  skew  symmetric  part 

has  only  occurred  while  synthesizing  y  (p) .  Note  that  even  though  a 

L 

minimum  number  of  capacitors  and  resistors  have  been  used  an  excess 
number  of  gyrators  occurs.  By  shifting  elements  through  the  gy-ators 
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we  can  easily  obtain  Fig.  VII-3a)  from  Fig.  VII-2,  or  by  direct  synthesis. 
Decomposing  this  latter  circuit  yields  '.he  resistive  circuit  of  Fig. 
VII-3b),  loaded  in  capacitors.  We  observe  however  that  this  latter 
configuration  possesses  no  admittance  coupling  matrix,  on  •  conclusion 
is  that  always  our  synthesis  of  PR  admittances  will  worn  but  that  in 
some  instances  more  than  the  minimum  number  of  gviators  will  be  UBed, 
though  never  more  than  the  minimum  number  of  capacitors  and  resistors 
is  needed. 


1/4 

a)  b) 

Minimal  Gyrator  Realization 


Figure  VII-3 
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F.  Discussion 


In  this  chapter  we  have  presented  a  method  oi  synthesis  of  positive- 
real  rational  admittance  matrices,  and  by  duality  impedance  matrices. 

The  method  is  based  upon  only  algebraic  operations  and  thus  is  readily 
programmed  on  a  computer.  The  key  point  of  Use  theory  is  the  proper 
application  of  the  PR  Lemma  to  obtain  the  appropriate  transformation. 
However  it  is  in  the  application  of  this  Lemma  whoie  tin  greatest 
difficulty  occurs  since  a  rather  complicated  factorization  of  the  para- 
Hermitian  navi  of  Y f p '  sometimes  must  be  undertaken  in  order  to  Obtain 

W(p) .  For  nonposi t ivo-real  matrices  or  positive-real  matrices  of  infinite 
dimension  similar  steps  appear  to  be  possible  but  as  ye:  have  not  been 
extensively  studied. 

The  ideas  of  the  method  can  bo  applied  '<  „  hybrid  coupling  matrix 

in  such  a  manner’  that  some  promise  holds  tor  minimal  gyvator  synthesis 
r4l.  That  is.  Y  can  be  interpreted  as  n  hybrid  matrix  if  some  ports 

L  J  -<  C 

are  loaded  in  inductors  in  place  of  capacitors;  in  such  a  case  one  still 

desires  Y  PH  when  Y(p)  is.  Alternately,  by  using  ihe  hybrid 

•»  c  — 

interpretation  one  can  give  a  synthesis  in  terms  of  the  cut  set  and  tie 
set  matrices  previously  studied,  at  least  in  the  lossless  (and  gyrator- 
less)  case  [5].  However,  as  with  the  minimal  gyrator  situation  improved 
conditions  are  still  needed  to  complete  the  theory.  Nevertheless  the 
nonlinear  theory  has  been  interestingly  investigated  1 6  . 

Because  of  the  situation  illustrated  in  Fig.  Vll-d,  where  no  coupl¬ 
ing  admittance  matrix  exists,  it  seems  important  to  extend  the  method  to 
scattering  matrices  where  partial  results  of  the  Pit  Lemma  type  are 

available  [7],  The  work  of  Youla  and  Tissi  represents  a  step  in  this 

direction  [8]. 

Since  it  was  possible  to  find  one  trans tor-mat  ion  T  taking  any 

minimal  realization  into  a  passive  one  it  is  of  intei  ■  t  to  find  all 

such  T.  As  yet  little  solid  theory  is  available  jn  this  direction  but 
the  theory  of  continuous  transformation  groups  seems  applicable. 
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H.  Exercises 

1.  Synthesize  the  PR  impedance  matrix 


1 

P 

1 

ri  "2| 

2  „ 
p  +2 

1  4p . 

p+1 

-  2  4  J 

bv  converting  to  the  admittance  matrix  ami  applying  the  methods  of 
this  chapter. 


108 


S EL- 6 7- 11C 


2.  Fill  in  the  steps  of  the  PR  Lemma  proof. 

*3.  Discuss  a  positive-real  lemma  for  nonrational  matrices  and  how 
this  night  be  used  for  synthesis. 

4.  Investigate  possible  methods  of  determining  W (p> ,  [3],  [1,  p.  168  j 
and  discuss  the  simplest  for  machine  calculation. 

5.  Show  that  the  synthesis  of  the  text  uses  both  the  minimum  number 
of  resistors  and  capacitors. 

*6.  For  the  example  of  the  text: 

a)  find  all  minimal  realizations  and  isolate  those  for  which 

Y  is  PR. 

•*c 

b)  investigate  possible  ways  of  accounting  for  Fig.  VII-3. 

c)  find  all  minimal  realizations  on  a  scattering  matrix  basis. 

7.  Discuss  the  various  methods  of  calculating  the  matrix  P  [9], 

*8.  Investigate  methods  of  synthesizing  bounded-real  rational  matrices 
by  the  techniques  of  the  text  [7]. 

*9.  Show  how  the  same  techniques  can  be  extended  to  cover  nonminimal 

synthesis  of  PR  matrices  and  discuss  how  such  may  be  of  importance 
for  minimal  gyrator  synthesis. 

10.  Apply  the  PR  Lemma  to  show  how  to  synthesize  through  the  equations 

[10] 


s  =  •i(A-A)  s  +  Bv - Lv* 

~  2  *»  w  "  —  J2  — 

i  =  Bs 

mm.  W- 


subject  to  i  =  -v*.  In  particular  show  that  a  network  realization 
occurs  by  terminating  the  gyrator  network 


Y 

**  c 


0 


0 


-5 


0 


v/2 


J_  L 
i(A-A) 


in  unit  resistors  and  unit  capacitors.  Show  that  the  minimum  number 
of  resistors  and  capacitors  are  used. 
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II s  semblaient  par  moments  nous 
regarder  du  haut  d'une  tour.  II  est 
vrai  que  rien  n'est  cach£;  et  vous 
tous  qui  me  rencnntrez,  vous  savez 
ce  que  j’ai  fait  et  ce  que  je  ferai, 
vous  savez  ce  que  je  pense  et  ce 
que  j 'ai  pensd. 

M.  Maeterlinck 

"Les  Avertis"  du  "Tresor  les  Humbles" 

CHAPTER  VIII 

LUMPED-DISTRIBUTED  LOSSLESS  SYNTHESIS 

A.  Summary 

Here  we  briefly  summarize  the  application  of  the  previously  discussed 
techniques  to  the  synthesis  of  networks  constructed  of  lossless  lumped 
circuit  elements  and  LC  transmission  lines.  The  theory  is  based  upon 
the  use  of  frequency  transformations  to  obtain  lossless  but  r'  onal  2- 
variable  matrices. 

B.  Introductory  Material 

We  first  review  some  properties  of  LC  transmission  lines  as  well  as 
a  method  of  treating  circuits  constructed  from  lumped  circuit  elements 
in  conjunction  with  the  LC  lines.  This  will  lead  us  to  positive-real 
and  rational  2-variable  matrices  and  their  synthesis.  As  we  will  see 
the  admittance  description,  which  we  adhere  to,  is  not  rational  in  the 
true  frequency  variable,  and  as  a  consequence  we  introduce  a  second 
frequency  variable  to  obtain  rationality. 

Let  us  first  consider  a  lossless  transmission  line  of  length  l 
and  inductance  L  and  capacitance  C  per  unit  length.  As  shown  in 
Fig.  VIII-1  this  line  can  be  treated  as  a  2-port  having  the  admittance 
matrix  [1,  p.  66] 

T  ctnh  yp  -csch  yp  1  y  =  J'hC  l 
Y (p)  =  y  (VI I 1-1 a) 

tm  O 

.-csch  yp  c  anh  yp  J  y  ^  =  ,/c/L 
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i 


1 


Lossless  Transmission  Line 
Figure  VIII-1 


We  observe  that  the  admittance  matrix  Y(p)  is  not  rational  in  p  but 
that  *he  positive- real  frequency  transf ormation 

X  =  ctnh  (yp/2)  (VTlI-lb) 


yields  a  rational  positive- real  admittance  description 


r  i+x2 

i-x2 

2X 

2X 

1-X2 

1+X2 

-  2X 

2> 

(VIII-lc) 


In  fact  we  observe  that  any  transmission  line  which  has  its  y  an 
integer  multiple  of  this  basic  line  also  has  an  admittance  matrix  which 
is  rational  in  X.  Since  given  a  set  of  transmission  lines  for  which 
the  r's  are  rationally  related  there  always  exists  a  smallest  y  for 
which  the  admittance  description  is  Eq.  (VIII-lc),  we  will  assume  that 
all  lines  under  consideration  arc  rationally  related,  that  is  have 
rationally  related  y's. 

If  next  we  assume  the  presence  of  lumped  capacitors,  inductors  and 
gyrators,  as  well  as  the  rationally  related  LC  lines  considered  in  the 
X  piano,  a  node  analysis  yields  branch  admittances  of  the  form 

yij  =  CijP  +  (1/Lijp)  +  cijX  +  (1/eijX)  +  Bij  (VIII-2a) 

and  for  i  ^  j 

y  .  .  =  C.  .p  +  (l/L.  .p)  +  c.  .X  +  (1A\  .X)  -  g.  .  ( VI II-2b) 

)i  ij  i.l  i.l  i.l  i.l 
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Such  a  network  we  will  call  lumped  distributed.  We  note  that  for  passive 
elements  a  lumped  distributed  network  has  an  admittance  matrix  y(p,X) 
at  any  ports  which  is  positive-real  in  both  variables  and  satisfies  the 
lossless  constraint 


y(p,>)  =  -y(-p,-X)  (VIII-2C) 

I* 

In  actual  fact  y(p,X)  satisfies  the  2-variable  posi ti ve-real  constraints. 
Thr.t  is,  by  definition  a  matrix  is  2-variable  positive-real  if  \2,  p.  252  | 


a) 

y(p,X)  is  holomorphic 

in  Re  p  >  0,  Re  X  > 

0. 

b) 

y(p,\)  is  real  for  p 

and  X  real  in  Re 

p  >  0 

Re  X  >  0, 

c) 

the  Hermitian  part  of 

y(p,X)  is  positive 

r«v*  • 

semi- 

definite  in  Re  p  n  0, 

Re  X  0. 

A  rational  2-variable  positive-real  matrix  will  also  be  called  PR. 

A  property  of  interest  for  synthesis  is  that  the  poles  on  the 
imaginary  axes  can  be  separately  extracted  to  yield  [3,  p.  34] 

y(p,X)  =  y,  (p)  +  y„(X)  +y  (p,X)  (vin-3) 

where  y  ,  y  ,  and  y  are  all  separately  positi ve-real  and  rational 

*»o  m.1  i—2 

when  y  is  rational;  here  y  has  only  poles  which  explicitly  depend 

«■»  >*•  o 

upon  both  p  and  X.  Of  primary  interest  is  the  fact  that  y  (p,X) 

«***  o 

has  no  polos  at  infinity  in  either  variable. 

To  head  toward  synthesis  it  is  of  importance  to  note  that  those 
lines  which  have  lengths  one-hall  the  basic  length,  called  unit-elements, 
are  described  by 


Y  (p) 
-uo 


/x5 


/?..l 


(VII 1-4) 


Although  such  a  description  is  not  rational  we  observe  that  when  loaded 
in  a  short  circuit  the  unit-element  appeals  as  a  capacitor  of  capacitance 
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y  in  the  X-plane  when  observed  at  the  input.  Similarly  loading  in  an 
o 

open  circuit  yields  a  X-plane  inductor  at  the  input .  These  relation¬ 
ships  can  be  depicted  as  shown  in  Fig.  VII1-2. 


25 

uni  *  - 

ol  onion  l 

V 

o 

0 

1  C 

a) 


•  -plane 


I  <r-  £/2  — >  | 

•~Q 


*1 


c  =  yn 


|<~  1/ 2  -)| 

° — o  yQ  o-* 

p-plane 


b) 


X  -pi  ant 


p  vs.  X- Plane  Elements 
Figure  VIII-2 


With  this  last  observation  we  see  that  a  synthesis  method  could 
possibly  arise  by  loading  a  p-plane  (n+c)-port  described  by 


Y  <P) 


J)2(i”' 

•y21(p) 

y22<P)- 

(VIII-5a) 


by  a  set  of  c  unit  l^-plane  capacitors  (which  are  p-plane  shorted  unit- 
elements),  as  shown  in  Fig.  VIII-3.  If  such  occurs  then  one  obtains 


y(p,x>  = 


Jn!p) 


i2<P)i'ic'-yV2 


(l>) 


J21(>5) 


( VI II-5b) 
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y(p  A) 

W* 


1 

mC 


CX)  =  xi 

*-c 


Possible  Configuration 
Figure  VI I 1-3 

We  observe*  tlie  following.  In  the  general  expansion  of  a  2-variable 
PR  matrix,  Eq.  (VIII-3) ,  the  matrix  y  ( X )  can  not  be  absorbed  in  Eq. 

ma2 

(VIII-5b)  while  y  (p)  can.  However,  both  y  and  y  can  be  synthe- 
sized  by  standard  methods  with  the  resulting  networks  being  placed  in 
parallel  with  that  for  y  Hence  we  really  need  only  consider  Eq. 

(VIIl-5b)  for  y  (p,X).  Now  Eq .  (VIII-3b)  is  in  the  form  of  previous 

—  o 

results  except  that  the  realization  matrices  vary  with  p.  Thus  we  are 

after  a  realization  R(p)  =  (A(p)  ,  B(p)  ,  C(p)  ,  D(p)}  =  ^22’  Zl2’  ~l21’  Jll^ 
in  which  case  the  previous  theory  should  hold.  In  fact  we  can  use  the 
methods  of  Chapter  IV  to  create  a  minimal  realization  R(p) .  However,  the 

transformation  to  bring  Y  (p)  to  be  PR  though  obtainable  in  theory  is 

*~c 

not  known  in  explicite  form.  Thus  we  proceed  by  directly  finding  a  PR 
coupling  admittance,  this  being  possible  because  of  the  lossless  nature 
imposed. 

C.  Minimal  Realization  Creation 

To  obtain  a  realization  R(p)  for  an  n>n  PR  y  (p,X),  for  which 

O 

Y  (p)  is  also  PR  we  will  simply  modify  the  previous  realization  theory, 

^  c 

presenting  the  method  of  Itao  f4],  in  some  places  omitting  the  details  of 
proof  which  can  be  rather  lengthy  for  their  content. 

As  before  we  write 

'’iWp)  *  *  •••  *  ■'V”  1  Vi” 

- 7 - TVrri - V - 

X  +  a^ (p) X  +  .  .  .  t  a  (p) 


1 1*1 
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( VIII-6a) 


V  -i(p) 

=  A  ,  (p)  +  N  — 7TT- 


where  the  latter  is  the  expansion  about  > 
matrix  is  chained  as 


rnyrn  companion 


U  (p)  = 

»n 


-a  (p) 1  -a  (p)l 
1  -n  2  — n 


-a  ''o)l 

r  *»n 


(VIII-6b) 


and  a  modified  Hankel  matrix  defined  by 


T  (p)  = 

r 


VP> 

-Ai(p) 


;VP)  &<p) 


-A  (p> 


(VIII-6C) 


Because  of  the  lossless  nature  of  y  .  T  ( )  is  •  anal  to  T  f-p) 

wO  **«r  ) 

[the  para-Hermitian  property]  and  it  is  nonneica  1 1  vv  S'-m  i  de  f  i  n  i  t  e  for 

p  =  jio .  Consequently  T  (p)  can  tic*  fac  tored,  in  tact  by  the  method 

*-r 

used  at  Eq .  (VI 1 1- 5a),  to  obtain 


T  (p) 
*~v 


U(p)  l'(-p) 


(VIII-6d) 


where  U(p)  as  well  as  its  left  inverse  b  (p)  are  holomorpnic  in 

. -  . 

Re  p  >  0;  this  factorization  preserves  Ui<>  real -i  at  i  onal  nature  of  T  , 

—r 

that  is,  I'(p)  is  also  rational  with,  real  coo  I  f  i  c  ;  cuts .  Further  the 
matrix  l!  can  be  taken  of  size  nr>.f .  where  is  the  rank  of  T  (p) 

iw  \  } 

and  then  partitioned  into  n>b  blocks  to  define  tiie  e  tries  in 

A 


SEL-67-1 1 0 


(VIII-6e) 


u(P)  =  ryo(p),  u^p),  u^tP)] 

Noting  that  A  (p)  =  -y, 0(p)y0, (p)  =  U  (p)U  (p)(-p)  we  see  that  as  we 
*«*o  1  .-»o 

desire  y  (p)  =  -y.„(-p)  because  of  the  lossless  constraint,  we  are 

w>  21  1 2 

led  to  define 

y,9(P)  =  -y91(-P)  =  U  (p)  (VIII-7a) 

Noting  further  the  previous  method  of  defining  A  by  Eq.  (IV-9)  some- 
what  justifies  the  definition 

y._(p)  =  l/Acpltl  <P)U  (p)  ( VIII-7b) 

2.22  ^  i 

Of  course  we  also  define 


2ii(p>  -  Jo<p>"> 


(VIII-7C) 


With  these  the  coupling  admittance  matrix  of  Eq.  (VIII-5a)  is  completely 

specified.  In  fact  Y  (p )  is  PR  and  satisfies  the  lossless  condition 

*»  c 

Yc(p)  =  -Y  (-p)  though  both  these  properties,  especially  the  PR  one, 
are  rather  delicate  to  prove;  the  interested  reader  is  referred  to  [4], 

Further,  the  degree  of  Y  (p)  is  the  minimum  possible  and  equal  to  the 

wc 

p  degree  oi  y  (p,X)  defined  as  ■  =  max  ciy  (p,X  1.  The  number  of 

++  O  P  \  mO  o 

‘n 

X-plano  capacitors  is  equal  to  c  where  in  fact  c.  =  max  t;[y  (p  ,X!  = 

-  A  n  m^O  O 

1  o 

rank  T^(p) .  Wo  comment  that  the  whole  process  could  have  been  undertaken 
hy  making  p-plane  capacitor  extractions  from  which  we  conclude  that  b 

P 

represents  the  minimum  possible  number  of  p-plane  reactive  elements, 

while  b  gives  the  minimum  number  of  X-plane  reactive  elements. 

A 

In  summary,  loading  the  PR  (n+b  )>(n+b^)  matrix 

A  A 


Y  (pi  = 

w*»C 


y0(p’M)  Jk(p) 


-i 


L-io(-p)  Sj  <p>“n(-p)iio(p)J 


( VIII-7d) 
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realization  in  l  unit  ) -plane  capacitors  (which  are  shorted  unit- 
elements)  yields  y  (p,>)  =  y  (p,  ctnhryp/2l)  at  the  n  input  ports.  A 

^  O  *mO 

synthesis  of  the  lossless  coupling  admittance  Yc(p)  by  a  minimum  number 
of  reactive  p-plane  reactive  elements,  which  is  readily  possible  [1, 

Chap  8],  yields  a  network  possessing  a  minimum  number  of  lumped  reactive 
elements  as  well  as  (p-plane)  transmission  lines. 

D.  Exampl es 

Let  us  synthesize  the  function 


y(p,\) 

X(p/2)  f  1 

(VTTI-Ha) 

X  +  (f2+p2l/2p) 

We  have 

T  -1.1 L 

r  “  2  4 

(VlII~Hb) 

with 

V» 

if1  -  -2  -  -  - 

«  "  P\A’  " 

Z  <vm-*c) 

i  >r  which 

Y  (p)  - 

|  r  " 

-p+/2  P  ‘  -v 

(  VI 1  J-bci ) 

I'sing  a  short,  circuited  transmission  lint  for  the  -plane  capacitor 
yields  the  circuit  of  Fig.  VIII-4, 

1/2 


Example  Synthesis 
Figure  VIII-4 
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To  illustrate  the  difficulties  of  the  more  general  situations 


consider  the  lossless  PH 


y(p,x>  = 


2 

X  p+2X 

_ 

X  +pX+2 


(VIII-9a) 


The  expansion  about  X  =  “>  gives 


,  (2-p2)  (p3-4p)  (-p4+6p2-4) 

y(p,X)  =  p  +  — r -  + - 5 -  + - =r- -  + 


(VIII-9b) 


Thus 


I2 


2-p 


3  . 

P  -4p 


3  .  4  2 

-p  +4p  -p  +€p  -4  J 


(VIII-9c) 


One  then  needs  to  factor  this  as  discussed  at  Eq.  (VIII-6d),  which  is  no 
simple  task.  Hence  we  drop  this  example  at  this  point  with  the  comment 
that  a  simple  factorization  to  produce  the  holomorphic  factor  would  be 
most  welcome. 


E.  Symmetrization 

As  we  saw  in  the  last  figure  the  method  may  use  gyrators  where 
actually  none  a rp  apparently  required.  Here  we  show  how  these  gyrators 
can  be  avoided  by  the  procedure  of  Koga  [3,  p.  44]. 

Given  the  PR  admittance  Yc(p),  of  Eq.  (VIII-7d)  for  example,  if 
it  is  not  already  symmetric  we  form  the  following  coupling  admittance 
matrix 


Zs(p) 


11 

2-  12S 

-v  ! 

*i2A  : 

12s 

»  22S 

~*22A 

12A 

^22A 

v 

«  22s  . 

(VTII-lOa) 


where 
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2Zl2S 

=  Zl2 

+  Z21 

2Zl2A 

*V12 

"  ^21 

(Vlll-lOb) 

2%22S 

=  122 

+  J~22 

2^22A 

=  122 

"  122 

(VIII-lOc) 

with  the  subscript  S  and  A  standing  for  the  symmetric  and  (skew) 
antisymmetric  parts.  The  matrix  Yg(p)  is  PR  and  lossless  with 

Y  (p)  ,  and  y  (p,X)  results  at  the  first  n  ports  of  a  circuit 

<•  c  *o 

realization  by  loading  the  final  2b  ports  in  unit  X-plane  capacitors, 

A 

as  we  will  discuss  below.  If  next  we  extract  a  (cascade)  gyrator  from 

each  of  the  final  6.  ports,  as  shov/n  in  Fig.  VIII-5a),  we  obtain  a 

A 

symmetric  coupling  admittance  matrix  Y  ;  for  example,  when  (as  is  the 
normal  situation)  y22q  nonsingular 


vc(P)  = 

-1  - 

Zll+Il2AZ22S^12A 

-1  - 

<wl2S+.*12A~22S^22A 

-1 

Il2A^22S 

(Vlll-lOd) 

-1  ~ 

Zl2S+Z22A£22s£l2A 

-1  ~ 

J22S 22A^  22S^22A 

-1 

Z22\l22S 

-1  ~ 

^22sil2A 

-1  ~ 

^22S^22A 

-1 

^22S 

The  extracted  gyrators  can  be  combined  with  the  loading  capacitors  to 
yield  s-plane  inductors  while  Y  ( p)  can  be  synthesized  by  a  reciprocal, 

***c 

passive,  lossless  p-plane  configuration.  The  overall  structure  is  then 
reciprocal  with  y  (p,X)  and  as  shown  in  Fig.  VIII-5b). 


a)  b) 

Procedure  for  Reciprocal  Synthesis 
of  a  Symmetric  y  (p.X) 

mlQ 

Figure  VIII-5 
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To  see  'Aliy  the  met  hoi!  works  lit  ns  reason  as  follows.  Since 

v  (p,>)  is  assumed  symmetric  we  can  write 

in 

=  y]l(p)"  I{Jl2(p)[2'2v(l,)f>icAr2'21(P>  (VIII-lla) 


4  j21(p)[Z22(p)+>i6>.;’  .?]2(p) 


for  which  a  realization  is  seen  to  come  from  the  coupling  admittance 
matrix 


~l(p)  Jll  J2  "12  J~2  -21 


(VIH-llb) 


^2-21  2  22 


0  v 

-  ~22 


That  is,  a  circuit  realization  for  Y  yields  at  the  input  when 

terminated  I>y  2t .  unit  '-plant  capacitors.  Next  we  lind  an  equivalent 
realization  usinp  i  q.  ( V-  2a  1  villi  f  l><  »ii  lv>  I  »  ?  •  1 1  ■  i  .  1  m.i  1  i  on 


T  —  r  1 


]  .  I  . 

****c  >  *-*c.  ^ 


(VIII-llc) 


Thus  we  obtain 


Y  .  ( p)  [1  +T  Y,  <|>)l  i  Yf! 

— -  n  -  ,~|  * n  — 


(Vlll-lld) 


which  nixes  I  <1 .  ( \1  I  !-!]:'>.  fin  I'll  propel  ly  ,1-  v  •  1  1  ,1-  I  •  •-.s!"«sn-,s«-'  i  • 

preserved  llirnllu'l  these  operations.  F  1 11a  1  1  \  W"  ‘ ' 1  'l'li’l f  'll '  I  tia  I  jl  yv9 
is  not  non  i  n  11 1  ,•  ;  I'li  t  q ,  (YllT-iod)  it  ean  t>e  11. ut  s',  tiv  .111  or  t  h'  n;oiuil 


t  l  ate;  I  o '  11  ■ 


•  I  :  "  1  •  i  .  1  1 


Tin  nii  i  ions  ■  xainpl  ■  -it  I ' f  1  ,  fYlll-Sdl  : 


I '  1  *  *  ' '  I  •  1 1 1 1  . 


IV  e  tia  V 


1  r*>  /  >  ' 


NOT  RNrxiU 


of  I  -('7-  I  ]e 


p 


(VIII~12a) 


"P 


'A 


Extraction  of  the  gyrator  at  port  three  yields 


Zc<p> 


2t> 

P  +  2 

-p  # 

p  +2 

p  +2 

_ri 

P 

p  +  p  u 

0  2 

p  +2 

p  +2 

(VIII-12b) 


Synthesis  of  Yc<p) ,  which  is  symmetric,  yields  yQ(p,X)  at  the  input 
when  the  second  port  is  loaded  in  a  unit  capacitor  and  the  third  port  in 
a  unit  inductor,  the  latter  two  being  p-plane  short  and  open  circuited 
LC  transmission  lines.  Note  however  that  four  p-plane  (lumped)  reactive 
elements  must  be  used  to  synthesize  Y  (p),  in  contrast  to  the  two  used 
at  Fig.  VIII-4. 


F.  Discussion 

Given  a  nonrational  admittance  matrix  in  p,  Y(p) ,  if  there  exists 
a  x  such  that  Y(p)  =  y(p,X)  is  rational,  PR,  and  lossless  in  the 

M  ftl 

two  variables  p  and  X  =  ctanh(]fp/2) ,  then  the  procedures  of  his 
chapter  can  be  used  to  obtain  a  synthesis.  In  particular  the  synthesis 
uses  both  lumped  and  distributed  LC  components,  a  minimum  number  of 
all  types  when  gyrators  are  also  allowed.  If  the  original  matrix  is 
symmetric  then  also  a  series  of  operations  can  be  used  to  eliminate  the 
gyrators,  but  an  excess  number  of  reactive  elements  is  needed  for  the 
given  procedure,  though  it  seems  that  other  methods  should  be  available 
to  reduce  this  number. 

In  the  treatment  given  we  have  extracted  X-plane  elements  as  the 
load  to  obtain  realization  matrices  which  depend  upon  the  other  variable 
p,  R(p)  =  {A(p),  B(p),  C(p),  D(p)).  Of  course  we  could  have  reversed 
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the  role  of  X  and  p  since  in  y(p,X)  there  is  no  real  preference. 

•»» 

The  only  difference  occurs  in  the  synthesis  where  the  extraction  of  the 
lumped  p-plane  elements  means  that  the  X-plane  eounling  network  needs  to 
be  synthesized  in  terms  of  distributed  elements.  This  latter  though  can 
be  com eiiiently  carried  out  in  terms  of  cascade  synthesis  methods  using 
the  unit-elements  [1,  Chap.  7]  and  is,  thus,  in  some  ways  superior. 

The  same  methods  can  be  used  for  the  synthesis  of  lumped-distributed 
RC  networks  of  considerable  interest  to  the  theory  of  integrated  circuits. 
For  such  one  introduces  a  different  variable  s  =  /p".  T  .-n  a  given  admit¬ 
tance  Y(p)  can  be  synthesized  by  a  synthesis  of  the  lossless  admittance 

[5] 

y(s,X)  =  4=  Y(p)  (VIII-13) 

-  /P  ~ 

Such  a  synthesis  can  follow  that  of  the  text  with  the  s-planc  elements 
replaced  by  resistors  (for  the  inductors)  and  capacitors  while  the  X- 
plane  elements  are  replaced  by  RC  lines  to  obtain  the  original  p-plane 
Y (p)  . 

In  the  case  where  there  are  nonrational ly  related  lines  the  methods 
discussed  can  be  extended  by  considering  v-variable  matrices,  with  v  >  2. 
Although  minimal  realizations  can  relatively  easily  be  given,  as  yet  it 
has  not  been  possible  to  obtain  a  PR  coupling  admittance  in  terms  of 
v— 1  of  the  variables. 
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H.  Exercises 


1.  Synthesize  the  lossless  PR 


a)  y(p,X)  = 


2\(2p^+l) 

4Xp+2p2+l 


b)  y(p,X)  = 


4Xp+2pz+l 

2X(2p2+l) 


2.  Prove  that  Y  (p)  of  Eq.  (VIII-7d)  is  PR  and  lossless. 

3.  Carry  out  the  steps  for  p-plane,  instead  of  X-plane,  extractions. 

4.  Show  that  the  gyrator  extraction  of  Fig.  VIII-5a)  yields  Y^p)  of 
Eq.  (Vlll-lOd) .  Carry  out  the  details  when  yoori  is  singular. 

n  ZZb 

5.  Obtain  a  realization  for  y  (p,X)  using  the  method  of  Chapter  IV 

t*»o 

and  show  how  to  obtain  the  realization  of  this  chapter  from  the  other. 

6.  Analyze  any  2-port  lossless  lumped-distributed  circuit  and  from  the 
resulting  y(p,X)  synthesize  the  network  by  the  methods  of  this 
chapter.  Compare  the  final  circuit  with  the  original  and  discuss 
the  problems  raised. 
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II  y  a  ainsi  une  part  de  la  vie,  —  et 
c’est  la  meilleure,  la  plus  pure  et  la 
plus  grande,  --  qul  ne  se  mele  pas 
a  la  vie  ordinaire,  et  les  yeux  des 
amants  eux-memes  ne  percent 
presque  jamais  cette  digue  de 
silence  et  d' amour. 

M.  Maeterlinck 

"Les  Avertis"  du  "Trdsor  des  Humbles" 


CHAPTER  IX 

TIME-VARIABLE  SYNTHESIS 


A.  Summary 

Using  similar  but  generally  somewhat  different  techniques  than  for 
time- Invariant  structures,  realizations  for  impulse  responses  can  be 
obtained  and  manipulated  to  yield  synthesis  results.  Of  primary  interest 
is  that  discussed  for  transfer  voltage  functions  and  that  for  special 
types  of  admittances. 

B.  Properties  of  Impulse  Response  Matrices 

If  we  are  given  the  state-variable  equations  with  time  variable 
coefficients 


i(t) 

=  A(t)s(t)  +B(t)u(t) 

(IX-la) 

y(t) 

m 

=  C(t)s(t)  +  D(t)u(t) 

(IX-lb) 

we  can  find  the  zero  state 

output  through 

00 


y(t)  =  /  [D(t)6(t-T)4C(t)<5(t,T)B(T)]u(T)dT  (IX-lC) 

-00 

where  <t>( t , t)  is  the  state  transition  matrix  satisfying 


124 


SEIr-67-110 


d^t,T)  =  A(t)$(t,T),  t>  T 


(IX- Id) 


$(t,t)  =  lk,  $(t,x)  ----  Ok  t  <  t 


(IX-le> 


In  actual  fact  since  <f  satisfies  the  differential  equation  of  Eq. 

tart 

(IX-ld)  it  can  be  shown  [1,  p.530]  to  be  the  product  of  two  matiices, 
one  in  t  and  one  in  t 


$(t,x)  =  ~(t)A(x)l(t-x) 


(IX-lf) 


where  l(t-x)  is  the  unit  step  function.  Further  the  number  of  rows  in 
A(t)  can  be  assumed  reduced  to  its  minimal  size  6,  this  size  being 
independent  of  t  for  reasonably  behaved  A(t)  [1,  p.530] 

As  a  consequence,  we  can  associate  with  the  state- variable  equations 
an  impulse  response  matrix  [here  6(t)  is  the  unit  impulse] 


T(t,T)  =  D(t)6(t-T)4C(t)H(t)A(T)B(T)l(t-T) 

m  •»  m  —  **  «•» 


(IX-2a) 


such  that 


y(t)  =  f  T(t,x)u(x)dx 

—  J  — 


(IX-2b) 


This  latter  can  be  conveniently  denoted  as 


y  =  T*u 

M 


(IX-2c) 


Since  T  contains  impulses  it  i3  often  referred  to  as  a  distributional 
kernel  defining  the  mapping  of  u  into  y,  y  =  T»u.  If  we  have  two 
such  mappings  defined  by  kernels  T  and  T„  we  can  apply  one  after 
another,  as  might  occur  in  a  cascade  of  voltage  transfer  functions.  This 
leads  to  the  definition  of  the  composition  T  oT  through 


l  - 


(IX-2d) 
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As  an  integral  this  conposition  takes  the  form 


Tj oJ2 ( t , t )  =  j  Ijft^JTjtX.TjdX  (IX-2e) 

.00 

Through  the  concept  of  composition  the  inverse  of  a  kxk  kernel  can  be 
defined  by 

T_1oT  =  ToT"1  =  &(t-T)l,  (IX-2f) 

mm  mm  — K 

Consequently  T  can  be  given  the  representation  alternate  to  Eq.  (IX-2a) 

** 

as 

T  =  D6+(C&)o[6,lk-A6]“1o(B&)  (IX-3a) 

Since  [&’  (t-T) ]  1  =  l(t--r)  we  see  by  comparison  with  this  last  expression 
that  if  we  are  given 

T(t  ,t)  =  H(t)6(t-T)+'r(t)0(T)l(t-T)  (IX-3b) 

■—  m  •* 

then  a  possible  realization  is 

A  =  0,  B  =  9,  C  =  T,  D  =  H  (lX-3c) 

This  is  minimal  if  the  number  cf  rows  in  6  has  been  minimized.. 

«•* 

If  we  make  a  transformation  on  the  state 

s  ( t )  =  ,T(t)s(t)  (IX- 4a) 

**  •» 

then,  since  the  transforming  matrix  must  now  be  differentiated  we  have 

&  =  B  =  flB,  C  =  Cff"1  (IX- 4b) 

Consequently,  the  freedom  of  using  time- variable  transformations  allows 
one  to  change  the  structure  of  the  A  matrix,  resulting  in  some  rather 
interesting  phenomena. 
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.  Passive  Voltage  Transfer  Function  Synthesis 

Let  us  consider  the  problem  of  synthesis  of  kernels  mapping  voltages 
into  voltages;  the  material  follows  to  a  large  extent  the  ideas  of 
Silverman  f2]. 

As  a  preliminary,  let  us  first  observe  that  if  we  define,  for  a 
given  £(t)  a*>d  a  fixed  tQ, 

t 

V(t)  =  /*(t,T)*(t,T)dT,  t>  t  (IX-5a) 

*0 

(which  is  positive  definite)  then  the  choice 

-1  1/2 

?  =  (V  )  (IX-5b) 

■"*  mm 

yields  on  using  Eqs.  (IX-ld,  4b). 

A  +  A  =  -V-1  (IX- 5  c) 

w*  X m  m* 

As  a  consequence,  from  what  we  previously  learned  at  Eq.  (V1I-11)  we 
should  be  able  to  use  this  transformation  for  a  passive  synthesis.  We 
comment,  however,  that  V(t)  varies  with  time  even  in  the  time-invariant 
case  so  that  slightly  different  procedures  are  preferable  when  A,  B,  C 

Ml  fm  Wm 

are  constant. 

As  the  next  preliminary  let  us  synthesize  a  voltage  to  current 
transfer  function  (kernel)  T,  i '  =  T«v, ,  where  i '  and  v,  are 
measured  at  different  ports.  Given  any  realization,  say  the  one  of 
Eq.  (IX-3c)  let  us  perform  the  transformation  of  Eq.  (IX-5b)  to  obtain 

§  =  AS  +  Bv  (IX-6a) 

to  to  to  to  to  J 

i '  =  es  +  Dv,  (IX-6b) 

to  2  to  to  to  toj_ 

Let  us  next  introduce  another  set  of  variables,  the  current  i. 

’  —l 

associated  with  the  first  ports  and  v^  the  voltage  associated  with  the 
final  ports  to  write 
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(IX-6c) 


s  =  As  +  [6,6]  v 

m  i%  m  fc  MiJ, 

.Si. 

'  i,  1  =  ils  +ro  -Dir  V,  (IX-6d) 

i '  6  D  0  v! 

-2  L-  *“  ~2 

Note  that  if  we  sm  (I  and  ignore  t  he  inpul  port  currents: 

then  the  origin  .1  description  is  returned.  However,  as  in  the  time- 
invariant  case,  Eqs.  (IX-6c,d)  define  a  coupling  (time- variable)  resistive 
network  through 


zc<t,T)  =  r  2 

-|(t)  b(t-T) 

(IX-6e> 

£(t> 

Q 

■* 

-C(t) 

.  B(t) 

C(t) 

-A(t) . 

Note  that,  by  virtue  of  Eq.  (IX-5c) 

Y  +  Y  s  0  +  V_1(t)&(t-T)  (IX-6f ) 

*•  C  •“*  C  *•* 

in  which  case  Yc  can  be  synthesized  by  time-variable  gyrators  and 
resistors  both  of  which  are  passive.  Termination  of  the  resultant  net¬ 
work  in  unit  capacitors  yields  Eq.  (IX-6c,d).  At  the  final  ports  we 
can  next  insert  unit  gyrators  to  obtain 

=  ~2’  h  =  2a  (IX-7a) 

Setting  v'  =  0  results  in  an  open  circuit  load  while  i'  =  v„  yields 

v  =  T» v  (IX-7b) 

•+  •"I 

As  a  consequence  the  procedure  results,  for  t  >  t  ,  in  a  passive 

realization  of  any  T(t,t)  of  the  form  of  Eq.  (IX-3b).  Since  practically 

such  constructs  are  only  used  after  a  finite  time,  the  t  >  t  restriction 

o 

is  of  no  practical  restriction;  but  in  some  cases  t  =  -°°  can  be  used 
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In  which  case  the  theory  of  Silverman  results  when  D  =  0.  The  synthesis 
is  summarized  in  Fig.  IX-1 . 
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M  C  t ) 


1 


(IX-8e> 


Circuit  for  V2  =  (3p/[p+2])V1 
Figure  IX-2 

are  time- variable  even  though  the  overall  network  is  time-invariant,  c) 

the  elements  become  unbounded  for  t  approaching  t  .  If  we  would  have 

o 

chosen  A  =  -2  and  t  =  -»  this  latter  (unboundedness)  could  have 
been  avoided  while  n  slightly  different  approach  (see  the  Exercises) 
would  allfw  n  time- in  variant  synthesis. 
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I 


.  Passive  Admittance  Synthesis 


Following  the  previous  ideas  we  can  form  the  ccupling  admittance 


matrix 


Y  (t,T)  = 

»c 


D(t) 

J(t)B(t) 


-C(t)J(t) 

■a 

-1  •  -1 

■J(t)A(t)J(t)-J(t)J(t) 


6(t-t)  (IX-9a) 


We  then  wish  for  a  passive  synthesis  to  be  able  to  choose  J  such  that 
the  symmetric  part  of  Y  is  positive  semidef inite.  On  evaluating  tnis 

m  C 

symmetric  part  we  have,  assuming  a  symmetric  J, 


Y  +Y 

»c  »c 


D+D  [B32-C]J-1 

j_1[j2b-c]  -j'Va+a^+Ij2]}?"1 

m  pa  m  I**  as  m  m  *■**  9 


( IX-9b) 


•  o  •  • 

<T  4  —  cnr^crcr 


where  we  have  also  used  Jz  =  Jif+JJ.  In  the  case  where  D+D-  which  is 

m  **•  —  am' 

twice  the  syimnetric  part  of  D,  is  positive  definite  and  A,  B,  C  have 
bounded  entries  the  (Riccati)  equation 


j2a+at2+[j2]  =  -[j2e-c][d+d]-1[bj2-c] 

mm I  ««  ««  «,  m  —  »%  ~  mm  amJ 


(IX-9c) 


is  known  [3]  to  have  a  solution  for  a  nonsingular  symmetric  J.  Con¬ 


sequently, 


Y  +Y  =[  (D+D) 1/2 

» c  »C  mm 


o  j”1(52b-c){d+d)"1/2  1 

m  m  m  m  m  m  m  J  L*»n. 


1  IF1  1  1 

-n  l|_ -n  -nj 


(D+D) 


o  (d+d)'1/2(bj2-c)J_1 

>■*  m.  «a9M  m  m 

(IX-9d) 


which  shows  that  Yc  can  be  synthesized  by  n  constant  resistors 
loading  time-variable  gyrators  (for  the  symmetric  part)  and  time-variable 
gyrators  (for  the  skew-symmetric  part) ;  here  n  is  the  number  of  terminal 
ports. 

We  conclude  that  if  a  given  nxn  admittance  kernel 


y(t,r)  =  D(t)5(t—r)  +  C(t)«(T)l(t-T) 


(IX-10) 
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has  the  symmetric  part  of  I)  positive  definite  (as  well  as  bounded 
minimal  C  and  II  matrices)  then  a  passive  synthesis  can  be  given. 

Such  results  by  solving  the  nonlinear  Riccati  equation,  (IX-9c),  for 
ft'(t) ,  and  forming  Y  which  then  yields  y(t,T)  by  loading  of  the 
passive  coupling  structure  in  unit  capacitors.  Several  observations  are 
worth  noting.  First  M  is  very  difficult  to  obtain,  if  not  impossible 
practically,  since  a  nonlinear  variable  coefficient  differential  equation 
must  be  solved.  Second,  Eq.  (IX-9d)  shows  where  difficulty  arises  if 
the  symmetric  part  of  D  is  singular;  hence  the  method  seems  hard  to 
extend  to  cover  more  general  cases.  Third,  the  presence  of  terms 
E(t)b'(t-t)  is  handled  bv  writing  E(t)L'(t-r)  =  J( t)J(T)0' (t-r)  - 

NM> 

J(t)J(t)6(t-t) ;  if  y  is  known  to  come  from  a  passive  network  this 
decomposition  is  always  possible  since  E  is  then  positive  semidef inite. 
Fourth,  although  the  passivity  conditions  on  _D  (and  E)  are  known, 
those  on  B  and  £  are  not,  except  in  the  lossless  case  where  B  =  C 
is  possible  and  an  alternate  synthesis  applies  to  cover  all  cases  [4] 

{  ;i  =  ln  holds  to  yield  a  skew- symmetric  Yc). 

An  alternate  and  interesting  method  results  from  the  following 
manipulation  [5],  Let 

A  =  liAT-1  +  *:i_1  (IX-lla) 

I**  •»*.*-*  M 

6  =  :ib,  8  =  cm-1  (ix-iib) 


then  from  Eq.  (IX-9c) 


-:r1r:i2D-cirD+51“1rffi'2-c1:t' 

***-'  4»- 

/v\ 

2LL 


(IX-llc) 


where  £  is  defined  as 


1 


(D+D)-1/2re:,2-cl:.-1 

#**  **■*  M  I**'  ***** 


If  further  we  define 


(IX-lld) 
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z 

=  (£^,1/2  -  i 
v  2  ’ 

(IX-lle) 

where  the  positive 

definite 

symmetric  square  root 

is  again  meant,  we 

obtain 

s  = 

-i  ( X-X)  s + (6+hz )  v-L  v  * 

2.  a*  tm  m  m  a*  ••  M 

(IX-12a) 

i  = 

am 

(%+Lz)  s+zv* 

m  ««*  mm 

(IX-12b) 

i*  =  ■ 

-V*  =  -ts-zv 

( IX- 12c) 

Here  direct  substitution  of  the  last  constraint  upon  noticing  that 
2zL  =  yields  the  original  set  of  equations 

m  m  m 

s  =  Xs+Bv,  i  =  &s+[»i£]v  (IX-12d) 

a*  ^  *  m  W  w>  ■*  Cd  !•» 

The  constraint  i *  =  -v*  corresponds  to  resistive  loads  at  the  v*.  i* 
ports.  As  a  consequence  we  consider  the  coupling  admittance  matrix 


«■»  c 


o  z  -8+z£ 

m  **  m  m  «***- 

-z  0  L 

MW  **4 

A  A"  A  I  A  A 

B+Lz  -L  -(A-A) 

M  m«“\  »■»  2 


(IX-13) 


which  is  skew-symmetric  and  hence  realizable  by  gyrators.  When  loaded 
at  the  final  ports  by  unit  resistors  and  at  the  next  to  final  ports  by 
unit  capacitors, the  input  admittance  y(t,t)  occurs  at  the  input  ports. 
In  this  manner  an  alternate  synthesis  results  when  D+D  is  nonsingular, 
for  a  passive  y(t,f).  It  should  be  observed  that  this  method  requires 
that  the  skew-symmetric  part  of  D  must  be  extracted  before  Eqs.  (IX-12) 
are  considered,  as  seen  from  Eq.  (IX-12d).  Of  course  the  skew-symmetric 
part  of  D  is  obtained  by  gyrators  connected  in  parallel  with  the  input 
ports.  Note  that  this  again  shows  that  all  time- variations  for  time- 
variable  circuits  can  be  placed  in  the  gyrators. 
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E.  Discussion 


Because  the  state- variable  equations  are  expressed  in  the  time 
domain  they  are  primarily  suited  for  obtaining  syntheses  of  time-variable 
networks.  Here  we  have  investigated  two  types  of  synthesis,  one  for 
voltage  transfer  and  the  other  for  n-port  admittance  impulse  responses. 

Emphasis  has  been  placed  upon  passive  structures  but  it  is  clear 
that  the  same  ideas  can  be  applied  to  synthesis  using  active  elements, 
perhaps  in  an  even  simpler  manner.  The  transfer  function  synthesis 
contains  relatively  simple  calculations  while  the  solution  of  a  nonlinear 
differential  equation  makes  the  admittance  syntheses  extremely  difficult 
to  carry  out.  Consequently  one  would  hope  for  a  simpler  admittance 
synthesis  and  in  fact  one  which  relaxes  the  unnecessary  constraint  of  a 
nonsingular  symmetric  part  for  D. 

In  the  time-invariant  case  the  methods  yield,  in  general,  circuits 
with  time-variable  components.  In  some  instances  these  can  be  combined 
to  obtain  time-invariant  components  but  the  result  does  show  that  perhaps 
some  other  synthesis  methods  exist  which  reduce  to  the  known  time- 
invariant  techniques  perviously  discussed.  It  is  worth  observing  though 
that  many  of  the  previous  concepts  discussed  only  for  time-invariant 
structures  do  extend  to  the  time-variable  situation.  For  example  it 
seems  relatively  simple  to  set  up  a  theory  of  equivalence  for  time- 
variable-  structures  from  the  discussions  in  Chapter  V. 

Although  the  n-port  synthesis  techniques  have  been  given  in  terms 
of  admittances  the  classical  synthesis  methods  in  terms  of  scattering 
matrices  can  also  be  extended  [6]  [7]  [8]  though  as  yet  these  latter 
time-variable  methods  have  not  really  applied  the  concepts  of  state- 
variable  theory  for  their  success. 

As  a  point  of  philosophical  interest  we  point  out  thfit  the  passive 
synthesis  of  Section  C  can  be  applied  to  non-stable  network  functions, 
such  as  T  =  l/(p-l).  Consequently  one  can  relatively  easily  construct 
passive  unstable  networks,  a  rather  paradoxical  situation  when  it  is 
realized  that  many  intuitive  deductions  concerning  passive  networks  have 
rested  upon  the  "stability"  of  passive  structures. 
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G.  Exercises 
1 


Synthesize  by  a  passive  structure  the  voltage  transfer  functions 

p+2 


a)  T(p)  = 


(P+3) (p+1) 


b)  T(t,  r)  -  tT2e'(t-T)l(t-r) 


2.  Synthesize  by  the  methods  described  the  time-invariant  admittance 


y(p) 


2p_ 

p+3 


From  the  result  discuss  various  simplifications  which  can  be  made, 
or  need  to  be  made,  in  the  theory. 


* 
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3.  Give  a  synthesis  for  time-invariant  voltage  transfer  functions  using 
ideas  similar  to  those  of  Section  C.  For  this  one  can  choose  any 
positive  definite  constant  matrix  V  and  solve  for  J  to  yield 

Am 

Eq.  (IX-5c) . 

’  4.  Synthesize  the  voltage  transfer  function 

T(p)  =  FI 

by  the  method  of  Section  C.  From  the  result  discuss  why  a  passive 
network  need  not  be  stable. 

5.  Discuss  means  of  solving  Eq.  (IX-9c)  for  fF. 

*6.  Develop  a  state-variable  synthesis  of  passive  scattering  matrices. 

7.  Extend  the  results  of  this  chapter  to  nonpassive  structures  and 
discuss  the  meaning  of  your  methods. 
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Nous  vivons  a  cote  de  notre  veritable  vie 
et  nous  sentons  que  nos  pens^es  les 
plus  intimes  et  les  plus  profondes 
meme  ne  nous  regardent  pas,  car  nous 
sommes  autre  chose  que  nos  pens£es 
et  que  nos  reves.  Et  ce  n'est  qu'a  certains 
moments  et  presque  par  distraction  que 
nous  vivons  nous-memes. 

M.  Maeterlinck 

"Les  Avertis"  du  "Tr^sor  des  Humbles" 


CONCLUSIONS 

Paradoxically  the  simple  expediant  of  introducing  a  set  of  first 
order  differential  equations  to  describe  high  order  ones  has  led  to  the 
solution  of  previously  unsolved  problems,  such  as  the  determination  of 
all  equivalent  active  structures  for  a  given  network.  As  we  have  seen 
there  are  many  areas  where  the  ideas  can  be  applied,  perhaps  with  a 
possibility  of  gaining  insight  into  the  behavior  of  a  system. 

Thus,  because  most  systems  of  practical  significance  possess  an 
identifiable  state,  the  state-variable  equations  give  a  general,  or 
universal,  means  of  observing  systems.  By  keeping  track  of  the 
solutions  of  the  describing  equations  in  state-variable  form  one  can 
keep  track  of  the  behavior  of  the  subparts  of  a  system  in  orderly  fashion. 
And  because  this  tracking  can  be  done  orderly,  the  theory  allows  readily 
for  the  computer  analysis  of  networks,  this  analysis  having  the  possibil¬ 
ity  of  proceeding  in  two  ways,  as  we  have  seen  in  Chapter  II  in  either 
the  topological  or  capacitor  extraction  form.  Once  a  computer  analysis 
is  set  up  in  this  manner  the  results  can  be  used  for  other  purposes 
than  keeping  a  record  of  voltages  and  currents;  for  example  Chapter  VI 
has  shown  how  a  sensitivity  analysis  can  proceed  from  a  state-variable 
analysis  program. 

But  the  most  striking  uses  of  the  theory  occur  when  synthesis  is 
considered.  Here  we  have  seen  that  minimal  degree  realizations,  that  is 
minimum  reactive  element  circuits,  result  for  general  transfer  functions 
by  the  theory  of  Chapter  IV.  Even  though  this  latter  is  somewhat 
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abstract  its  significance  should  not  be  overlooked.  Because  of  its 
form  it  allows  convenient  integrated  circuit  constructions  as  well  as 
analog  modeling  for  simulation  and  preliminary  testing  of  designs.  Also 
because  of  its  algebraic  form  the  realization  technique  allows  for  the 
complete  computer  design  of  a  system,  though  as  yet  such  a  program 
remains  to  be  carried  out.  In  the  area  of  classical  multiport  synthesis. 
Chapter  VII  has  shown  that  the  introduction  of  state-variables  can  lead 
to  a  contribution  since  a  minimal  resistor  and  minimal  capacitor  circuit 
results  by  application  of  the  given  method. 

Still  it  is  by  way  of  generalization  of  the  positive-real  admittance 
synthesis  where  the  most  significant  contributions  of  state-variable 
theory  seem  to  be  made.  We  have  illustrated  this  in  two  different  ways. 
The  first  is  through  the  introduction  of  a  second  variable  to  allow  for 
design  with  both  lumped  and  distributed  elements,  as  covered  in  Chapter 
VIII.  The  second  generalization  is  that  of  Chapter  IX  for  the  synthesis 
of  time- variable  circuits.  Though  this  latter  is  as  yet  not  completely 
finished,  to  us  it  represents  a  beautiful  application  of  the  theory  which 
in  almost  all  parts  is  carried  out  in  the  time  domain. 

Once  a  circuit  has  been  designed  the  material  of  Chapter  V  on 
equivalence  shows  how  many  other  circuits,  in  fact  almost  all,  with  the 
same  terminal  behavior  can  be  found.  To  complete  the  picture  any  of 
these  can  be,  in  turn,  analyzed  by  the  methods  of  Chapters  I  and  II  to 
check  its  performance. 

In  summary,  the  theory  of  state-variables  has  allowed  an  almost 
complete  picture  of  the  theory  of  networks,  in  fact  within  the  larger 
framework  of  scientific  systems.  It  has,  however,  raised  many  fascinating 
problems,  some  of  which  we  have  tried  to  point  out  along  the  way.  Thus, 
though  the  theory  may  offer  little  to  some  people  it  can  offer  an  immense 
amount  to  those  who  would  allow  it  —  so  is  it  with  almost  all  that  we 
rneai;. 


Quel  jour  deviendrons-nous  ce  que  nous  sommes? 

Nous  nous  ecartions  sans  rien  dire  et  nous 
comprenions  tout  sans  rien  savior, 

M.  Maeterlinck 

"Les  Avertis"  du  "Tr^sor  du  Humbles" 
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LIST  OF  PRINCIPAL  SYMBOLS 
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•» 

b,  number  of  branches . ....26 

input-state  matrix . ..7 

C,  state-output  matrix . ...7 

c,  capacitance  matrix . 21 

cut  set  matrix . ...27 

D,  canonical  equation  matrix... 7 

E,  coefficient  matrix . .7 

g(p),  least  common 

denominator . 30 

h,  impulse  response . 8 

H,  transfer  function . 9 

i,  port  currents . 27 

£[  ]>  Laplace  transform.. . 9 

J,  number  of  links . 26 

m,  input  size . 4 

n,  output  size . 4 

p  a  o+Ju,  frequency . 9 

P,  observability  matrix. ..... .64 

PR  lemma  matrix . 96 

Q,  quality  factor . ,48 

Q,  controllability  matrix . 64 

R,  realization . 7 

resistance . 6 


s,  number  of  separate  ports... 26 

Pk 

8^  ,  pole  position 


sensitivity . .86 

s,  state . 4 

|M  * 

sensitivity . 84 

S  ,  Hankel  matrix . 62 

t»r 

t,  number  of  tree  branches ....  26 

time . 3 


t  ,  initial  time . 3 

o’ 

T,  state  transformation . 16 

T[,],  system  transformation . 4 

T  .  modified  Hankel  matrix . 115 

*r 

‘J ,  tie  set  matrix . 28 

state  transformation . 126 

u,  input . 3 

U,  factor  of  T  . 115 

U~  ,  left  inverse  of  U . 115 

v,  number  of  vertices . 26 

v,  port  voltages . 27 

V,  transition  energy  form . 127 

»» 

W(p),  para-Hermitian  factor . 97 

W  ,  right  inverse  of  W . 97 

"r 

y,  output . 3 

y(t,T>,  admittance  kernel . 131 

m 

Y c,  coupling  admittance. ,, .20,  128 
Y_,  symmetric  coupling 

v 

admittance . 119 

Y  ,  gyrator  coupling  matrix . 51 

5,  impulse,  =  6[  ] . 9,  10 

6.,  X-degree  =  rank  T....115,  116 

X,  frequency  variable . Ill 

9,  zero  state . 4 

undamped  natural  frequency. .48 
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0,  zero  matrix . 4 
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l(t),  unit  step . 9 

1  ,  kxk  identity . 9 
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*p,r 

+  ,  direct  sum . 51 

transpose . 23 

•  ,o,  compositions . 125 
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